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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. ‘The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 
ty using double spacing and leaving a margin of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 
margin. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it. will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 1g, 239, 1938. 

H. Jeffreys, Theory of Probability, and edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation,—-For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (7yrans. 
1A.U., Vol. VI, p. 345, 1938). Council has decided to adopt the I1.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, p. 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London: Royal Society, 1951) except where these conflict with 1.A.U. practice. 


4 —— These should be designed to appear upright on the page, drawn 
about twice the size required in print and prepared for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should appear. Blocks are retained by the Society for 10 years; unless the author 
requires them before the end of this period they are then destroyed. 


5. Tables.-These should be arranged so that they can be printed upright on 
the page. 


.6, Proofs.--Costs of alteration exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society's interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts.—-When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MSS. 


Reading of Papers at Meetings 
8. When mpeg _ authors are requested to indicate whether they will be 
t 


willing and able to re € pe er at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


9. Postcards giving the programme of each meeting are issued some days before the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF t955 MAY 13 


Professor Sir Harold Jeffreys, President, in the Chair 


‘The election by the Council of the following Fellows was duly confirmed : 

Vivian Mary Bradley, B.Sc., Department of Astronomy, University of 
Manchester (proposed by Z. Kopal) ; 

Gilbert Fielder, B.Sc., Department of Astronomy, University of Manchester 
(proposed by Z. Kopal) ; 

Vera Hewison, B.Sc., Department of Astronomy, University of Manchester 
(proposed by Z. Kopal) ; 

George Edwin Stone Kendall, Ph. C., F.N.A.O., M.P.S., 23 High Street, 

Ascot, Berkshire (proposed by R. L. Waterfield) ; 


Peter Yorke Millns, B.Sc., Department of Astronomy, University of 
Manchester (proposed by Z. Kopal) ; and 


Gordon Ernest ‘Taylor, Celestia, London Road, Horsebridge, near Hailsham, 
Sussex (proposed by D. H. Sadler). 


The election by the Council of the following Junior Members was duly 
confirmed ; 


Robert Owen Bishop, University Observatory, Oxford (proposed by 
H. H. Plaskett) ; and 


Francis Wykes, 1 Brookmead, Hildenborough, Kent (proposed by P. Moore) 


One hundred and eighteen presents were announced as having been received 
since the last Meeting, including : 
Sir Richard Phillips, Wonders of the Telescope (presented by Sir Eric Miller); 
Roger Boscovitch, Nova methodus adhibendi phasium observationes (presented 
by Sir Eric Miller) ; 
J. Cassini, Tables astronomiques du Soleil, de la Lune, des Planétes, etc 
(presented by Sir Eric Miller) ; and 
International Astronomical Union, Symposium No. 1.—Co-ordination of 
galactic research (presented by the Cambridge University Press). 
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THE SYMMETRY OF THE CORONA OF 1954 JUNE 30 
T. Gold 


(Communicated by the Astronomer Royal) 


(Received 1955 March 2) 


The corona seen during the total eclipse of 1954 June 30 was extraordinary 
for its regularity. ‘che time of the eclipse was as nearly the sunspot minimum 
as can be determined, and no spot was seen at the Solar Observatory at 
Herstmonceux for 8 days preceding or 12 days following the eclipse. The 
corona did not possess any streamers extending out in odd directions as is normally 
the case, but was of a high degree of symmetry. In the equatorial direction 
there was the furthest extension, a little more to the east than to the west. In 
higher latitudes there were polar plumes fairly regularly spaced. ‘There was 
a sharply defined position of changeover on the Sun’s limb in each quadrant 
between the striated polar plume zone and the smooth zone associated with the 
equatorial extension. 

The Royal Greenwich Observatory expedition to Syd Koster island, Sweden, 
(T. Gold, A. Hunter, E. W. Bastin, D. W. Sciama) had intended to photograph 
the eclipse star field with a telescope of 21 ft focal length, for the purpose of a 
measurement of the Einstein deflection of light. Thin cloud prevented this, 
but a good picture of the corona was obtained, marred by cloud only in the north- 
west quadrant. ‘I‘he negative contains a very large range of density and this 
renders it unsuitable for direct reproduction in print. ‘The print shown was 
obtained from the negative by a process involving the subtraction of intensity 
by an amount which was a suitable function of radius only. The real isophotes 
are much more nearly circular than appears in the reproduction; but the 
reproduction serves to show some of the coronal pattern. On the negative the 
following measurements were carried out. 

1. The direction of the equatorial extension.—A line was drawn through the 
equatorial extension splitting it as symmetrically as possible into the northern 
and southern parts. On the east side this required drawing it through the 
furthest extension that could be seen on the plate (about 4 solar radii from the 
limb) while on the west side the extension was clearly and quite symmetrically 
bifurcated, and the line had to be drawn through the middle of that. This line, 
defined solely by the outer corona, passed accurately through the centre of the Sun. 
Its direction seemed indeterminate by not more than +4 degree. The angle 
that this line made with the normal to the projection of the Sun’s axis was o-8 
degrees. 

2. The changeover between the polar and the equatorial type of corona.—The 
lines separating the striated polar regions from the smooth equatorial ones can 
be seen right in tothe Sun’slimb. The heliographic latitudes of these changeover 
points on the four quadrants of the limb are very similar and do not appear to be 
uncertain to more than about one degree. ‘They are: 

N.E, 60° 
N.W. 59}° 
5.W. 60° 
S.E. 634° 
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The symmetry of shape of the corona was considerably more marked between 
north and south than between east and west, both so far as the inner features are 
concerned, and also the furthest equatorial extensions. ‘This symmetry is known 
to be absent at other phases of the sunspot cycle, when disturbances on the surface 
of the Sun appear to influence or even to determine the coronal appearance. 
The corona at sunspot minimum would therefore appear to have the special 
interest that it reveals the shape that the system settles down to when not subjected 
to the erratic disturbances marking periods that are richer in sunspots. 

The implication of this high degree of symmetry would seem to be that the 
equatorial regions of the corona, even out to a distance of 4 radii, are very much 
under the control of the Sun. ‘The relaxation time in which the corona settles 
down to this symmetrical condition is not likely to be more than a few months, 
judging from the evidence of previous near-minimum eclipses which still showed 
great disturbances. ‘The mechanism which transfers this precise information 
about the rotation of the Sun out to that distance is not clear, but it seems likely 
that any explanation will require an outward flow of material in the equatorial 
region. Any suggested pattern of flow involving an inward flow in the equatorial 
zone would greatly augment the difficulties, for it would then be necessary to 
find a mechanism capable of pulling the incoming material into a shape determined 
by the Sun, while it is still at such a great distance. 

The polar plumes, on the other hand, do not require an outward flow for theit 
explanation. Most patchy distributions of material at some distance from the 
Sun would be drawn into such shapes by falling under the combined effect of 
gravitation and a magnetic field, whose lines of force would thereby be made 
visible. If, as seems very likely, the polar plumes are indicators of magnetic 
lines of force, then another consideration is of interest. ‘The shapes of the 
plumes form an ordered sequence (as would be expected from lines of force 
resulting from electric currents in the Sun) and the boundary of the plume 
zone, namely an edge of the equatorial extension, is clearly an adjacent shape in 
that sequence. If the plumes are the lines of force, then the outlines of the 
equatorial extension, at any rate the parts within the range where the plumes 
can be seen, must also be lines of force. In searching for the physical difference 
that underlies the clear distinction in the appearance of the two zones we ar 
again led to a hydrodynamical explanation, not only because of the difficulty of 
invoking another physical effect, but also because any streaming has to occur along 
lines of force in a system such as this, where ohmic dissipation plays no essential 


part, and a changeover in the direction of flow must therefore be on a surface 
defined by lines of force. 

The greatest departure from the shape of a dipole field that can be seen in 
any of these outlines is in the edge of the equatorial extension and the nearby 


plumes. While there would be no need to suppose the currents in the Sun to 
be just those that make an external dipole-like field, the departure in shape is in 
the sense of making the loops of the lines of force stretch out very much further 
than a solar dipole field would. ‘This is very unlikely to result from currents 
in the Sun (where sets of nearby opposing current rings would be required) 
but would result from the outward convection of the field due to an outward flow 
in the equatorial zone. 

Estimates of the surface field of the Sun and of the conductivity of the corona 
at different ranges from the Sun will enable an estimate to be made of the strength 


24° 


‘ 
‘ 


342 T. Gold Vol. 115 


of the magnetic forces that the coronal material has to bear to uphold this 
deformation of the field; and with a knowledge of the densities this would 
specify the order of magnitude of the streaming velocities. 

I should like to acknowledge the skilled work of Dr Hunter, and the valuable 
assistance of the volunteer members of the expedition, Dr E. W. Bastin and 
Dr D. W. Sciama. 


Royal Greenwich Observatory, 
Herstmonceux 
1955 February. 
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INTERFEROMETRIC MEASUREMENTS OF THE GREEN CORONA 
LINE DURING THE TOTAL SOLAR ECLIPSE OF t954 JUNE 30 


A. H. Jarrett and H. von Kliiber 
(Received 1955 March 1s) 


Summary 


A detailed description is given of a successful attempt to observe inter- 
ference fringes with a Fabry-Perot interferometer from the 5303 A corona 
line during a total eclipse. 

Low-absorption coatings for the interferometer and a narrow pass-band 
interference filter (equivalent width of ~25 A of the continuous spectrum, 
peak transmission 25 per cent) have been employed. 

In spite of thin cloud during totality, which was estimated to have reduced 
the available intensity by a factor of 2, it proved possible to measure fringe 
half-widths for various position angles and distances from the solar limb 

The measurements seem to indicate a decrease of the line widths with 
increasing distance from the Sun. Interpreted in terms of kinetic tempera- 
tures, T'=2°2 10% to 5 10° deg. K, with most values near 2°5 10°, for 
regions situated 1°05 to 1°3 radii from the solar centre 

The paper concludes with suggestions for improvements to the equipment 
for a future occasion 


The possibility of getting interference fringes of the green corona line with 
a Fabry-Perot interferometer during an eclipse, by using the well-known method 
that Buisson, Fabry and Bourget (1) applied at one time to the Orion nebula, 
was considered and tried as long ago as 1918 (2), but so far has always failed for 
one reason or another. W. P. Wright, and later H. D. Curtis (2), developed and 
tried specially constructed instruments for this purpose, at no less than four 
eclipses. At two of these eclipses (1926 and 1929) one of us had an opportunity 
of becoming acquainted with these experiments. Later J. A. Carroll (3) made 
three further attempts on these lines, His last attempt, in 1947, was frustrated 
by damage to the interferometer plates in a serious aeroplane accident. 
A. A. Kaliniak (4), with a method differing in principle from the one considered 
here, using a slit spectrograph combined with a Fabry-Perot interferometer, 
obtained interference fringes with the green coronal line at the eclipse of 1941, 
but only with a low resolving power of about 5000. In recent years we have 
learnt more about the width and intensity of the green line, while considerable 
technical progress has been made in preparing narrow-band light filters and 
efficient coatings for interferometer plates. In a short report, P. J. ‘Treanor (§) 
recently pointed out the conditions under which it should be possible, with 
modern optical equipment, to get interference fringes of the desired type during 
aneclipse. ‘This report, together with the fact that one of us (6, 7, 8) is acquainted 
with the technique of producing low-absorption interferometer coatings, were 
important reasons why we were led to make another attempt on this problem. 
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‘The aim and purpose of this present experiment was not, as the application 
of an interferometer may suggest, to achieve very high resolving power for 
investigation of the green line. ‘This, if desired, could be done better with 
optical arrangements such as those used successfully for this purpose by Lyot 
and his collaborators (9, 10, 11, 12) without an eclipse. Because of the large 
width of the green line a resolving power of ~50 000 (~ 0:1 A) for most purposes 
would be quite sufficient. Since the line is rather broad, even less powerful 
instruments than those used by the authors just mentioned can give quite good 
results, as was shown by Waldmeier (13, 14, 15). 

Were it possible, however, to photograph the corona during an eclipse with 
a Fabry-Perot interferometer (as was done in (1) for the Orion nebula), and in 
this way to obtain superimposed over the whole image of the corona a set of 
interference fringes of the green line, with one exposure only, the result might 
provide valuable information on the behaviour of the width of the green line in 
different parts of the corona, and should perhaps also give information with regard 
to Doppler shifts. 

As the 1954 eclipse took place almost exactly at a very low minimum of solar 
activity, the intensity of the green line was expected to be very low, and the 
successful outcome of our experiment was, therefore, considered to be rather 
doubtful. ‘To avoid as far as possible the loss of money and time which would 
arise from a failure, no elaborate equipment was used and the experiment was 
arranged rather as a first trial with future eclipses in mind. Although very 
thorough preliminary tests at Cambridge were essential, the optical arrangement 
as actually used during the eclipse was very simple. 

In view of the history of this type of eclipse work one of our aims in this paper 
has been to describe rather completely our experimental experience, as we 
think it should be useful to others wishing to pursue further work of this kind. 

Optical arrangement.—-'The optical requirements of the apparatus fall under 
two headings: (a) speed, (+) fringe contrast. We first deal with (a). 

In view of the known intensity distribution between the maxima and minima 
given by a Fabry-Perot interferometer, and considering the fact that the width 
of the green line is of the order of 1 A, the range of the interferometer should be 
of the order of 3 to 4A. If the range were smaller, which would be desirable 
since it would give greater resolving power, there would be much danger of 
losing contrast at the “ feet” of the intensity distribution of the line (as was, for 
instance, experienced (16) in a similar problem recently by H. Briick and 
D. A. Jackson). After some tests we decided to use a range about 4A, which 
fixes the spacer in air at 0°35 mm, giving the order of interference about 1320. 
The correctness of this choice was subsequently confirmed, as is shown by Fig. 9. 
Then if a resolving power of ~ 50000 is required, one would need a reflectivity 
of ~g2 per cent. From this figure, using the valuable measurements of Kuhn 
and Wilson (17), one can find immediately that the transmission of one of the 
etalons, using silver coatings of the highest quality, would be o-o45. If instead 
one takes multi-coated etalons, as we did, the transmission may be increased by 
a factor 1°5 or even 2, and for the present purpose this increase is not merely 
appreciable, it may be decisive. 


We have furthermore to deal with an emission line superposed on a continuous 
spectrum, so that a very effective narrow-band light filter is needed to achieve 
the necessary contrast between interference fringes and the background caused 
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by this continuum. ‘The method of diluting the continuum by a dispersion 
arrangement (prism or grating) did not seem entirely satisfactory, for reasons 
mentioned below. We decided therefore that we must use either a filter of the 
type developed by Lyot (18, 19) or a narrow-band interference filter of the 
Fabry-Perot type, as described for instance in (21, 22, 23). A Lyot filter using 
only quartz and calcite would transmit ~ 40 per cent. Narrow-band interference 
filters, if made with silver layers, transmit 20 to 30 per cent, and if made entirely 
with dielectric surfaces can transmit as much as 75 per cent. Interferometer 
and filter together would therefore pass not more than ~ 4 per cent of the incoming 
radiation of the corona. Near sunspot maximum the green line radiates ~ 1 per 
cent of the visible light of the corona (24). Experience by one of us at the 1952 
eclipse in Khartoum showed that in good seeing a white light photograph of the 
corona on a plate of medium sensitivity will extend to about o-5 solar radii from 
the limb, with an exposure of 20 seconds at //120. ‘Taking all these facts together 
a camera working at / 2 or {3 seemed necessary to photograph the corona in 
5303 A as envisaged, with an exposure time of 20 to 60 sec, an estimate which 
proved subsequently to be right. On account of the nearness of the eclipse to an 
exceptionally deep sunspot minimum, the intensity of the green line was expected 
to be unusually low, although no reliable information on its intensity was available. 
But with this knowledge in mind we did our best to achieve high speed and high 
contrast between continuum and fringe intensity. 

It may have been supposed that some gain of light would result from the well- 
known formula 7* (1 — R*), giving the peak intensity of an emission line when 
observed with a Fabry-Perot interferometer. In our case, using multi-coated 
flats with R~ 0°87, 7’~ 10 per cent ( the actual values for the eclipse observations), 
this would lead to a peak intensity of not less than 60 per cent, whilst the continuous 
background from the continuous light which has passed the filter should have from 
T?/(1 — R*) a theoretical intensity of only ~ 4 per cent. But this makes the common 
error of supposing that ail light of the line is concentrated within the instrumental 
line width of the interferometer, which is not the case here (resolving widths 
~o-1 A, actual line widths ~ 1 A), so that a certain decrease of the contrast will 
take place. An example of a similar case can be found in (25). In the present 
instance the interferometer will work more or less like an ordinary spectrograph, 
giving a peak of somewhat higher intensity for the emission line, but in the absence 
of a filter it has the disadvantage that the continuum is diluted only by the very 
small amount of the range of the interferometer, here ~ 4A. Furthermore this 
background is not uniform but depends on the intensity distribution of the 
continuous light of the corona. 

Apart from the necessary aperture ratio of ~1:2 we wished to avoid too 
small a corona image, say of less than 1 cm diameter. In the simplest arrangement 
this would call for a camera with aperture not less than 30 cm and of 60 em focal 
length. But then the dimensions of the narrow-band filter and the interferometer 
would become far too big. We decided therefore to use an optical arrangement 
with intermediate imaging, as sketched in Fig. 1. L,, a normal astrographic lens 
of 30 cm aperture and 340 cm focal length, images the corona in the plane G, 
The focus of the lens L, (aperture = 7-5 cm, /,= 27 cm) lies in G,. At the same 
time L, images the objective L, on the plane J. The interferometer 1s at J 
immediately preceded by an appropriate interference filter F. At J the cross- 
section of the beam is only 2:4 cm, so that the interferometer plates and the filter 
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need be only comparatively small. After J the light is collected by the lens L,, 
a Zeiss Sonnar (aperture = 2-5 cm, f/,=5 cm) of excellent quality. This forms 
simultaneously images of the corona and of the interference fringes on the 
film P. The size of the image at G, is reduced by /,:/, = 1: 5°4, so that the solar 
diameter at P is about6mm. ‘The effective aperture ratio of the whole system 
is ~ 1:2; some small losses will of course occur by reflections, at L,, F, J and L,. 
This arrangement has the advantage that by changing /, and/or /, one can alter 
the image scales of the corona and/or the fringes independently. For adjustment 
purposes and to obtain test fringes conveniently immediately before and after 
the eclipse, the light of a mercury lamp could be reflected by way of the mirrors 
Mi, and Mi, (Mi, being removed for the actual eclipse observation) on to the 
interferometer. ‘The mercury lamp was imaged so that the interferometer was 
illuminated by the comparison light in very nearly the same way as by the corona 
during the eclipse. When using a Lyot or an interference filter in such a set-up 
it is important to remember that by re-imaging with the two supplementary 
lenses L, and L,, reducing the image scale between G, and P by a factor 5-4, all 
angles beyond L, (in the sense of deviation of rays from the optical axis) are 
increased by the same factor. 


G, 


biG. t.--Geneval sketch of optical arrangement. 
Full line : optical path. Broken line ; path for mercury comparison light. 


Narrow-band filter.-For the investigation of a rather faint emission line 
superimposed on a continuous spectrum with a Fabry Perot interferometer 
it is imperative to use a filter with the band width as narrow as possible, because 
contrast, and therefore the visibility of the fringes, in such a case depends entirely 
on the filter. Before modern narrow-band filters were available the obvious 
thing to do was to dilute the background by dispersion with a prism or grating, 
as was actually done in (2) and (3). But unless one is able to use very high 
dispersion, which could be obtained only with difficulty, one can easily see that 
the dispersion method will give worse results than for instance a narrow-band 
light filter of ~ 20A half width, which to-day can be produced in various ways 
without too much difficulty. The dispersion method uses essentially a slitless 
spectrograph, and each half crescent of the corona can be regarded as a very wide 
and not uniformly illuminated slit. A good approximation to the contrast to be 
expected between emission line and continuous spectrum can be obtained easily 
by a graphic integration, taking into account the radial decrease of the light of 
the corona and the dispersion available. A dispersion of 20A per mm can be 
obtained under reasonable conditions with two large glass prisms, but this would 
give lower contrast, even near the limb of the Sun, than light filters of 20 A half 
width. Furthermore, since in the final image the continuous spectrum gives 
a uniform background while the intensity of the emission line decreases rapidly 
with increasing distance from the Sun, the conditions for contrast are becoming 
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more and more unfavourable as we proceed outwards from the limb. For 
geometrical reasons the situation will be still worse near the lower and upper 
limb of the spectrum. 

In preliminary experiments at Cambridge attempts were made to use an 
ordinary Lyot filter already existing at the Observatories. ‘This consisted of 
eight quartz plates and was designed to transmit a band width of ~2A at He. 
It does not seem necessary to have quite such a narrow band for the present work, 
and by shifting the maxima of the different steps in the well-known manner 
(18) it was quite possible to get a band width of ~15A centred on A5303A 
with only three quartz plates, using an additional interference filter with a rather 
wide band and high transmission to remove the unwanted maxima. But since 
only polaroids were available for polarizers, the transmission of this filter was 
found to be far too low. Another difficulty was caused by the rather large field 
required by our double imaging system, which demands a wide-angle type of 
Lyot filter, as described in (18), and one using only calcite instead of polaroids. 
This turned out to be far too expensive for the present experiment. But it ts 
still our belief that probably the best filter for this purpose would be a wide-angle 
L.yot filter, constructed entirely of quartz and calcite. A system with only four 
quartz plates should do quite well. It would probably be possible without 
much difficulty to convert an Ha filter of this type so as to transmit at A 5303 
under reasonable temperature conditions. 

We then tried interference filters of the type described in (1a) and elsewhere. 
Several firms kindly prepared samples of such filters for us which were all tested 
and measured very thoroughly with respect to peak wave-lengths and transmission 
in many different parts of the filter, using the large solar spectrograph (/ = goo em) 
of the Cambridge Observatories. Amongst the filters, those of Baird Associates 
Inc., Cambridge, Mass. and of Bausch and Lomb (26), both of dielectric type, 
showed a remarkably high peak transmission of 60 to 70 per cent, and in parallel 
light when only a very small area of the whole filter was used, half widths of 
30 to 40A. But there were two difficulties: (1) either the wave-length for peak 
transmission was not quite correct and could not in this special case be sufficiently 
adjusted by inclining the filter; or (2) the wave-lengths of peak transmission of 
different parts of the filter were widely different, a defect (sometimes known as 
the wedging effect) from which most interference filters suffer. After some 
trials, however, Hilger and Watts Ltd produced a special interference filter with 
semi-transparent silver coatings, the size being 25 mm square. Using only a 
very small area of a few mm* in parallel light and at normal incidence, the trans- 
mission curve of Fig. 2 was found, with a very small asymmetry on its short-wave 
side. ‘The peak transmission was ~ 25 per cent, the half width as little as 15 A, 
the transmission dropping down in a relatively short distance from the peak to 
well below 1 per cent. Integrating the whole transmitted light, as indicated by 
Fig. 2, and allowing some correction for the low-intensity wings of the curve, 
the total transmission could be estimated to have an equivalent width of ~ 25 A 
of the continuous spectrum. A second and broader transmission maximum was 
found at A4045A. But this excellent performance was again partly spoilt, as 
with the other filters, because the peak wave-lengths for different parts of the 
filter were found to vary as indicated in Fig. 3. This figure shows that use of the 
filter with its full aperture in a telecentric beam would smear out the actual half 
widths over a range of ~60A. ‘This scattering of the peak wave-lengths does 
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not increase the total background intensity allowed through the interferometer, 
but cuts down the transmission for the A 5303 line and hence must have caused 
some loss of contrast because each point of the image at P is formed by rays which 
have passed through a considerable area of the filter. ‘The way in which the 


10 0 10 20 A 


Fic, 2.—Relatiwe transmission of a small area of the Hilger interference filter for A 5303 A 
Absolute peak transmission about 25 per cent, 


Difference of peak tranomission for different parts of the interference filter of Fig. 2, 
in angstrom units, for normal incidence 


~30 -10 OA 
Pia. 4 Shift af peak transmission of the interference filter of Fig. 2 with the angle of incidence, 


peak transmission depended upon the inclination of the transmitted light relative 
to the normal to the filter is shown in Fig. 4. A careful test, taking into considera- 
tion the measurements illustrated in Figs. 3 and 4, led to the conclusion that a very 
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slight inclination of the filter would’ give the best result during the eclipse 
observation. ‘This inclination could be set very reliably before the eclipse by 
using reflected light from one surface of the filter. To cut out the secondary 
transmission maximum near A4045A the filter was combined with an OY 6 
Chance glass filter, whose transmission curve as found by actual measurement in 
Cambridge is given in Fig. 5. 


4000 4500 §000 5500 6500 A 


Absorption curve of the auxiliary glass filter which was 
to be combined with the interference filter of Fig. 2 


Photographic emulsion.\t seemed desirable to provide an effective cut-off 
also on the long-wave side of A5303A._ If the interference filter transmits some 
very small fraction of the light over a long wave-length range, the background on 
the photograph can accumulate so as to diminish the contrast of the fringes 
appreciably. As no suitable glass or gelatine filter seems available for this 


purpose the colour sensitivity of the photographic emulsion used was adjusted 


so as to give a sensitivity as near as possible to A 5303 and then a sharp fall to the 
red. On the other hand, the emulsion had to have just sufficient sensitivity to 
record the green mercury line A 5461 A of the comparison fringes. We are much 
indebted to the Research Laboratories of Kodak Ltd at Harrow for producing 
for this work several samples of such emulsions, sufficient both for tests with the 
large Cambridge spectrograph and for the actual eclipse observations, In a case 
like this the choice of the right emulsion is most important. Eventually an 
emulsion was selected similar to the Kodak Ila-J, a plate (27) of fairly high speed, 
moderately high contrast and of a graininess which seemed to give a good com- 
promise between speed and photographic resolution The resolving power 
of the actual emulsion was tested by very high enlargement and by microphoto- 
metric records of very sharp interference fringes for this purpose. A medium 
yellow dye was added to the emulsion to reduce halation his had been proved 
most effective in earlier experiments, for instance in the photometric work of the 
1952 eclipse, carried out by one of us. Such a dye, of course, may slightly affect 
the spectral sensitivity. Fig. 6 shows the relative colour sensitivity of the emulsion 
from high-dispersion spectrograms taken with the Cambridge equipment. 
For several reasons it was impossible to get the peak sensitivity quite at the right 
wave-length, but the speed proved nevertheless to be suffictently high for our 
work. 
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Interferometer.—Some early tests were made with a single-plate quartz 
interferometer with a thickness of o-5 mm, giving a range of 1-8 A, but this proved 
insufficient for the present purpose. Subsequently a conventional Fabry—Perot 
interferometer was selected; actually the one which had been used for certain 
tests in 1952 (7). The plates were of crystalline quartz, figured to A/40 by Hilger, 
and coated by one of us with four quarter-wave layers of zinc sulphide and three 
layers of cryolite deposited alternately. ‘This coating had been deposited in 
1953 April and the plates had after that been in use for some considerable time 
for high-resolution spectroscopy in the Cambridge Solar Laboratory. As they 
were found to be still in quite good condition in 1954 March we decided to use 
the same coating also for the eclipse. When freshly coated their reflectivity was 
about 94 per cent and their transmission 5 per cent for wave-length 5461 A. 
For separators steel balls of diameter 0-35 mm, produced by Miniature Ball 
Bearing Co., Bienne, Switzerland, were used. They worked excellently and 
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bic. 6.—Spectral sensitinty of the special Kodak emulsion. 


once the interferometer was properly adjusted it could be kept in this adjustment 
for days. For test exposures it was only necessary to remove the interference 
filter F (Fig. 1) and to insert the already adjusted auxiliary mirror Mi,. For 
comparison the light source was a water-cooled mercury isotope lamp produced 
by the Research Laboratories of the General Electric Company, Wembley. 
From interferometric measurements with a resolving power of over 1000000 
it was known (7) that the line width of A5461 from this lamp does not exceed 
0005 A. No correction for the mercury line width or for hyperfine structure 
was therefore necessary. 

The final camera was a standard Zeiss Contax II with the Sonnar / 2 lens 
already mentioned. Although the plate adaptor of this camera with its ground 
glass screen proved very convenient for adjustment purposes, film was used at 
the eclipse, partly because it allowed a much more rapid succession of exposures 
than could be made when changing plate-holders, and partly because the focus 
for the Contax miniature plate-holders was found to be not as stable as with the 
film. 

In a long series of test exposures made before the expedition left Cambridge 
the fringe definition was found to be excellent, the instrumental line width being 
j of the range (~o'15 A). The resolving power was only slightly below the 
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theoretical value. The range at A5303 was determined with the solar spectro- 
graph using so-called channelled spectra, and found to be 4:04 A (28, 29, 30, 31). 
Such spectra can be produced easily by a continuous light source with the inter- 
ferometer in front of the large Cambridge solar spectrograph. Furthermore, 
before leaving Cambridge the whole optical arrangement was tested thoroughly 
in the laboratory with all parts exactly in the position in which they were used 
afterwards during the eclipse, using a kind of “ artificial corona” as well as the 
full Moon. For the laboratory test the large coelostat system (f= 1800 cm) 
of the Cambridge solar spectrograph was used as a collimator. An “ artificial 
corona” of correct size was set in the focal plane of the 1800 cm lens. The 
outgoing telecentric beam was then reflected from the auxiliary mirror back on 
to the optical arrangement of Fig. 1. ‘This artificial corona could give continuous 
light and superimposed, when required, the green mercury emission line also. 
These tests were found most valuable as a check on vignetting, false light in the 
apparatus, and on exposure time. 

Contrast.—-A general idea as to exposure time was obtained by tests at 
Cambridge on the full Moon. However, in an experiment of this kind it is also 
most important to estimate the contrast to be expected between the fringe and the 
continuous background. Here lay the largest unknown factor in the experiment, 
for not all the necessary data were known sufficiently accurately. A rough 
estimate could, nevertheless, be made in the following way. ‘The intensity ratio 
of the A 5303 line to the continuous spectrum of the corona was measured carefully 
by Grotrian (24) on slit spectrograms taken at the 1929 eclipse (maximum of 
solar activity in 1928-4). He obtained a ratio about 27:1. We assume the ratio 
to be of the order 20: 1, which agrees sufficiently well both with Grotrian’s work, 
with the measurements of Lyot, and with estimates from other eclipses (32, 33). 
Our interference filter had an effective band width of ~ 25 A, and taking the range 
of the interferometer to be 4A there is a decrease in contrast by a factor ?~6. 
As we have seen that the wave-length of peak transmission of the interference 
filter is not quite constant over its surface there is an additional decrease of contrast 
which may be estimated at perhaps a factor of 2. For an ideally sharp line the 
contrast between minimum line intensity and continuous background should be 
1~R*/(1—R)*; with R=o-87, the value we may assume for the eclipse 
observations on Koster, this gives an increase of the contrast of ~15. But 
since the corona line is actually much broader (~1 A) than the instrumental 
width (~ 0-2 A) and has a characteristic Doppler broadening shape, the contrast 
is again decreased by a factor of say6. ‘Taking all together this brings the contrast 
of line to background to about 20 x 15/6 x 2x6=4:1. But this value relates to an 
average phase of the solar cycle, whereas we were working at a very low minimum 
in an unusually unfavourable cycle. No measures were available to show how 
weak A 5303 really had become, but it seemed reasonable to expect about 4 
of the value above (34, Fig. 12). So we had to be prepared for a contrast as low 
as 1:5, and this subsequently proved to be of the right order. Using an artificially 
broadened emission line, it was confirmed by special tests in Cambridge that the 
photographic emulsion and the method of developing employed later at Koster 
would give reasonable results for rather low contrasts of this order. 

Observations on Syd Koster.—Vhe Cambridge eclipse party was located on 
the small granite island of Syd Koster about 15 km off the west coast of the 
Swedish mainland not far from Strémstad, We considered it unwise to use the 
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sensitive interferometer optics in a temporary hut or tent as is usually done with 
apparatus at a solar eclipse. Fortunately a private house with a substantially 
built cellar on granite ground was available where we could set up the apparatus, 
shielded beautifully from weather and temperature changes, and where in fact 
we had most of the essential features of a good laboratory. Water and electricity 
supplies were available in the house and the cellar was easily blacked out to give 
a photographic dark room and a darkened room for optical work. ‘The Sun’s 
light was reflected through one of the cellar windows by a 40 cm coelostat and an 
auxiliary mirror. ‘The slight loss of light by this double reflection was a small 
price to pay for the good laboratory conditions under which we were able to work. 
The whole arrangement proved very satisfactory and can be strongly recommended 
for similar purposes, provided seeing requirements are not too critical. The 
temperature changes over the whole day were always lower than 1 deg. C and it was 
a great advantage to be able to make adjustments and tests quite independently 
of weather conditions, ‘The whole equipment was set up quite simply as in the 
Cambridge laboratory on standard optical bench rails and fittings (reproduction 
in (3§)). Instead of pillars of brick or concrete some of the expedition’s bigger 
cases were used, loaded with rocks, and these proved to be very solid. 

We had some doubts as to whether to use the interferometer adjusted with 
the fringe system centred on the Sun or in an eccentric position. ‘The centre 
position would have the advantage that one might obtain perhaps one or two 
fringes right round the corona at more or less constant photographic density, 
and could obtain measurements of possible radial velocities much more easily. 
On the other hand, it is well known that it is not practicable to adjust an inter- 
ferometer so that a given fringe has exactly some desired diameter. ‘Thus it 
could easily happen that the first visible fringe would just miss the most interesting 
partofthe corona. In any case not more than about two fringes could be expected 
to show. ‘The eccentric position on the other hand gave far more chance of 
covering the whole corona with a system of fringes, without much loss in resolving 
power. After experience one of us had gained in 1953 with similar observations 
on the coronagraph of the Pic du Midi, and since our exposure times were 
relatively short, we decided to try both arrangements and to change from one 
adjustment to the other during totality itself. 

The fringe system was centred on the solar image a few minutes before second 
contact, and then, after totality had started and the first three frames had been 
exposed, the fringes were moved to the eccentric position by inserting a prepared 
wedge under the interferometer mounting. ‘Two more frames were taken with 
this adjustment. 

During totality, the coelostat and auxiliary mirror outside the house were in 
charge of Mrs von Kliber, whilst the observers inside the cellar took the eclipse 
exposures and the calibration and test exposures. ‘The corona was photographed 
through a thin high layer of irregular cirrus cloud under which was another thin 
layer of fairly fast-moving broken-up cumulus. The clouds interfered rather 
seriously with the first and second exposures; for the rest the seeing must have 
improved considerably. 

Shortly before totality a set of intensity marks was photographed on the 
beginning of the film using an accurately calibrated Hilger rhodium step filter 
(using L, for imaging) with an artificial continuous light source, and through an 
interference filter similar to the one used in the eclipse itself. ‘The wave-length 
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range of spectrum and the exposure times were quite close to those of the actual 
eclipse photographs. Then a few minutes before the eclipse, several sets of 
comparison fringes from the mercury lamp were exposed on the same film in 
the centred fringe position. During totality the following set of exposures was 
made on the corona: 


Frame 1 exposure 10° 
2 20° } centred fringes 
3 60* 

Frame 4 exposure 25° 


eccentrically placed fringes 

Immediately after the eclipse another set of mercury test fringes was 
photographed on the same film, this time in the eccentric position, before the 
interferometer had been touched. ‘The interferometer was then brought back 
to the centred adjustment and additional mercury fringes were taken exactly 
as at the beginning. Lastly a duplicate set of calibration marks was printed 
on the end of the film. As the contrast of the A 5303 line against its background 
was expected to be very low, the film was developed for high contrast in a tank 
with continuous agitation, for five minutes at 20 deg. C in undiluted Kodak Digb 
developer. ‘The photographic quality of the film (and of a second control film) 
was found, after the development, to be very good. The eleven different 
photographs of the mercury fringes showed that the optical adjustment had 
remained good throughout, while the intensity marks testified to the general 
photographic quality. 

Films. Frame 1.— This frame was waken at the beginning of the eclipse with 
10° exposure time and seems to have been considerably affected by cloud, It is 
very under-exposed and far weaker than frame 5 taken later with the same exposure 
tume. The innermost part of the corona is clearly visible but its image nowhere 
reaches the distance of the third interference ring, the first which would have 
been visible, i.e. no fringe could be expected to show and none is visible. 

Frame 2.—The density of this frame, exposed for 20%, is still rather weak 
compared with the fourth frame, which had a similar exposure time, and clouds 
must still have interfered. ‘The prominences P, and P, (Figs. 7, 8) can be clearly 
seen and offer a good method for correctly orientating the photograph. On the 
east side the corona extends just to the third fringe which can be traced very 
weakly but quite clearly, as indicated in Fig. 7. ‘The fringe seems to be sharpe: 
in the N.E. quadrant and more diffuse in the 5.E. quadrant. On this frame 
the west side of the corona is much weaker and only a trace of the third fringe can 
be seen in the N.W. quadrant. ‘The corona nowhere reaches the fourth inter- 
ference ring. It does not seem worth while to measure this film. 

Frame 3.—This frame had the longest exposure, 60%, and as transparency 
of the sky had now become much better it shows the highest density. The 
centre of the fringe system coincides well with the Sun’s centre. The first 
and second rings were calculated to fall well within the solar image, the second 
being very near to the Sun’s limb. In spite of the rather high density, the third 
fringe can be seen distinctly on the east side of the corona, as indicated in Fig. 7, 
which was drawn from this film. ‘This fringe appears to be rather sharp in the 
N.E. quadrant and a small kind of knot can be seen at P.A. 60°, of the kind well 
known to observers with the Lyot coronagraph. In the region of the equator 
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the contrast to the background density becomes very low, so that the fringe can 
be seen only under low magnification. In the S.E. quadrant the contrast increases 
again to fairly good visibility; the fringe seems here to be definitely broader and 
more diffuse than in the N.E. quadrant. Also the widths seem to vary a little 
from place to place and these variations are apparently not caused by slight 
intensity variations of the underlying corona. As all the comparison fringes from 
the mercury isotope lamp taken shortly before and just after the eclipse on the 
same film and with the same adjustment proved to be excellent, this behaviour 
of the fringes is probably real. 


Fic. 7.--Shketch indicating the centred interference fringes visible on Frame 3. Dotted line : 
outer contour of the corona image on this frame. Thin line : position of the third and fourth 
interference fringes. Heavy line ; visible parts of interference fringes from 45303. F : centre 
of interference system, P,, Py, Ps: positions of prominences. C,, Cy: positions of second and 
third contact, 


In the surroundings of the north and south pole no fringe is visible, but it 
cannot be expected, because at this distance from the Sun’s limb the film is 
practically clear, as can be seen from the isophotes indicated (dotted line) in Fig. 7. 
‘The density of the corona image on the west side is higher than on the east side. 
In the region occupied by the third interference ring the corona is definitely 
over-exposed. Nevertheless the third ring can be traced quite well in the N.W. 
and S.W. quadrants if a low magnification is used, although with low contrast. 
‘These parts of the fringe have about the same widths as in the S.E. quadrant but 
since the density is too high and the contrast too low they cannot be measured 
here with sufficient accuracy. 

The fourth ring at the east side is clearly visible at the outer edge of the corona 
image where the density of the film is already too low for reliable measurement. 
On the west side this ring is hardly to be seen at all in spite of the fact that the 
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density there is still quite sufficient. ‘This is in agreement with our other frames ; 
all of them show the fringes far more distinctly on the east side than on the west 
side of the Sun. 

Frame 4.—Fig. 8 was sketched from this frame. It is well exposed; indeed the 
inner parts of the corona are rather over-exposed. The eccentric fringe system 
is easily and clearly visible on the east side of the corona. ‘The density on the 
west side is higher and the (photographic) contrast of the fringes therefore much 
lower. From micrometric measurements of the visible parts the whole system 
of eccentric fringes can be easily constructed as indicated on Fig. 8, their centre 


b1G. 8.--Sketch indicating the eccentric interference fringes visible on Frames 4 and 5. 
figures at the right-hand side indicate the fringe numbers, Other notation as for Fig. 


being at F. ‘The well-known law that the squares of the radii are in arithmetical 
progression is satisfied by these measures extremely well. In Fig. 8 all parts of 
the interference fringes which could be seen and measured without difficulty 
on Frame 4 are drawn in heavy line. ‘The dotted curve again indicates roughly 
the outer isophote on this film. It is clear from this frame, and also from the next 
Frame 5, that the fringes avoid the region near the equator where as far as photo- 
graphic density is concerned they should show up easily. ‘The same is true for 
the polar regions. This behaviour is just what one would expect if one takes 
into account what we know of the average intensity distribution of the green 
corona line in this phase of the solar cycle (14). Near the Sun’s limb on both 
east and west sides the density of the film is very high, the (photographic) contrast 
therefore low, and it gets progressively more difficult to trace the fringes as the 
limb is approached. This frame has produced most of the measurements of 
Table I, 

Frame 5.-——This is also a very good exposure, its density being much higher 
than that of Frame 1 which was exposed for the same length of time (10*) but which 
apparently was taken through clouds. All the fringes indicated in Fig. 8 can 
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clearly be seen; but here they can be traced practically down to the limb of the 
Sun. Again the density of the west side is higher and the contrast of the fringes 
is there diminished. Measurements of this film are also included in Table | 
and form a valuable supplement to those determined from Frame 4. 
Measurement of the films.—-The various sets of mercury comparison fringes 
taken on the same film before and after the eclipse were measured visually with 
the micrometer and then with the recording microphotometer. Both pre- and 
post-eclipse images proved to be of the best quality, in keeping with those obtained 
in the tests in Cambridge and on Koster before the eclipse. ‘The interferometer 
therefore was in perfectly good condition during the eclipse itself. ‘The half- 
widths of the fringes, which were examined specially in the neighbourhood of 
the fifteenth fringe, where most of the later measurements of the corona line were 
made, was found to be o-21 A, i.e. about g of the range. This is equivalent 
to a resolving power of 26000 and, from resolving power = 3.m ,/R/(1— R), in 
good agreement with an effective reflectivity R of 87 per cent. It is also in good 
agreement with the theoretical half-width (36) which from AA(1—R)/7.R, 
where AA is the range of 4°04 A, would be expected to be 0:18 A. These extremely 
sharp mercury isotope interference fringes give the instrumental profile 
immediately without further correction. Fig. g shows two of them drawn to 
correct scale and in correct relation to their range. 


T 
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Fic. 9.--Example of two line contours of 4 5303 from Frame 4 drawn in correct relation to 
the range (4°04 A) of the interferometer. The dots indicate actually observed values. The full 
curve indicates the line contour after correcting for the limited resolving power of the photometer. 
The line contour of the mercury isotope comparison fringe is given by the thin line. 


From Frame 3 the widths of the third fringe of A 5303 could be measured 
with the micrometer in different parts of the corona, ‘The results are in Table I. 
‘The density of this frame is too high to give reliable microphotometric records, 
so that only the micrometric measurements are given here. 

On Frames 4 and § micrometric measurements of the fringe widths could be 
obtained without difficulty for many selected points. ‘hese measurements are 
to be found also in Table 1. ‘They give immediately, without further reduction, 
fairly reliable values for the widths of the green line at different parts of the corona. 
Lyot (9, 10) made his first measurements of this line with a slit spectrograph in 
a similar manner. As the range of the interferometer was known accurately 
ana the separation of the fringes can be measured also very well, the scale value for 
converting the measurements into angstroms provided no difficulty. On the 
other hand, the high density gradient within the image of the corona was a source 
of difficulty in the measurements. ‘This was of course expected, but since we had 
to face so many other uncertain factors concerning exposure time and contrast 
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we had deliberately avoided taking special precautions against this effect, so as 
not to introduce more complication than necessary into the experiment. It was 
for this reason that we had planned to get frames of different exposure times. 

But even if one takes into account the density gradient in comparing Frames 4 
and 5 (which have altogether different densities) with each other there is not much 
doubt that the line appears to be much broader quite near to the limb of the Sun 
than further out. On the average this is confirmed by tracings made on several 
fringes with the Cambridge self-recording microphotometer. Owing to the 
fact that the coronal image was rather small, the photographic density high and 
the fringes narrow, it was found difficult to obtain these tracings without loss of 
resolving power. ‘The graininess of the emulsion on the other hand did not cause 
much trouble. ‘The many intensity calibrations on the film permitted reduction 
of the photometer tracings to intensities in the usual manner. As mentioned 
already, the fringes are superimposed on a background with a rather steep density 
gradient. A difficulty was, therefore, encountered similar to that in the photo- 
metric investigation of faint components very near to a strong spectrum line. 
But in spite of this, by carefully taking into account the slope of the background, 
fairly good intensity curves could be obtained of suitably selected parts of the 
fringe. ‘The possibility was even considered as to whether to apply a correction 
for the density gradient with relation to width of the photometer slit. ‘The 
necessary procedure was described long ago by Rayleigh (37). ‘This correction 


for x» can be written in the form 
(1 d*r 
3 dx* Le 


where « is the half-width of the slit of the spectrograph, and r(x) the coefficient 
of transmission of the emulsion at the point x. But the correction was found to 
be very small and not very certainly determined, so it was omitted. Besides 
recording the fringes of the corona line, tracings were made of the corresponding 
eccentric fringes of the mercury comparison, with the same setting of the 
microphotometer, which gave the instrumental width as produced by the 
photometer. ‘The intensity distribution in the fringes of both corona line and 
mercury lamp could be represented sufficiently well by an error curve, 


aexp{ 


Using the instrumental half-width, the observed half-widths of the corona line 
were corrected in the usual way (38), by taking the difference of the respective 
squares of the denominator in the exponential expression. Because of the 
relatively high resolving power provided by the interferometer these corrections 
were rather small, and in good agreement with corrections derived for the special 
case of a Fabry-Perot interferometer by Minkowski and Bruck (39). The 
corrected half-widths are in Table I, ‘This procedure for the reduction is open 
to some criticism. It would have been better to approximate to the observed 
curves by a two-parameter type of expression, to take better account of the wings. 
‘This more elaborate reduction definitely becomes necessary if one has to deal with 
the intensity distribution of a Fabry-Perot pattern for reflectivity <80 per cent 
where the wings are no longer negligible. But in our case the wings cannot be 
determined with sufficient accuracy anyhow and an application of a Voigt function 
as suggested by v. d. Hulst and Reesinck (40) would not have been justified. 
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The photometrically derived widths are on the whole in fairly good agreement 
with the purely micrometric measurements. But there seems to be a slight 
tendency for the photometric values to be larger than the micrometric ones. 
Fig. g shows an example of the fringe contours in two successive fringes from 
Frame 3, drawn in correct proportion to the range of the interference system 
and normalized to peak intensity 100. The profiles of the mercury isotope 
comparison fringes as produced by the interferometer are shown by the thin 
curve. The observed points for the corona fringes are given by dots and the 
heavy lines denote the contours after correction had been applied for the 
instrumental profile. ‘Taking into account the rather small scale of our image, 
its high density and the limited contrast of the fringes, photographed under most 
unfavourable conditions, we do not feel justified in putting too much weight on 
these preliminary measures or in drawing any far-reaching conclusions from 
them. Also it seems inadvisable to attempt to deduce macroscopic movements 
by Doppler shifts from this material, for this really needs a larger image, and 
probably would be better done with the centred fringe system. The present 
measures however seem to indicate a decrease of the line widths with increasing 
distance from the Sun, an effect also indicated in coronagraph measurements 
by Waldmeier (13, 14, 15). On the other hand, other observers’ measures of 
the red line 46374 seem not to show such a decrease (12), though this would also 
be expected. ‘The line widths appear rather large on some parts of the fringes 
and this behaviour is revealed simply by visual inspection of the fringes under 
alow power. In spite of the not very favourable photometric conditions 
mentioned above, it seems improbable that there would be any large systematic 
falsification of the widths by photographic effects; and there seems no doubt 
whatever that the interferometer worked perfectly well. ‘The mean error of the 
values given in ‘Table I must be estimated as at least o-1 A, and even larger for those 
values where the contrast was very low. But in any case even where the values 
seem to be somewhat doubtful there always seems to be the tendency towards the 
larger value. 

If one supposes that the width of the line is due entirely to temperature 
broadening, the temperature can be derived easily from the expression 7'= uf? 2k 
with the molecular weight 56 (Fe), the gas constant and the velocity derived 
from Doppler effects from the half-widths. ‘This would lead to a temperature 
T of ~ 2:2 x 10° to 5 x 10° with most of the values round about 2:5 x 10° in the 
range of 1°3 to 1-05 solar radii from the Sun’s centre. In a recent paper Allen (41) 
has pointed out that temperatures derived by direct spectroscopic measurements 
of line widths seem to be slightly larger than those derived by other methods 
This may point to high turbulent velocities in the corona on a macroscopic 
(but still rather small) scale. But with most of the previous spectroscopic values 
a rather high correction had to be applied because of the large instrumental widths, 
while a simple Gaussian approximation was used for the observed curve, so that 
some of them may have been under-corrected (40). 

Spectroscopic measurements of the width of the green line with a coronagraph 
were obtained as early as 1931 by Lyot (g) in the second order of a Rowland plane 
grating, giving a dispersion of 7-5 A/mm, but a resolving power of only 8000. 
Further measurements with a coronagraph have been published by Waldmeier 
(14, 15). He found the half-widths to be about 1-3 A, but used a resolving power 
as low as 4000 (~ 1°3 A), so that his measures required rather large corrections. 
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In an active region of the corona he found in 1945 (1§) half widths between 1-56 
and 0o-96A,but again the resolving power was low (~1 A), and the necessary 
corrections therefore were still rather high. In 1942 Lyot, using a grating 
spectrograph attached to the Pic du Midi coronagraph, obtained resolving power 
20000 (0°27 A) at a dispersion 4A/mm (12). Measurements of the green line 
with this instrument have not yet been published, but for the red line (6374 
Dollfus has given a number of measurements with a resolving power of 0-3 A, 
giving a temperature of 2 x 10° at 1-1 radii from the Sun’s centre. Measurements 
of Righini (42) at the eclipse of 1936, with a slit spectrograph giving a width of 
3°8A for the green line, are not reliable, as the author himself points out, because 
of the low resolving power of the instrument. At the 1941 eclipse Kaliniak (4) 
obtained some measurements from a rather faint interference pattern from the 
green line, using a Fabry~Perot interferometer in front of the slit of a prism 
spectrograph, ‘The range of the interferometer was ~ 8 A, reflectivity ~ 70 per 
cent, and the resolving power, therefore, could not have been larger than ~ 5000 
or approximately 1A. In this case again a rather large correction for the observed 
profile was needed which actually could have been obtained only by extrapolation 
as no suitable exposed reference line was available. He has derived line widths 
of 1°4A from his observations, Apparently no other measurements of the widths 
of the green line during an eclipse have so far been published. 

General comments.—The main aim of the present work was to establish the 
practicability of this method of investigating the corona during an eclipse. With 
the experience outlined here there is no more doubt that valuable observations 
of this kind can be made at most eclipses, assuming exposure times of 1 minute 
or more are possible. Desirable improvements, some of which we had omitted 
deliberately to avoid too many complications for the first test, should be outlined 
here: 

1. A factor of 2 in exposure time could be gained if an emulsion could be 
obtained with its peak sensitivity centred exactly at A 5303 (see Fig. 6), but no 
emulsion with bigger graininess should be chosen. ‘This might be achieved, 
perhaps, by sensitizing dyes of the cyanin or carbocyanin type, but the emulsion 
has to be carefully balanced for spectral and general sensitivity and account 
has to be taken of the dye added to reduce halation. 

2. Another factor of between 2 and 3 could be gained if a special dielectric 
filter were used. Alternatively a wide-angle Lyot filter would be very useful and 
might give an improvement by a factor of 1-6. 

3. The corona should be imaged on at least double the scale used here, the 
easiest method being to use a longer focal length for L,. 

4. It is most important that the steep gradient near the innermost part of the 
corona should be smoothed out to facilitate the attainment of higher accuracy 
in the photometry. ‘This could be done by introducing near G either a specially 
designed circular absorption filter as pointed out long ago by Michailov (43) 
or alternatively a rotating sector of suitable form. 

5. Our observations were made at a most unsuitable phase of the solar cycle 
when the green line was extremely weak. At the time of the eclipse we knew 
from routine measurements, kindly transmitted by radio from the Wendelstein 
(44), that even under exceilent sky conditions great difficulty was encountered 
in detecting the green line with the Lyot coronagraph. ‘The average intensity 
variation of the green line with the solar cycle is well known from many thousands 
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of observations made with Lyot coronagraphs (34). In favourable circumstances 
the line may be 10 or more times brighter than at the time of our observations. 
The continuous light of the corona (K-corona) seems also to alter with the phase 
of the solar cycle, the brightness at maximum may be about twice that at minimum 
(45, 46). This would mean that near maximum solar activity not only the 
intensity but also the contrast of the green line with respect to its continuous 
background should rise appreciably, making interferometric investigations 
correspondingly easier. It is difficult to estimate how much we may have lost 
on Koster by thin cloud, even during the exposures of Frames 3, 4 and 5, but the 
experience of other observers nearby supports our belief that we may have lost 
a factor of at least 2 in intensity. Altogether, therefore, there should be good 
chances for observations of this kind at a phase of higher solar activity. 

6. The intensity of the red corona line 46374 according to (34) should be 
sufficient to give fringes during an eclipse, provided a suitable interference filter 
were available. We had hoped for such observations on Koster, but had 
unfortunately been unable to obtain the right filter in time. The observations 
might be made with the same arrangement as described above by simply changing 
from a green to a red interference filter during totality itself. 

7. If an interference filter is used it should be coated, if possible, directly on 
a suitable coloured glass filter to cut out undesirable secondary maximum. ‘This 
would avoid having an extra filter in the optical arrangement, which may always 
give trouble by unwanted reflections. As there is a tendency for these inter- 
ference filters to change their wave-length of peak transmission with time, they 
must be checked by a spectrograph of rather large dispersion shortly before 
totality. It should be possible now to obtain such filters with an equivalent 
width of about 10 A, a peak transmission ~ 50 per cent and a peak to off-peak 
transmission ratio of 10* (42). 
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Summary 


Radial velocity determinations have been made on a uniform system for 
55 cepheid variables in the Southern Hemisphere. Velocity curves have been 
drawn from 830 observations secured at phases well distributed over the light 
curves of the variables, and the normal velocities have been determined 
graphically by the method of equal areas 

The plates for radial velocity were obtained with the Cassegrain spectro- 
graph of the 74-inch Radcliffe reflector, the dispersion of the spectrograms 
being 49 A/mm at Hy. The velocity system to which the determinations 
refer has been examined by means of 39 spectrograms of 9 standard velocity 
stars, mostly giants, and 46 spectrograms of 9 supergiants obtained during 
the cepheid programme. ‘The Radcliffe velocities of the standard velocity 
stars were found to be in very satisfactory agreement with velocities on the 
Lick System 

The radial velocity observations presented in the paper complete an 
extension of Joy’s work to southern cepheids, with photographic magnitudes 
brighter than 10-5 at minimum light, contained within a hitherto largely 
unobserved quadrant of galactic longitude 


Introduction.— The results of a comprehensive programme of spectroscopic 
observations of cepheid variable stars, carried out at the Mount Wilson Obser- 
vatory by A. H. Joy (1), were published in 1937. Joy's programme completed 
a survey of cepheids north of declination — 40° and brighter than the fifteenth 
photographic magnitude at minimum light. Radial velocity curves were 
obtained for 106 stars brighter than 14-5 at minimum light from observations 
conveniently distributed in phase, and the normal velocities were found by the 
method of equal areas from schematic velocity curves drawn through the 
observations. Joy also gave provisional normal velocities for twenty-two additional 
cepheids for which velocity curves could not be drawn, either through lack of 
sufficient data, or uncertainty in the ephemerides of the light curves. ‘This 
brought the number of cepheids with known radial velocities to one hundred and 
fifty-five. In a discussion which followed the publication of the observations, 
Joy (2) gave a table of velocities for 156 cepheids, one additional star having been 
brought into the list. Apart from some redeterminations by Sanford (3), Struve 
(4), and Herbig and Moore (§), it appears that no further contributions have 
since been made to the list of cepheid velocities. 

Joy’s list contains very few stars in the sector of galactic longitude extending 
from about /= 230° to the longitude of the galactic centre near /=425°. ‘This 
sector contains many cepheids, but only a few of the brightest had been observed, 
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and in some cases the observations from which the normal velocities were calculated 
had been obtained at phases inadequately distributed over the periods. It has 
not been possible until recently to secure further observations in this part of the 
galaxy on account of its southern declination. In this communication an account 
is given of a programme of observations of the radial velocities of fifty-five cepheid 
variables, with photographic magnitudes brighter than 10-5 at minimum light, 
contained in this largely unobserved quadrant of galactic longitude. 


1. Spectroscopic equipment and observational procedure 


The observations were made at the Cassegrain focus of the 74-inch Radcliffe 
reflector, using the two-prism spectrograph (6) designed by Professor Redman. 
‘The prisms of the spectrograph are of light flint, and are set for minimum 
deviation at A 4200 A which is the central ray through the camera system. ‘There 
are four camera combinations giving linear dispersions of 21, 29, 49 and 86 A/mm 
at Hy, designated for convenience the a, b, c and d cameras respectively. The 
spectrograph is provided with a curved slit for use with all cameras. A sensitive 
thermostat ensures that the temperature of the spectrograph is adequately 
controlled. ‘The specifications for the prisms and the a, b and c cameras were 
chosen so that the spectrograms would closely resemble those obtained with the 
Cassegrain spectrograph of the Dominion Astrophysical Observatory, Victoria, 
B.C, 

| 


Approximate exposure times for Radcliffe spectrograms with dispersion 49 A/mm at Hy 


Plates Photographic Exposures 
magnitudes in minutes 


Ilford Process 


Ilford Special Rapid 
Ilford Special Blue Sensitive 


Ilford Zenith Astronomical 


Kodak 


The c-camera was used throughout the observing programme, the standard 
slit-width being 0°075 mm, projecting to 0-018 mm on the plate. ‘The stars were 
trailed to and fro along the slit in right ascension to give a spectrum of width 
o20mm. Approximate exposure times for the photographic plates used in the 
programme are given in Table I for these standard widths. ‘The limiting photo- 
graphic magnitude for cepheids in the observing list was set somewhat arbitrarily 
at1o™-5. Increased accuracy could have been gained by using a higher dispersion 
for the brighter stars, but this could not have been achieved without loss of 
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observing time incurred by changing from one camera to another. Accordingly, 
it was decided to use the same camera for all stars, regulating the exposure times 
for the brighter stars by the use of slower photographic plates. In this way the 
velocities of the brightest cepheids were automatically determined with the 
greater accuracy made possible by the fine grain and high contrast of the slower 
emulsions. At the same time, uniformity was maintained in the measurements 
and reductions by using spectra with the same dispersion throughout the 
programme. 

An iron arc comparison spectrum was used, the iron arc lines being impressed 
on each side of the stellar spectrum both before and after each stellar exposure, 
and in some instances during the exposure as well. ‘The separate exposures on 
the iron are were of nearly equal duration. ‘The practice was regularly followed 
of checking the collimation of the iron arc beam with due care before observations 
were made with the spectrograph. In making the observations an attempt was 
made to minimize the effect of atmospheric dispersion across the slit of the spectro- 
graph by guiding on the blue core of the stellar image on the slit by means of a 
filter clamped over the guiding eye-piece. ‘This was only possible for the brighter 
stars in the observing list, and was only necessary when the excellence of the 
seeing conditions demanded its use. 

Most of the plates were calibrated for intensities by means of a tube sensito- 
meter employing a filter for average calibration in blue light. 


2. Measurement and reduction of the spectragrams 

The spectrograms were measured for radial velocity on the projection 
micrometer of the Radcliffe Observatory. ‘This machine, the design of which is 
due to Redman (7), gives a total magnification from plate to the projected image 
of about forty diameters, the projected field being about 2-5 mm of the plate. 
The plates were measured in the usual way, both in the direct and reverse 
positions. 

Iron arc lines.—A total of thirty-five arc lines were included in a selected list 
of comparison lines in the wave-length range 3920A to 4529A. Both strong 
unblended lines and weak lines were included in the list to ensure that the Hartmann 
correction curve would be sufficiently well defined for all likely densities of the 
are spectra. ‘The arc exposure times were chosen to give a measurable density 
for the line A452gA. Weak, rather than strong, arc densities were favoured as 
they were found to be more suitable for measurement by the projection method. 

Stellar lines.—The wave-length standards recommended by the Dominion 
Astrophysical Observatory (8) for the determination of radial velocity for solar- 
type stars from low dispersion spectrograms have been used. ‘These wave-length 
standards give correct dynamical velocities for objects in the solar system, and 
stellar velocities on the Lick System for spectral types 4 to K8 

Reduction.— The reduction of the plate measurements was carried out in the 
usual way. The constants in the Hartmann formula were found from the 
measured positions of iron arc lines on a sample of spectrograms taken under a 
variety of conditions. ‘The Hartmann formula for the average conditions was 
found to be entirely satisfactory for the reduction of all the spectrograms. ‘The 
following procedure was adopted. ‘The wave-lengths of the iron arc lines, and 
the stellar lines for zero velocity, were first transformed into micrometer readings 
by means of the Hartmann formula. The difference between the observed 
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and computed micrometer readings for the arc lines then gave a correction curve 
for the reduction of the dispersion on the plate to the standard dispersion given 
by the formula. ‘The Doppler displacements of the stellar lines from the tabular 
values of the micrometer readings for zero velocity were then converted into 
velocities by means of factors computed for each line by differentiation of the 
Hartmann formula. Finally, the velocities were reduced to the Sun by correcting 
for diurnal rotation whenever necessary, and for annual motion by means of the 
24-hourly changes in the rectangular coordinates of the Sun given in the Nautical 
Almanac. No correction for curvature was applied to the velocities as an investiga- 
tion by Wesselink (6) showed that the curved slit of the spectrograph compensated 
adequately for the curvature introduced by the optical system, the residual 
correction for curvature being less than o-1 km/sec in the photographic region. 


3. Standard velocity stars and standard wave-lengths 


The extent of the list of standard velocity stars for inclusion in the observing 
programme was determined by the following considerations. In the first place, 
the observations were begun soon after the spectrograph had been installed at the 
Cassegrain focus of the 74-inch reflector, and a periodic check on the performance 
of the spectrograph was necessary, particularly after adjustments had been made 
to the spectrograph itself or to the secondary mirror of the telescope. Secondly, 
the suitability of the Victoria wave-length standards for radial velocity work 
with the Radcliffe spectrograph required verification on account of the dependence 
of the effective wave-lengths of blended lines upon the dispersion of the spectro- 
graph and its diffraction pattern. ‘Thirdly, the wave-lengths had not previously 
been shown to be entirely satisfactory for supergiant stars throughout the range 


Il 
Standard stars 


Radial velocities and differences 


General | Radial velocities Differences 
Spectral 
Catalogue | 


number | Radcliffe Lick No. | Gen. Cat 


tyre 


dF3* | 10083 

|| 933! 
gGr | 33031 
32401 
suKko | 6628 
wiki 73981 
| 53441 
| Ss53t | 


* This star was originally selected on the basis of the HD classification Fs 


+ L.A.U, standard velocity stars in equatorial regions 


only one star,« Persei with spectral type cl’ 4, having been examined when 
the Victoria wave-length standards were originally determined (8). A list of non- 
variable c-stars covering the entire range F4-K8 was therefore necessary in order 
to ascertain if the recommended wave-lengths gave satisfactory velocities for 
these cepheid-like spectra. ‘Two sets of standard velocity stars, primary standards 
and supergiants, with spectral types distributed as well as possible throughout 
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the range F4~-K8 and located in both equatorial and southern regions were chosen 
in order to meet these requirements. 

3.1. Primary standards.—The radial velocities of nine standard velocity 
stars covering the range F5-K4 were determined from 39 spectrograms. No 
evidence was found for flexure effects in the individual plate velocities, in 
agreement with other tests for flexure which have since been carried out (6). 
The mean velocities are given in Table II, with the velocity differences Radcliffe 
minus Lick (R—L) and Radcliffe minus General Catalogue (R—C). The nine 
stars are listed in the Catalogue as having a-quality velocities. 

The mean difference R—L for the standard stars considered is o-okm sec, 
the mean differences R—C being —o-zkm/sec. Since the Mount Wilson 
General Catalogue gives velocities on the Lick System, two sets of differences 
would be the same for a sufficiently large sample of stars. ‘The stars in the 
sample are almost entirely giants, and the inference is that Radcliffe velocities 
for giants do not differ significantly from velocities on the Lick System. ‘The 
velocity differences for the stars examined are shown in Fig. 1. 
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Fic. 1.—Velocity differences Radcliffe minus Lick (R—L) and Radcliffe minus Mt Wilson 
General Catalogue (R—C) for standard velocity stars. 
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Radial velocities of supergiant stars 
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General Radial velocities 
Catalogue 
number Radcliffe 


Spectral 
type 


Lick No | Mt Wilson No | 
| 


cF3 11022 5 3 
cF4t 7287 4 | 
cF6 7734 5" | 
14547 4 
eGs 12659 25 
cKo 11657 2 | I 
13 
cKs 

| 

| 


cKs 41°9 13 


42 


* The Mount Wilson velocities for No. 14847 (HD 222574) are tt, ~6,+8, +1, probably 
variable (Ap. J., 111, 221, 1950). The Radcliffe velocities are +6, 4, +2, +4, 4 6, and 
do not support much variation 

cF4 in Ap. J., 70, 207, 1929; in General Catalogu 

tcKr in Ap. 7., 111, 221, 1950; BKo in General Catalogue 

§ Lick Catalogue gives velocity for No. 18 (HD 225212) as variable. ‘The Radcliffe velocities 
are —41, —41, ~39, ~42, ~38, ~—43, the last value being the mean of two plates taken at the 
same epoch. 
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3.2. Supergiant stars.—The radial velocities of nine non-variable c-stars 
with spectral types ranging from cF3 to cKs5 were determined from 46 spectro- 
grams. ‘The velocities are given in Table III together with the Lick and 
Mt Wilson velocities, and the velocities in the General Catalogue. From the 
velocity differences we obtain the following mean values with sign: Radcliffe minus 
Lick, o'tkm/sec ; Radcliffe minus Mt Wilson*, +0°-2km/sec ; Radcliffe 
minus General Catalogue, 0-o km/sec 

‘Two of the nine supergiants have velocities listed in the General Catalogue 
as a-quality, six being b-quality, and one c-quality. ‘The limits of probable 
error and the mean probable errors corresponding to these designations in the 
Catalogue are respectively as follows:—a, +0°5; 5, +2°0, 1°25 ¢, + 5°0, 
t2°skm/se. The velocity differences Radcliffe minus Lick and Radcliffe 
minus Catalogue are satisfactorily small for the two stars with a-quality velocities, 
also for the only b- ~quality velocity for which Lick observations are available. 
For the remaining six stars, the comparison Radcliffe minus Catalogue rests 
entirely on corrected Mt Wilson velocities. ‘The Radcliffe velocities for three 
of these, Nos. 540, 7287 and 11657, lie within + 2 km/sec of the Catalogue values 
and are in satisfactory agreement within the limiting probable error for b-quality 
velocities. For the remaining three stars, Nos. 7734, 10998 and 11022, the 
velocity differences are greater than one might have expected in view of the 
generally good agreement with velocities on the Lick System. The velocity 
differences are shown in Fig. 2 together with the differences Mt Wilson 
minus Lick for eight northern supergiants observed by Joy (1, p.355) for the 
standardization of the Mt Wilson velocities of cepheid variables. 
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Fic, 2.—-Velocity differences for supergiant stars. 


3.3. Standard wave-lengths.—In the calculation of the radial velocities of 
primary standards and supergiant stars, the velocity given by each line in the 
Victoria list of standard wave-lengthst was determined without reference to the 
spectral type of each star. ‘The mean plate velocities were then obtained from 
the recommended selection of lines. In some instances, additional lines were 
either rejected or given half-weight in the calculation of the plate velocity, the 

* The discordant velocity difference for No. 14847 has been omitted in finding the mean 


difference Radcliffe minus Mt Wilson. See Table III for a note on this star. 
+ The discordant line A 4482 A has not been used. 
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criterion being principally the difficulty experienced in obtaining a satisfactory 
measurement on the line as recorded when the plate measurements were made. 
Rejections based only on the size of line residuals were avoided as far as possible. 
In the case of the supergiant stars, the criteria for selection and rejection of lines 
were the same as for the primary standards, except that rejection based on the size 
of residuals was considered advisable in some instances, particularly for lines 
suspected of dependence upon luminosity class. 


Taste lV 


Line residuals for standard giants and supergiants 
i 


Standard giants Supergiants 


Residuals Plates | Residuals Plates 


km/s + p.e. km/s + p.« No 


33 
4°07 23 
63+0°7 23 
66+06 23 


The velocity residual for each line in the wave-length table was obtained by 
finding the average, over all spectrograms for a given star, of line velocity minus 
mean plate velocity from selected lines. ‘The line residuals were then averaged 
over stars within the valid intervals of spectral type. It was found that there 
were only four lines with velocity residuals of 4km/sec or greater. The residuals 
for these lines and their probable errors are given in ‘Table IV in which is also 
listed the number of stars within the recommended intervals of spectral type and 
the total number of plates involved in the calculations. 

The largest residuals for both giants and supergiants are shown by the lines 
A 4325 and A4340. ‘These two lines are blends for which the recommended 
wave-lengths are empirical, that is, wave-lengths which are adjusted to give a 
velocity in agreement with other lines, the adopted value being a wave-length 
which reduces to zero the average residual for stars within the recommended 
interval of spectral type. For these two lines the intervals are F4-G7 and 
F4-Gs respectively. 

An examination of the Mt Wilson spectral types for the stars used in the 
Victoria work on standard wave-lengths reveals that within the interval F4-G7 
the stars are a mixture of luminosity classes, comprised of three dwarfs, two 
giants and one supergiant. Now the dependence of the effective wave-length of 
the line A 4325 A on luminosity was suggested in the Victoria work on the super 
giant « Persei, but it does not appear to have been recognized that the effectiv: 
wave-length as given has been derived from three dwarfs and two giants. Accord- 
ing to Evans (9g) the recommended wave-length for this line appears to be 
satisfactory for dwarfs. ‘The present work suggests that it is not suitable for 
giants and supergiants. It is clearly desirable that in setting up the system of 
standard wave-lengths due regard should be paid to the luminosity classes of 
the standardizing objects. 


Line 
A(A) | Stars | 
No | 
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The effective wave-lengths for the interval F4~—G5 apparently depend upon 
two dwarfs, two giants and one supergiant. ‘Iwo blended lines are given for 
this interval, A4307°994A and A4340°441A. The former wave-length is the 
weighted mean position of solar lines due principally to Fet, Catand Tit. The 
Vet line is coincident with a line at 4307-05 A due to Titi, which, presumably 
being stronger in the spectra of supergiants, will reduce the effective wave-length 
of the blend. This is confirmed by the negative sign found for the velocity 
residual for this line. 

The only blends in the Victoria wave-length table which involve lines used 
in spectroscopic absolute magnitude work are A4215 A containing lines due to 
Fet and Sri, and A4340A containing Hy. ‘The latter is, however, only used 
for absolute magnitudes for K and M type giants, and is recommended for radial 
velocities only for types F4~G5. ‘The velocity residual for A4215 A was found 
to be satisfactorily small for both giant and supergiant types in spite of the fact 
that the intensity ratio AA 4215/5250 is a well-established measure of luminosity. 
The line A4314A due principally to Fe1t, Secu and Tilt was suspected of 
dependence upon luminosity in the Victoria work on « Persei, but the average 
residual for this line for stars in the range from cF4 to cK8 has been found in the 
present investigation to be small. ‘The residuals for individual stars in the range 
may, however, be quite large. 

These considerations suggest that for supergiant spectra there is some 
justification for rejecting measures of the lines A 4307 A and A 4314A which differ 
markedly from the plate velocity. For these lines the wave-lengths are the 
weighted mean positions of the blended components, but the lines A 4325 A and 
A 4340 A have wave-lengths which are adjusted to give the correct velocity. It was 
considered impracticable to make a satisfactory adjustment of these wave-lengths 
for cepheid variables in view of the fact that the line spectrum is abnormal (10) 
on the ascending branch of the light curves although closely corresponding to a 
supergiant of luminosity class Ib at minimum light. It was decided to discard 
the lines AA 4325, 4340 rather than attempt to make further adjustments to their 
effective wave-lengths since such wave-lengths, even if correctly adjusted, cannot 
add much weight to a plate velocity, although they might be useful for normal 
spectra when some of the calibrating lines are for some reason not accessible. 


TABLE V 


Effective wave-lengths used for the determination of 
radial velocities of cepheid variables 


| Wave-length | Spectral Wave-length | Spectral 
| in angstroms interval in angstréms interval 


3901°537 4215686 *4-Ko 

4005 244 4254348 
4030°045 4260°4588 
4934°494 43°7°994 
4045 °705 "4-K 4314°6053 
4071°735 4404°745 
4092°478 4415°137 
4132°230 4461810 


J 
4 
: 


No. 4,1955  cepheid variable stars in the Southern Hemisphere 371 


4. Cepheid variables 

The programme of observations of cepheid variables was carried out with 
due attention to the need for securing spectrograms at phases well distributed 
over the light curves. In order to facilitate the computation of the phases, and 
hence the selection of stars for observation, tables were constructed which gave, 
for each star on the programme, thirty-one multiples of the reciprocal period. 
The phase of the light variation on any night in a given month could then be 
written down at once by adding the phase shift for the night in question to the 
phase at midnight on the last night of the preceding month. ‘These predicted 
phases were used to estimate the photographic magnitudes from the light curves, 
and hence the exposure times required. ‘The light elements given in the Russian 
Catalogue of Variable Stars(11) were used, Elements for the star CPD — 60°4739, 
which is not contained in the Catalogue, were kindly supplied by the Riverview 
College Observatory. 

The observing programme was limited to cepheids brighter than photographic 
magnitude 10™"*5 at minimum light. Fifty-five cepheids were observed on 112 
nights, a total of 830 spectra being secured. ‘The published magnitudes were 
found to be not entirely dependable, and the estimation of exposure times was 
frequently more a matter of judgment than of calculation. In some cases, the 
ephemerides were not accurate enough for the prediction of phases, notably 
for the large-amplitude cepheids ‘T Antliae, U Carinae, VY Carinae, AQ Puppis 
and RZVelorum. In these cases provisional corrections to the computed 
phases were made with the help of rough velocity curves deduced from plate 
measurements made during the observing programme. 

4.1. Radial velocities and velocity curves. —The radial velocities of the cepheids 
were determined from the effective wave-lengths given in ‘Table V, using approxi- 
mate spectral types for each phase. ‘The spectral types given in the Russian 
Catalogue, mostly from Harvard classifications, were used. In cases where the 
data are insufficient, spectral types at maximum and minimum light were taken 
from the period-spectrum relations of Code (12), the approximate types at 
intermediate phases being then obtained from the range in type and the phase, 
‘The radial velocities are given in Table VI in which the following data are 
tabulated: (i) the Radcliffe plate number with the designation C for Cassegrain 
and c for the camera used, (ii) the Julian Date of observation, (iii) the geocentric 


Taste VI 
Radial velocity observations of cepheid variables 


Epoch Phase Velocity Plate Epoch Phase Velocity 
JD P km/sec Ce ID P km/sec 
T Ant JD 2434000 V496 Aql JD 2433000 

122°28 249 871°44 12°7 
130°25 ‘12 36°3 261 875°43 40 64 
133°26 63 276 ‘28 4°3 
137°23 "30 +47" 290 897°43 63 18-2 
145°24 66 300 900'44 ‘ 26 
167°25 308 61 
174°20 "57 319 06°33 
197°19 “45 338 913°37 o"7 
200°20 0°96 352 g28°41 


Plate 
694 
720 
73° 
744 
763 
835 
918 
93° 
26 
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Epoch Phase Velocity Epoch Phase Velocity 
Jb P kin/ sec JD P km/sec 


V496 Aql (continued) V Car (continued) 
937°3° + 988-58 
943°27 996°53 
955°24 995°53 
g62°23 - 1000°54 
965°25 10° 1018°46 
1157°66 L 6 1020°56 
1183 °61 1022°50 
1197°61 ‘ 1025°49 
1219°50 - 1027°58 

1223°4! ° 1029°54 
iCar JD 2434000 

1098 °27 

1110°28 

1113°23 


ee su on 


Y Car JD 2434000 
57°5° 
85°49 78 
88°49 
92°50 

110°39 62 
122°32 ‘go 
130°33 "10 
137°26 ‘or 
145°27 ‘21 
191°27 


owns wus 


nN 
“I @ 


ER Car JD 2433000 

$22°28 0°78 

827°26 42 

829°26 68 

831°24 04 

834°28 ‘33 

840°28 

844°30 63 

1057°56 26 

1077°58 “BO 

1085 °59 ‘go 

1092°54 ‘Bo 

34 . 6 1102°48 08 

40°34 1117°40 02 

3°35 1122°35 66 

29 11§7°32 "19 

5°35 ' 1163°33 ‘97 

4°32 03 + 1183°29 o's6 

116730 Gli Car JD 2433000 

1174°30 +25" 822°30 
V Car JD 2433000 827°29 
973 °60 829°29 
984°59 6° 831°26 


Ce 
359 +7 
37° 
383 
391 + 3 
be +24 
920 + © 
977 
1012 
499 TIS 
504 
$12 
516 
530 -36 5°3 
539 44-44 ‘73 + 547 
543 53°38 98 1°5 602 
546 57°45 09 18°5 610 
575 77°40 65 626 4 
585 80°39 ‘74 669 
603 85°38 ‘88 +15°7 694 
| 671 110°41 58 720 
681 113°32 66 + 16°38 744 
685 117°34 78 }22°7 763 
862 175°16 “41 - 899 
874 
—32 
179 
185 
192 17 
577 6 
587 
604 
720 20 
73° 
744 
764 re 
825 
$38 
860 
440 
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Plate Epoch 

Ce JD 
Gl Car 
834°30 
840°30 
1080°54 
1113°39 
1117°39 
1122°35 
1133°39 
1145°33 
1157°35 
1167°31 
1174°33 
1191°30 
1200°30 


IT Car 
133°: 
137°: 
157°: 
163° 
174° 
183° 
191° 
197° 


UX Car 
815°23 
829° 
831° 
534° 
844° 
1077° 
1035" 
1104": 
1122" 
1130” 
1137" 
1145" 
1167" 
1191" 
1197°22 


VY Car 


818-26 
$27°24 


$29°24 
831°23 


$40°26 


Phase 
P 


(continued) 
0-99 
"34 
"57 
“BS 


“15 
‘74 


JD 2434000 
0°56 
‘73 
‘39 
19 
64 
‘84 
“69 
‘09 


JD 2433000 
0°35 
‘43 
"79 
"34 
‘88 


JD 2433000 
0°30 
‘88 
0°46 


Velocity 


km/sec 


20°3 
32°1 
16°7 
24°3 
14°4 
22°6 
21°2 
27°4 
12°6 
31°3 
12°9 
18-2 
16°0 


Plate 


Ce 


227 
547 
577 
588 
626 
675 
633 
744 
763 
860 


Epoch 
JD 
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Phase 
P 


VY Car (continued) 


844 °27 
10§7°52 
1077'54 
1080°48 
1092°52 
1113°34 
1117°32 
1137°27 
1145°25 
1174°25 
1197°23 


V Cen 
827°40 
$29°36 
531°32 
334°37 
340°36 
544°39 
871 
375° 
SS1- 
1088 
1101 
1137 
1157 
1191" 
1226 
1233 


AZ Cen 
o8o0'56 
113°42 
117°43 
122°36 
130°36 
133°4! 
137°32 
145°30 
163°40 
181°29 
197°27 


UZ Cen 
113°46 
117°§2 
122°40 
145°35 
174°35 
151°24 
200°43 
206°25 
219'23 


0°67 
‘92 
‘97 
13 
07 
‘12 
"54 
07 
0°28 


JD 2433000 


o'l4 
49 
‘Bs 
‘40 
49 
23 
‘tl 
93 
‘58 
23 
“58 
"39 
72 


JD 2434000 
21 
40 
‘99 
“43 
65 
‘77 
‘34 


JD 2434000 
"26 
73 
62 
‘31 
35 
"190 
0'”77 


373 


Velocity 
km/sec 


28°6 
10°2 
2°3 
33°! 
15°6 
13°9 
26°38 
16°9 
3°4 


‘50 
19 
‘90 
35°9 
678 19°7 186 12°83 
73° 61 193 30°7 
744 22°5 202 16°8 
806 8-9 214 11°2 
825 9°1 228 
860 22°8 246 jo"o 
886 15°3 259 44°60 
19'0 275 39°7 
899 611 2°2 
gi8 19°6 645 32°0 
930 200 9'4 746 5°3 
9045 206°25 19°4 16°7 
950 207°23 17°6 888 ars 
goo 13'0 
989 18-4 
155 1022 316 
167 + 20°2 
179 
18s 588 16°5 
192 + 78 679 10°7 
216 "59 20°6 8. 
194 
227 “41 +10 
721 6's 
576 “76 730 
4 744 
) "10 2 i 
+m 99 4 918 o's 
763 "16 1 
825 "04 
$38 I 
899 65 + 26 679 
918 0°27 + 6 654 
694 r 
765 8-8 
155 O73 860 25°5 
179 17'S 886 20°7 
185 | 31'S 930 
192 27°53 945 6-2 
213 | TI3°2 960 4°2 
26* 


Epoch 
JD 


XX Ce 
$25" 


426 


$27" 


428 
429 
530 
$41 
$32 
433 


534° 


1122 


1130" 


1133 
1227 


n 
15 
21 
15 
19 
18 
15 
15 
18 
15 
18 
‘21 
25 
25 
27 


V339 Cen 


818 

$27 
1153 
1117 
1122 
1130 
1133 
1137 
1145 
1187 
1163 
1167 
1151 
1183 
1197 


39 
35 
64 
65 
62 
60 
5! 
56 
53 
49 
46 
“44 
‘41 


V378 Cen 


130 
133 
137 
145 
157 
167 
174 
rst 
183 


223° 


231 


56 
‘53 
40 

53 

48 
“41 
41 
40 
43 
29 


V381 Cen 


825 
826 
827 
825 
S29 
S31 
534 


1130" 


1133 
1137 


"19 
19 

19 
19 
19 
19 
20 
19 
60 
57 
“48 
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Phase 


JD 2433000 
04 


0°55 


JD 2433000 
0°87 
82 
‘06 
‘Ol 
“16 
‘58 
‘43 
‘69 
42 
‘74 
"22 
‘43 
o’9! 


JD 2434000 
‘42 
67 
“$2 
‘06 
14 
22 
‘54 
7 
JD 2433000 
19 
‘39 
"59 
77 
02 
71 
0°48 


Velocity 
km/sec 


24°7 
23°3 
25'1 


US 


260° 


19° 


2 
37°53 


16°9 
16°9 


Vol. 115 


Plate Epoch Phase Velocity 
Ce JD P km/sec 
V381 Cen (continued) 
919 1197°39 0°28 37°53 
93° 1200°41 0°88 39°6 
V419 Cen JD 2434000 
543 57°61 14°5 
577 77°63 ‘94 
538 80°62 ‘43 25°0 
628 92°57 ‘65 24°9 
684 117°47 ‘17 14°2 
694 122°38 ‘06 10°2 
721 130°39 "52 19°9 
73! 133°39 ‘07 
745 137°34 “78 
899 191°32 "$7 
gi8 197°29 66 22°5 
950 207°26 0°46 24°5 
CPD -—60° 4739 JD 2434000 
549 57°63 13°4 
560 62°64 ‘BI m9 
5387 80°65 ‘Ol 21°5 
604 85°64 “go 16°9 
609 58-62 43 12°8 
627 92°65 “15 15°35 
647 10164 ‘75 
653 104°60 14°! 
672 110°42 31 14°58 
686 117°61 59 4'2 
695 122°63 45 
807 157°50 ‘68 61 
539 167°44 0°45 9°3 
RCru JD 2433000 
155 7°3 
180 527°33 05 25°9 
186 $29°33 "39 13°5 
193 831°29 73 5°7 
202 534°34 25 23°7 
214 340°34 28 
227 344°34 97 30°9 
602 1085 °63 39 14°5 
610 108859 ‘go 7°6 
645 1101°§2 ‘12 29°1 
679 1113°54 "18 23°9 
695 1122°48 “71 
731 1133°44 ‘60 4°53 
765 1145°47 0°66 
Cru JD 2433000 
156 $18°33 0°70 4°4 
180 827°35 62 
186 829°35 05 23°"! 
193 S313! ‘47 
214 840°35 "40 
225 $44°358 25 15°4 
647 110162 27 


374 | 
176 = 
179 
ft 132 12 30°3 
185 ‘22 19°! 
189 
192 ‘40 11°6 
195 49 2°5 
200 58 
201 ‘67 68 
695 ‘00 26°1 
722 "72 17°2 
731 ‘00 33°3 
156 
180 
680 
686 
695 
722 
732 
is 740 
766 
307 
526 
439 
888 33°2 
893 
721 24°60 
731 10°7 
745 22°6 
765 
807 
839 15°3 
861 23°2 
SSS 24°2 
593 13°0 
975 12°7 
175 49°6 
176 38°7 
179 284 
182 13"! 
Ss 19’! 
201 
274 49°9 
722 
730 
740 —23°5 
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Plate Epoch Phase Velocity Plate Epoch Phase Velocity 
Ce JD Pp km/sec Ce ID P km/sec 
S Cru (continued) AG Cru (continued) 

731 1133°5I 899 191°34 0°96 22°6 
745 1137°44 ‘74 930 200°40 "32 6°3 
807 1157°45 ‘Ol g6o 219°30 25 

839 1167°39 "13 975 223°25 0°28 
893 118341 ‘54 
oss 
990°44 
T Cru JD 2433000 se 43 > 
818-34 
995° ‘ol 
9938 ‘21 
834°33 
840°32 
844°33 1015 ‘24 
1085-61 24 
1088°57 44 
1092°59 1022 
1101°50 1025 ‘96 
sine de 1027 16 
1104°58 1027°33 
1113°53 ‘37 
1122°46 
. 1033": "77 
1145°460 1033" 
1080 "$3 
1092" "76 
1092" 


o 


ONO 
< 


WN 


o-< 
nw 


x 


~ 


X Cru JD 2434000 
122°57 
130°52 "13 
133°45 "60 
137°42 "24 
145°5° 54 
157°43 "46 
167°37 
174°35 15 
151°35 
133°35 63 
197°360 “BS 
227'22 0°65 


0 


R Mus JD 2433000 
818°; 0°95 
$27": "14 
‘41 
S41” ‘67 
‘08 
840 “BR 
844° ‘41 
1088" 
66 


AG Cru JD 2434000 
101°58 $7 
70 1137" 43 
17°58 "4 1183°3: 
122°52 03 
130°48 ‘10 ’ S Mus JD 2434000 
133°40 55°56 0°08 
137°40 3° 66 
145°45 137 13 
1§7°40 145 ‘97 
183°36 206 0°27 


741 
73! 27 
745 35 
765 28 156 
807 33 180 — 56 
839 20 186 
860 35 193 t 27°3 
886 34 202 6-4 
919 12 228 + 3°5 
989 21 610 ~ 9°5 
646 +17°9 
628 + 27°2 
640 r 32 
654 r 40 
6905 
721 
73! 
745 + 28-0 
705 
$07 2°2 
835 24°4 
886 
893 t+ 66 


Plate Epoch 
Ce Jb P 
RT Mus JD 2434000 
695 122'44 0°26 
721 130°45 86 
745 137°36 ‘10 
795 145°4! 71 
$07 157°37 ‘58 
586 ‘34 
891 183°31 ‘99 
197°32 53 
930 200°36 51 
959° 207 °32 ‘77 
975 223°21 0'92 
S Nor JD 2433000 
247 71°26 0°83 
260 $75°25 ‘24 
275 $81°27 85 
248 897 °22 49 
298 00°20 ‘79 
307 10 
317 go6'19 ‘41 
328 
339 
3452 g25°21 ‘66 
359 93420 ‘28 
37° 937°22 ‘79 
$87 1181-52 “64 
gol 65 
931 1200°44 58 
975 1223°33 "92 
1009 1231°34 ‘75 
1022 1233°24 94 
1022 1233°25 0°04 
BF Oph JD 2433000 
259 $75°28 
274 "42 
288 807°24 34 
298 Q00'22 ‘o8 
308 Bi 
317 55 
338 913°24 28 
352 928°22 96 
333 943°21 64 
893 1153°54 ‘72 
goo 1191°46 67 
919 1197°40 ‘14 
945 1206°33 32 
g60 1219°47 0°56 


«Pav JD 2433000 
245 571°39 0°48 
262 S75 ‘41 ‘92 


881°40 


Phase Velocity 
km/sec 


Plate 
Ce 


259 
297 


497 
506 
517 
§2!1 
53! 
544 
570 
536 


720 
73° 
763 
Kob 
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Vol. 115 
Epoch Phase Velocity 
JD P km/sec 
« Pav (continued) 
$97°39 0°34 31°0 
Yoo’ 32 "66 45°2 
9°3°33 ‘99 17°9 
906°35 "33 35°7 
913°40 ‘10 26°4 
925°33 ‘75 53°2 
934°32 ‘4! 
955°26 §2°5 
959°26 "15 25°! 
962°25 43 43°5 
965°27 
967°24 0'03 24°5 
AP Pup JD 2433000 
93401 25°2 
937°61 "14 4°4 
943°61 32 
955°59 ‘87 32°5 
62°55 o's 
965 °S4 63 23°38 
967 °53 02 
g65°51 22 
973°50 ‘20 1°2 
976°50 ‘79 36°3 
934°51 "36 10°3 
986'50 76 34°7 
988°51 15 
992°54 ‘94 22°0 
994°50° "33 
996°44 ‘71 
9938°46 13 
1000°42 ‘49 16°83 
03 7°s 
1020°37 2 11°3 
1022°39 82 32°9 
1110°26 10 2°8 
1163°17 19°3 
AQ Pup JD 2433000 
18°38 31°3 
22°43 ‘35 
27°40 
29°47 ‘58 
31°42 "65 62°0 
53°42 35 
77°32 “15 *37°9 
80°28 ‘28 6 


‘95 
05 
‘45 
‘Bs 


376 
309 
+14°2 318 
12°4 340 
4 352 
19°2 360 
73 398 
409 
14°9 423 
432 
+ 9°5 
= 
10°5 362 
374 
+-23°6 333 
22 394 
10°0 412 
+22°2 422 
3°5 43! 
12°9 435 
9°1 440 
+22°6 443 
+18°0 446 
+184 452 
462 
19's 468 
+ 3°5 474 
454 
497 
500 
~47°5 $00 
23°4 669 
52°6 
21°6 
16°5 
15°2 
11 6 
40"! 
40°58 
130°19 
145°20 39°9 


No. 4, 1955 cepheid variable stars in the Southern Hemisphere 377 


Plate Epoch Phase Velocity Plate Epoch Phase Velocity 
Ce JD P km/sec Ce JD P km/sec 
AT Pup JD 2433000 W Ser JD 2433000 
965°59 248 871°; 0°77 3°3 
973°59 ‘61 09 
976°57 06 
986°53 
988°53 ‘Bs 
992°56 ‘46 
996°46 04 
995 °48 ‘35 
1000°44 “64 
1022°46 ‘04 
1029°36 ‘98 
1057°34 
1062°33 
107735 
1080°32 
1088-30 
1130°22 "12 
1163°18 
1167°19 66 
1174°24 ‘72 
1181°18 0°76 


Ww WN 
nes 
wm 


897 19 42°! 
goo": 58 

903°: 97 

906": ‘37 

gio” 

913 

937° 

1191" 

1197°58 


wu w 
- n+ 


w 
on 


X Ser JD 2433000 
0°36 
875°3: ‘93 
goo 
903° 
gob 
913°33 
943°23 

WX Pup JD 2434000 

1206°39 

1219°41 

1223°35 

‘7 1231°35 

‘20 

65 

“66 

‘66 

‘99 

gl 
0°03 


NN SUVS & 


RV Sco J 
527°460 
429°50 
$31°37 
534°43 
871°27 
$75°26 
881°28 
597 °23 
993 
906 


Aauu acd 


JD 2433000 
4! 0°57 
5°40 “16 
39 "05 
32 
‘Bs 


Aw 


913 
1197 
1223 
1227 


WS 


aS 


RY Sco JD 2433000 

575 

51°36 
597 

goo’28 
03°28 
906°27 
913°26 


~ 


D 


s 
506 " 
517 25 
3 
558 
576 
608 | 
626 
645 
678 35 : 
806 9 
838 33 
838 
259 77 2 
U Sg 274 76 I 
- 308 38 23 
276 5 317 87 
288 87 328 
298 got a2 
308 93°32 30 5 
317 42 “74 12 989 
338 913°36 79 
359 934'°25 ‘89 24 
37° 937 °27 33 5 
353 943°25 ‘22 261 
955°22 19 276 9°6 
408 g62°22 4 258 
= 95°22 “45 7 295 
888 1151 55 "$5 308 
93 317 32°6 
ras 335 —21°9 
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Epoch Phase Velocity Plate Epoch Phase Velocity 
JD P km/sec Ce J.D. P km/sec 
RY Sco (continued) V636 Sco (continued) 
925°26 33°60 1197°460 
934'25 ‘59 1200°46 23 
1197°55 0°56 1206°35 ‘10 
1219°39 ‘OI 
V482 Sco JD 2433000 1223°34 o’6o 
S71°31 5 
75°30 RTrA JD 2433000 
03 827°41 0°48 
897°27 “56 $29" 06 
goo ‘21 $31" 
903 ‘87 534 ‘54 
4 871 ‘40 
928 875 ‘58 
937 7 1130°6: ‘95 
1130" 1133°67 ‘83 
1137" ‘ 1137" ‘99 
1133 1197 65 
1227°32 1233 0°20 
1234°15 
1234°19 


SK HK BHU 


S'TrA JD 2433000 
$29°3 
Vs00 Sco JD 2434000 831°, "32 
133°06 0'04 3° 834": ‘Bo 
137°63 47 2° 54442 ‘39 
145°606 43 5° 871°21 63 
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54 
96 19°7 
‘Bo + 
‘06 +25°0 
68 38°5 
‘57 45°2 
o’61 + 


phase of the light variation at the time of mid-exposure, calculated from the 
ephemerides in the Russian Catalogue, and (iv) the radial velocity. 
are arranged in alphabetical order and the observations for a given cepheid in 


chronological order. 
next to the star name. 
epoch column. 


VII 


Observational data for cepheid variables in the Southern Hemisphere 


Star 


Ant 
V496 Aql 
/ Car 
Car 
Car 
Car 
Car 
Car 
Car 
Car 
Car 
Cen 
AZ Cen 

UZ Cen 

XX Cen 

V339 Cen 
V378 Cen 

V381 Cen 
V4i19 Cen 
CPD —60" 4739 
R Cru 

Cru 

T Cru 

X Cru 


Myg 
Bo 
9g’! 
8-7 
9°4 
9°4 
g'orto'l 
g°6- 
9°3 
7'5 
7°4- 
BO 
9°4-10°0 


Period 
(days) 


§ "9071 
63069 
35°532 
35°749 
6°6967 
3°6308 
7°7178 
4°4306 
7°5356 
30323 
18-961 
4940 
3°2127 
3°3343 
10°950 
9°406 
6°45387 
§°0789 
§°§103 
5 6218 
53256 
46900 
6°7332 
6°2199 


Spectral 


type 


(52 


Normal 

velocity 
km/sec 
29°5 
5°2 
14 
m1 
21°4 
14°0 
10°4 
22°3 
15°7 
24°3 
17°7 
30°8 
17°2 
12°58 
13°5 
6°6 
—25°0 


Velocity 


range 
km/sec 


45 
14 


The cepheids 


The thousands part of the epoch in Julian days is given 
To this should be added the number of days given in the 


No. of 
plates 


380 
4°°3 586 
14°0 608 
2°8 626 
ae 642 
669 
+ 26°7 678 
10°0 683 
763 
$06 
9°7 825 
886 
10 899 
| 
= 9 
F8-Ko 40 21 
Fo-K3 52 16 
Go-K2 32 21 
Fs 24 10 
FR 4 18 
20 17 
F7-K2 18 13 
F4-Gs 34 19 
F9-Ko 59 16 
Fs5-Gs5 39 17 
F3-Go 19 
F2-G5 38 9 
F6-Ko 35 14 
Fs5-Gs 17 I 
F6-G7 39 13 
Ko 17 12 
F8 16 13 
F6-G7 34 14 
Fs5-Gs 36 14 
Gi-Gs 34 17 
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. Period Spectral Normal Velocity No. of 
(days) type velocity plates 
km/sec 
AG Cru 3°3374 Fs 4’ 
Dor 5° 9°8426 F2-F9 8: 
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WX Pup ‘9393 Go 53° ‘ 12 
U Ser "7449 F7-Gs5 17 

W Ser $947 F2-G6 

X Ser 0122 Fs5-—Go 1s 
RV Sco 0613 Fs5-G5 15 
RY Sco "314 F9-G7 ‘ 10 
V482 Sco 5278 Fs i8 
3163 ove 12 
7963 (55 21 
3893 F6-G4 3! 14 
3236 F6—G8 16 
5683 F4-G2 14 
6398 20 
"3710 17 
2272 F8p 17 24 
'§925 27 15 
‘9235 G5 30 18 
"396 G8 57 7 
097 F8-—Ko 49 10 
'§499 32 1§ 


19 
14 


Vs500 Sco 
V636 Sco 
R‘TrA 
TrA 

TrA 

T Vel 
Vel 


CN 


‘The observations were plotted according to phase, and schematic or smooth 
curves, whichever seemed appropriate, were drawn through the velocity points. 
The radial velocity curves are given in Figs. 3, 4, 5,6 and 7. Normal velocities 
were found from the curves by the method of equal areas. Using a planimeter, 
the residual area was found for figure of eight loops for three assumed, and 
approximately correct, values of the normal velocity, the velocity for zero residual 
area being then found graphically to the nearest o-1 km/sec. ‘The uncertainty 
in the normal velocity depends mostly upon how well the velocity curve is defined 
by the observations. ‘The estimated average probable error of a single normal 
velocity is of the order of 1km/sec. ‘The normal velocities are given in Table VII 
with other relevant data. 

4.2. Comparison with Lick.—Five cepheid variables for which definitive 
velocity curves had been obtained by the Lick Observatory from high dispersion 
spectrograms, taken at a much earlier epoch, were included in the programme 
in order to check the constancy of the normal velocities and, if possible, the shape 
of the velocity curves. ‘The stars in question are the five brightest cepheids in 
the Southern Hemisphere, each having noteworthy features : / Carinae for its long 
period, 35°5 days, and for suspected irregularities in the velocity curve (13); 
8 Doradus for its high galactic latitude, —32°-3, and the unusual shape of its 
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velocity curve (14); « Pavonis for its high galactic latitude, — 26-3, and its 
recently proposed identification (15) as a type II cepheid; W Sagittarii for the well- 
established variability in the shape of its velocity curve (16); and X Sagittarii 
for the early epoch, 1904-1905, of the one and only set of observations (17) made 
for its radial velocity. 
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Taste VIII 


Comparison with Lick 


Lick | R—L 


km SCC km sec 


Star Radcliffe 
km sec 
/ Car 
Dor 
« Pav 
W Ser 
X Ser 


The normal velocities for these stars are given in Table VIII] with the Lick 
velocities for comparison. ‘The velocity difference for / Carinae is probably 
significant in view of the fact that Jacobsen (13, p. 261) obtained evidence for a 
small variation in the normal velocity from separate solutions for two epochs. 
For B Doradus, the normal velocity differs inappreciably from the value previously 
found by Applegate (14), and the remarkable shape of the velocity curve in Fig. 3 
was found to be practically identical with the Lick curve from high dispersion 
spectrograms taken in 1920. ‘The present observations support Applegate's 
conclusion that the normal velocity of 8 Doradus is constant, and provide adequate 
evidence that the shape of the velocity curve is still much the same as it was some 
thirty years ago. For W Sagittarii, the hump in the velocity curve at phase 
o”-4in Fig. 6 is probably real. It coincides with the phase at which Jacobsen 
(16) found a similar feature which was considered to be variable. ‘The hump 
is slightly more pronounced in the present observations but the normal velocity 
is practically the same. No unusual features were noted about the velocity 
curves for « Pavonis and X Sagittarii. For the latter, the scatter in the observa 
tions is larger than for the other stars, possibly on account of the poor quality of 
the lines in the spectrum. 
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4.3. Comparison with Mt Wilson.—The cepheid programme included five 
variables between declinations — 19° and — 34° for which Joy obtained schematic 
velocity curves. The velocities for these stars are listed as b-quality in the 
Mt Wilson General Catalogue. The normal velocities and the number of observa- 
tions used in their determination are given in Table [X. For the purpose of 
this discussion, Joy’s phases have been recalculated with the ephemerides used in 
the present work. ‘This results in a small shift in the velocity curves without 
changing the normal velocities. 

The velocity differences Radcliffe minus Mt Wilson for BF Oph, RV Sco and 
RY Sco are satisfactorily small. The two series of observations for BF Oph 
and RV Sco combine to give well-defined composite curves, but for RY Sco a 
reduction of o” 18 in the computed phases of the Radcliffe observations is 
necessary for coincidence between the two curves. Apparently the ephemeris 
for RY Sco, which contains a squared term, predicts the light maximum about 
3°7 days too early at the present epoch. 

The velocity differences for WX Pup and U Sgr are larger than one might 
expect for b-quality velocities for which the maximum probable error is 
+ 2°0km/sec. A comparison of the individual observations in both series is 
necessary before it can be decided if there have been any significant changes in 
the velocity curves and the normal velocities. Joy’s velocity curves for these 
two stars are given in Fig. 8 with the Mt Wilson and Radcliffe observations. 
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‘The phases of the Radcliffe velocities for WX Pup have been increased by 0°16 
in order to bring the two series of observations into coincidence on the descending 
branch of the velocity curve. Evidently the ephemeris for WX Pup is not 
accurate enough for the present epoch on account of an uncertain period given 
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only to four decimal places. The present shape of the velocity curve is definitely 
not the same as Joy's curve which, in the interval of phase between o*-15 and 
os, lies about 10km/sec below the Radcliffe observations. The difference in 
the shape of the velocity variation at the two epochs depends almost entirely 
upon three of Joy’s observations with dispersion 120A/mm for two phases 
within the interval in question. It would appear that there is not much evidence 
to suggest that there has been any appreciable change in the normal velocity 
or the shape of the velocity curve for WX Pup in the interval of 17 years between 
the mean epoch of the two series of observations. ‘Turning now to the case of 
U Sgr, it is evident from Fig. 8 that Joy’s velocities combine satisfactorily with the 
Radcliffe velocities. Joy's schematic velocity curve gives a very satisfactory 
fit to the Mt Wilson observations, but it is clearly inadequate for the combined 
observations. ‘The difference is greatest in the interval of phase from o” 2 to 
o”-6 in which there are no Mt Wilson velocities. ‘The Radcliffe observations, by 
contributing seven velocities in the same interval, confirm the evidence from the 
shape of the light curve that a schematic curve is not likely to give a satisfactory 
representation of the velocity variation. 

Conclusion.—\t may be concluded from the discussion on the standardization 
of the velocities that the observations described in the paper have given a homo- 
geneous set of velocity determinations in a system which differs inappreciably 
from the Lick System of radial velocities. ‘The normal velocities of cepheid 
variables presented in the paper may be introduced with confidence into other 
radial velocity data, such as Joy’s data for cepheid variables in the northern 
hemisphere, which also differ inappreciably from the Lick System, 
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BRIGHT DOUBLE-LINED ECLIPSING VARIABLI 
Sergei Gaposchkin 
(Communicated by Cecilia Payne-Gaposchkin) 


(Received 1955 March 21) 


Summary 


The star HD 77464 found by Feast to be a spectroscopic binary is 
shown by photometric observations to be an eclipsing variable with the 


same period as that established by Feast The essential elements are 
Min = 2434091 "302 + 6% Max 6™-61 m,=7™-25; 
-o8; r, 0'126 4°260; re 0'123 = 4'160; 5618 
Me 5'4180; Mogi 1-78, V ge 1-67, 
Ve ™-45; Sp, Sp,=B3-2 


The dimensions of this new system may support the idea that there exists 
a group of progressively lower mass all along and below the main sequences 
beginning with the Wolf Rayet stars and ending at the late A spectral class, 


1. The bright star HD 77464 (8" 57™; 51 10’, 1900) of spectral class Bz 
(Harvard) has been announced by M. W. [east (1) as a spectroscopic binary 
with a period of about seven days and with similar components, Although the 


spectroscopic elements have characteristics which do not suggest that the orbit 
may be so inclined as to produce alternating eclipses as seen from the Earth, | have 
examined the star's brightness. 


2. From inspection of several hundred plates in the Harvard Collection, 
| have estimated the brightness of the star. ‘his photographic material suffers 
from heterogeneity; the location of the star in a thick nebulous region of the 
Milky Way, and the relatively long exposures of the plates, produced a halo around 
the image, but in spite of these difficulties the variation of the star's brightness 
is easily perceived. 


3. The comparison stars are collected in ‘Table I. 
I 
Comparison stars for HD 77464 


TABLE 


Hb Pg mag Steps Source of mag. 


77653 5°32 15°0 (2) 


76230 6.35 7°38 (2) 

76968 7°23 2‘0* (2) 

77305 7°52 oof (3) 
* Definitely variable + Has a faint companion. 


All my observations are given in ‘lable II] expressed in steps. ‘The relation 
between the steps and the magnitudes is expressed in the form 


Ipg = 6-30 + of 15(8 — x), 
where x is the number of steps. 

4. Feast’s period seems to be correct. | have observed few minima but 
observations at the beginning and end of the Harvard series fit the predicted 
minima. For the middle of the interval | have many plates taken with AX, AY, 
and AK telescopes, which, however, have appreciable differences in colour 
sensitivity. This factor, together with the brightness of the variable and the 
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Il 
Photographic observations of HD 77464 


JD Ini Jb Units JD Units 
2 400 000 + 2 400 000 + 2 400 000 


d 
14979304 17 342°534 20 620 
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Il (continued) 
JD Units JD Units JD 
2 400 000 2 400 000 | 2 400 000 + 
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surrounding galactic absorption, makes the plates of little value. Nevertheless, 
they seem to indicate that the observed minima in this intermediate interval 
deviate from the predicted ones. ‘This observation, if real, would suggest the 
existence of a third body. 

5. Fig. 1 represents my observations reduced according to the period given by 
Feast. ‘The maximal brightness, the depths of the minima and the duration of 
the eclipse were first expressed in steps and then transformed into magnitudes. 
The narrowness of the minima, the limited number of the plates, and the relative 
length of the period account for the fact that few observations fall in the interval 
of both minima. Nevertheless, if the photometric scale is correct, the light curve 
is quite sufficient to obtain important dimensions of the system. ‘The elements 
of the light curve are : 

Min = JD 2434091% 302 + 6°8g2« 
Max = 6-61 

= 7-28 

mM, 

D =o" -08 (computed) 
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Comparison with the spectroscopic orbit reveals that the slightly heavier 
and larger component is eclipsed during the deeper minimum (phase ~o*:s). 

6. The relative elements are the following : 

Ratio of the radii of the components 0°98 

Ratio of their surface brightnesses 1°034 

Ratio of their luminosities 1°10 

Inclination of the orbit 88° 14' 

Relative radius of the larger component 0°126 

Relative radius of the smaller component 0°123 


Fic. t.——Light curve (above) and radial velocities (below) of HD 77464 (HV 11109). 


The abscissae are phases in terms of the period The ordinates are magnitudes (above) and 
velocities in km/sec (helow) The dotted line refers to the velocity of the system The crosses 
represent the intersection of the orbital velocities and phases of predicted minima 


By using some of the spectroscopic results of Feast, the minimal masses and 
the distance between the centres, the final absolute dimensions of this system 
may be obtained. Assuming the Sun's radius to be 0-6951°10° km, then 

Radius of the larger component is 4°26 
Radius of the smaller component is 4°16 
Mass of the larger component is 5°618, 
Mass of the smaller component is 5°418 


Bite 
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If the spectrum of the system during its maximal brightness is Bz, then the 
individual spectra of the components may be computed from the observed 
difference of their surface brightnesses, as B2-8 for the larger component, and 
B3-2 for the smaller. 

And finally, taking the corresponding temperatures (4) as 17530 deg. and 
17 090 deg., their absolute photographic and visual magnitudes may be computed : 


Absolute photographic magnitude of the larger component 
Absolute photographic magnitude of the smaller component 
\bsolute visual magnitude of the larger component 
Absolute visual magnitude of the smaller component 


6. Conclusions..-HV 11109 1s interesting in the definiteness of its dimensions, 
which is the result of the two deep minima and the lack of a pronounced ellipticity. 
There are some ten systems having approximately the same spectral class but 
none of them has both of these advantages. 

The new system is remarkable in that its masses and to some extent the radii 
are smaller than the majority for the same spectral class. If Feast's spectral 
classification, B2, is accepted, the low masses will be even more striking. Perhaps 
this system demonstrates the existence of a group of low mass all along the main 


sequence whose upper end is formed by the Wolf Rayet stars (about ten solar 


masses), and whose lower part goes through the system of UX Monocerotis with 
ts spectral class of A7 and its mass below that of the Sun. 

When this paper was in proof, the writer's attention was called to the dis- 
covery of the star’s variability by van Houten (§), who, however, did not give the 
correct period 

Harvard College Observatory, 

Cambridge 348 
Vass., U.S.A 
1955 March 15 
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NOTE ON A SET OF MODELS FOR UNMIXED STARS 
A. E. Roy 
(Communicated by the Director, University Observatory, Glasgow) 


(Received 1955 March 24) 


Summary 


An erroneous equation in a previous paper is corrected and leads to the 
conclusion that the radii of the unmixed stars are smaller and the effective 
temperatures higher than was formerly indicated. 


1. Introduction.—In a previous paper (1) the present author derived a set of 
point-convective models of a star with a range of discontinuities of molecular 
weight between core and envelope. It was found that the radius of such a star 
could increase by a factor of ten owing to the exhaustion of hydrogen during 
evolution. 

It has been pointed out to the author by Dr Donald E. Osterbrock of Mt Wilson 
and Palomar Observatories that the use of equation (3.10) is erroneous in the 
present case where a discontinuity in molecular weight arises at the interface 
between core and envelope. Equation (3.10) was adopted from equation (14) of 
a paper by C. M. and H. Bondi (2). Equation (14), correct in that paper, requires 
modification when used with discontinuous model stars. 

The error, in fact, arises earlier. In equation (3.03) of the previous paper (1) 
the contribution from the discontinuous jump in density across the interface 
between core and envelope was forgotten. 

The equation should be 


N 
In =| dut+in (1.01) 


This changes expression (3.10) to the following : 


Me 


=InQ—u+ | du du +in (1.02) 
‘This correction leaves the set of tables calculated from the homology equations 
unchanged but changes the final reductions by which R/R,, L/L, and 7,/T., 
were found. If R,, and 7,, refer to the corrected values, then they 
are given by the following expressions : 


R’ R, R R, (@)" 


(2) 


< 
ae 
i; 
L'/L,= LiL, (1.03) 
T/T ty A 
Na 
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where a, B, ay, as, etc. are defined as before and R/R,, L/L, are given in Fig. 1 
of the previous paper(1). ‘The changes in R/R,, L/L,, T,/T.,, for a discontinuity 
in molecular weight of ,/u, = 2-64 are given in Table I below. 

Tasie 


Mel He = 2°64 


Uncorrected Corrected 
(paper 1) from (1.03) 


10 
1°35 
1°24 


An unmixed model for y,/4,=2°091 was interpolated from the corrected 
values and compared with a similar model kindly sent to the present author by 
Dr D. E. Osterbrock. ‘The results agreed, showing that the numerical part of 
paper (1) is accurate, viz. the set of models for unmixed stars calculated from the 
homology equations. 

2. Conclusion.—The fact that the radius of such a model does not increase 
by a factor of ten and the luminosity by a factor of 1-35 during evolution changes 
the H-R diagram of such stars obtained in the previous paper (1). ‘The changes 
in log 7, T,, and m—m., are shown below for y,/, = 2°64 where suffix a refers to 
the Cowling model. 

Il 


2°04 
Uncorrected 
Corrected 
(paper 1) 


log T,/T,4 0°45 


Thus as p, », grows from 1 to 2°64, the star does not move into the red giant 
region of the H-R diagram but remains relatively close to the line AB. ‘The 
previous conclusion that the two populations of mixed and unmixed stars do not 
mingle in the H-R diagram remains unchanged. 


My thanks are due to Dr Osterbrock for drawing my attention to this matter, 


University Observatory, 
Glasgow W.2: 
1955 March 1%. 
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NOTE ON THE OBSERVED DIFFERENCES IN SPOTTEDNESS 
OF THE SUN’S NORTHERN AND SOUTHERN HEMISPHERES 


H. W. Newton and A. S. Milsom 
(Communicated by the Astronomer Royal) 
(Received 1955 April 14) 


Summary 


The variation with time of the relative spottedness of the Sun’s two 
hemispheres, north and south, is examined by means of a simple index 


over a period embracing eleven 11-year sunspot cycles. ‘There are changes 


from cycle to cycle which, although not random, appear to have no definite 
period. ‘There is the suggestion of a dependence of the pattern of the excess 
spottedness, north or south, upon the degree of activity of any particular 


cevcle The results indicate no dependence upon a 22-year cycle 


Introduction. —Sunspot observers have long been aware of persisting periods, 
extending for months and occasionally to years, when there has been a marked 
inequality in the spottedness of the two hemispheres of the Sun, north and south.* 
The presence of any systematic or periodic behaviour of this phenomenon, if 
established, would be of interest as a new observational feature of sunspot 
behaviour and of immediate interest to current theories of sunspots. ‘he time 


covered by the data hitherto considered has been too short to suggest other than 
a random behaviour in the difference already noted between the spottedness 
of the two solar he mispheres. 

There are, however, now available from the Royal Greenwich Observatory 
the northern and southern annual mean daily areas of sunspots (corrected for the 
effect of surface foreshortening and expressed in millionths of the Sun’s 
hemisphere) from 1874 to 1954, thus comprising seven complete sunspot cycles.t 
It seems therefore worth while pointing out in this note certain features now 
apparent in this N— 5 spot occurrence over the period of 80 years covered by the 
Greenwich observations, and seeing to what extent the features can be recognized 
in the older records extending back to 1833. ‘The Greenwich observations will 
be discussed first. 


* Newcomb, 5., Ap. J., 13, 13, 1901, who drew attention to a general excess of southern spots 


from 1866 to 1898 Royal Observatory, Greenwich, VW.N., 63, 452 


, 1903. It is noted that for the 
two cycles ending 1902, the northern spots predominated during the rising phase and the southern 
spots towards the end of the eyck Brunner-Hagger, W. and Liepert, A., Astronomische Mit 
teilungen, No. 140, p. $56, 1941. The authors have studied the distribution of log N 5S within each 


cycle irrespective of time (1853-1999). Maunder, FE. W., W_N., 82, 534, 1922. Kiepenheuer, K. O 
The Sun, ed, G. P. Kuiper, p 


136, Chicago, 195% 


t The annual volumes of the Greenwich Photoheliographic Results have not yet been issued by 
H.M.5.0O. after 1945. A special volume that will contain the collected data 1874 
of publication, entitled Greenwich 


1954 is in course 
Sunspot and Geomagnetic-Storm Data, 1874-1954". The 
requisite tables 1946 to 1952 inclusive have also been published in Monthly Notices 
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(a) Yearly values of (N—S)/(N + S) (in terms of area) for each of the seven cycles 
1879-1954. Points where N + S~< 100 units of area are joined by broken lines. (bh) Means of 
adjacent cycles shmon in Fig. 1 (a). (ce) Means of adjacent cycles shown in Fig. 1 (a), alternate 
(d) Mean of the seven cycles shown in Fig. 1 (a). 
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Data.— Using the Greenwich values of the annual mean daily areas of sunspots, 
the difference E, of the ratios N/(N +S) and S/(N +S) has been taken for each 
year and plotted in Fig. 1 (a) for the seven sunspot cycles. In each case, the 
time origin has been adopted as the year of sunspot minimum, as given by Ziirich 
sunspot “numbers” in agreement with the Greenwich areas and the frequency 
of spotless days, Any epoch of minimum is admittedly partly dependent on the 
old cycle as well as the new but is better for our purpose than the epoch of maximum, 
which is sometimes rather indefinite as in cycle No. 3. At minimum, there 
is the well-known overlap in latitude of the old-cycle spots in low latitudes and 
the new-cycle spots in latitudes usually above 20°. ‘The spots in the period of 
the overlap have therefore been allocated to their appropriate cycles. Although 
these points in Fig. 1 (a) at or near minimum appear significant in giving trends, 
the points are numerically of low weight because of low sunspot frequency. 

Greenwich areas of faculae.——\t is of interest to compare the yearly values of 
EE, with those derived from faculae. Although faculae in the general sunspot 
zones are closely related to the spots, their dark streaks or masses (on the negative) 
present to the measurer a very different feature from the sunspots themselves. 
Moreover, faculae are longer-lived than their associated sunspots, but are usually 
visible only within three or four days of the Sun’slimb. Areas of faculae measured 
at Greenwich, divided according to the northern and southern hemispheres, 
are available* for deriving a similar ratio to that for the spots for the years 1874 
to 1916. ‘The agreement of the values given respectively by sunspots and by 
faculae is remarkably close (Fig. 2), and indicates that the departures of the points 
from a smooth curve in Fig. 1 (a) are not due to observational errors but are of 
solar origin. 

Counts of sunspots 1832 to 1874.—-Prior to the Greenwich observations 
commencing in 1874 April, there are no published areas of sunspots divided 


according to hemisphere, but there are the following dependable series of 
drawings: 


Schwabe’s records 1826 to 1867 preserved at the Royal Astronomical Society. 
Carrington’s published observations 1853 November to 1861 March. 
Sporer’s published observations 1861~1879.] 

The Kew series of photographs§ extend from 1862 to 1873 but the whole of 


the observations do not appear to have been published in the form necessary for 
this investigation. 


Before using any of these older records it was first ascertained from the extensive 
Greenwich records whether the yearly values of FE, derived from sunspot areas 
were satisfactorily represented by similar values Ey, obtained by counting each 
spot group once, whether it was seen for a few days only or for the whole disk 
passage. A recurrent group was counted as another group (i.c. as in the Greenwich 
results, in which the return group has a separate number). All groups seen only 
for one day were excluded. 

As indicated by the high degree of correlation exhibited in Fig. 3 between 
EF, and Ey over the period 1874 to 1954, it appears well justified to extend the 

* i. W. Maunder, Memoirs of British Astronomical Association, Vol. XXI11, p. §7 


1 R. ©, Carrington, Observations of thi Spots on the Sun from November 9, 1853, to March 24, 
made at Rediull, 1864 

1G. Spérer, Publicationen der Astronomische Gesellschaft, 13 (Leipzig, 1874); 13, Fortsetzung 
(Leipzig 1876); Pub. Astrophys. Obs., Potsdam, No, 1, 1, Part 1; No. 5, 2, Part I. 


§ Phil. Trans., 1§9, 1, 1869; 160, 389, 1870, 
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data backwards from 1874 by means of spot counts, north and south. Schwabe’s 
published observations in Astronomische Nachrichten give counts for the whole 
disk only, so that counts of spot groups, according to hemisphere, have had to 
be made on the original drawings. From about 1832, these drawings were made 


(N—S)/(N+S) 
(faculae) 


(N—S)/(N +S) 
(spot areas) 


Comparison of (N S)/(N + 8) (faculae) (full curve) with 
S)/ON +S) (spot areas) (broken curve) 1875-1916 


3 
° 
y 4 
a 
Z 


(N—S)/(N-+-S) (areas) 
Fic, 3.—Scatter diagram to show corre- 
lation of +S) (spot counts) and 
(N—S)/ON + 8S) (spot areas). 


(N—S)/(N+S) (in terms of spot counts 


(Ean 

Fic. 4.-— Yearly values of + S) (in 
terms of spot counts) for cach of the four cycles 
1833-1878. Points where so are joined b 

broken lines 
on a scale that was substantially constant, the solar disk being represented by a 
circle 2:1 inches in diameter. Copies of the Stonyhurst Disks reduced to the 
same scale were used to read off the latitudes of the spots. Permission to examine 
the records was kindly given by the Royal Astronomical Society, ‘The counts 
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were made wholly by one of us (A. 5. M.) but sample checks were made. As 
in the Greenwich counts, a group, whether of short or long duration, was counted 
once, group numbers having been assigned in the records both by Schwabe and 
Sporer. Counts were similarly made from the published records of Carrington, 
Sporer, and the Kew publications. After finding that there was satisfactory 
agreement of Ey over the few years common to these last three series of data, 
those adopted were 


1833 to 1867 Schwabe, 
1866 to 1878 Sporer (these overlap with the Greenwich records). 


‘The Schwabe data from 1826 to 1832 inclusive were not used, owing to uncertainty 
in the position of the north point of the image on the drawings. 

The yearly values of Ey thus derived from 1833 to 1878 are plotted in Fig. 4, 
the time origin of each of the four cycles being adopted, as before, as the year of 
sunspot minimum. ‘The separation of old- and new-cycle spots has been made, 
according to their respective latitudes, as was done for the Greenwich data. 

Discussion.-Although it is believed that a reasonably homogeneous series 
of data is here presented from 1833 to 1954, we prefer to discuss first the last 
seven 11-year cycles covered by the Greenwich observations for which the daily 
record is virtually complete. 

Examining Fig. 1 (a), we see that using an index derived from areas (E,): 

1. The values are subject to rather large fluctuations. It has been shown, 
however, that these are real since they are reproduced by the faculae and numbers 
of spot groups. 


2. If for each cycle a straight line be drawn that best fits the points 
(unweighted), a progressive change of slope from negative to positive (the sense 
being defined by taking north as positive) is apparent from Cycle 1 through to 
Cycle 7. 

The reality of this change of slope is confirmed by Fig. 1 (6) and (c). The 
former displays the means of adjacent pairs of cycles (1 and 2: 2 and 3...6and 7), 
and it will be seen that the progressive change of slope is still preserved. In 
Vig. 1 (c), however, the procedure has been repeated with alternate cycles inverted 
(i.e, rand —2: zand3: ... ~6and 7), and a progression is no longer apparent. 
If these data did not represent a significant variation in the behaviour of each of 
the two hemispheres, then this procedure of inverting alternate cycles could not 
be expected to alter the character of the results obtained in Fig. 1 (4). Fig. 1 (d), 
which gives the mean curve for all seven cycles, is informative, especially in 
showing that over the period of nearly eighty years the amount of spottedness 
is very nearly the same in the two hemispheres throughout the averaged cycle. 

3. The slowly changing pattern of E, from cycle to cycle implies that the 
relative behaviour of the two hemispheres is maintained in successive cycles, 
with no indication of an alternation as might have been suspected from the 
22-year cycle of magnetic polarities. 

With these fairly definite features shown from Cycles 1 to 7 covered by the 
Greenwich observations of spot areas, it is obviously desirable to extend the 
observations backwards in time by means of the data provided by the earlier 
sunspot observers. ‘The results of the differences in the ratios of spot counts 
N/(N +S) minus 8/(N + 5), Ey, are plotted in Fig. 4, which is directly comparable 
with Fig. 1 (a). 
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It seems that the progressive change of slope of E, is not continued back in 
time as might have been expected if the progression in Fig. 1 (a) had been part 
of a periodic change. ‘There is, however, one feature of interest concerning the 
whole series extending over 120 years, exhibited in Fig. 5. Excluding points 
for which N +5 was less than 50, numerical values of the slope of each cycle 
using FE, throughout have been derived by calculating the slopes from each point 
to the centre of the curve and taking the mean, each point being weighted in 
proportion to its distance from the centre (this method, for the present purpose, 
is sufficiently accurate). Each of the eleven values thus obtained is plotted in 
Fig. 5 using full circles. In the same diagram is plotted for comparison the 
annual mean daily sunspot number (R) of the maximum year for the appropriate 
cycle. ‘There appears to be a marked correlation between the two curves. “Thus 
a very active cycle as defined by a high maximum (either of R or mean daily spot 
area) tends to be associated with a positive slope, i.e. the southern hemisphere is 
the more spotted at the beginning of the cycle and the northern at the end, 


Sunspot number (Zurich) 


Change per year (N—S)/(N+S) 


Cycle number 
Pic. s.—Full curve : slopes of UN — S)/CN + S) (counts) for each of the elewen cycles 1644-1954. 
Hroken curve annual mean sunspot number for the maximum year of each cycle, 


Conclusions.—Vhe observational results discussed above suggest that the 
relative amount of spottedness of the northern and southern hemispheres of the 
Sun follows some pattern which cannot be accounted for in terms of the random 
fluctuations that would be expected in the data. ‘This pattern does not appear, 
however, to be periodic. ‘The spot areas of the homogeneous series of Greenwich 
photographic observations extending now over seven complete 11-year cycles 
suggested at first that a periodic change might be present. But on extending the 
observations backwards to 1833 by means of the records of the visual observers 


of sunspots, Schwabe, Carrington, and Sporer, the promising pattern given by 
the Greenwich results 1874-1954 does not appear to be continued, though the 
limited period (four cycles, 1833-1878) leaves this a little inconclusive. 
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A definite result that emerges is that over the whole period of observation 
(120 years) there is no suggestion of any effect related to the 22-year cycle of the 
magnetic polarities of sunspots. ‘The apparent correlation between the computed 
slopes given by Ey values for each cycle and the height of the peak of the cycle is 
noteworthy. 

A detailed statistical analysis of the data would be desirable to establish the 
degree of significance of each of the effects suggested here, and perhaps also to 
enquire into other possible effects given by these data. It is hoped that this note 
will stimulate such work, 


Our thanks are due to the Astronomer Royal for permission to publish this 
note. We are also indebted to Mr T. Gold, Chief Assistant, for his interest and 
helpful discussion. 


Royal Greenwich Observatory, 
Herstmonceux Castle, Sussex : 
1955 March 
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INTERFEROMETRIC MEASUREMENTS OF WAVE-LENGTHS 


II. MEASUREMENTS AT 6500 A AND A GENERAL DISCUSSION OF THE 
SOLAR RED SHIF1 


(Communicated by the Director of the University Observatory, Oxford) 


(Received 1955 April 18) 


Summary 


It has been suggested that the largest wave-length displacement of solar 
lines to the red will be observed for resonance lines. ‘To test this prediction 
the wave-lengths of a number of lines, including the calcium inter-system 
resonance line 4'S-—4"P, have been measured in the 6s00 A region by the 
method of circular channels. Contrary to the theoretical prediction the 
resonance line has the smallest red shift which has yet been observed 

This resonance line is very faint and it appears to exemplify the pheno 
menon (frequently referred to by Burns) of a variation in line displacement 
with line strength In a discussion of the red shifts at 6500 A, together 
with the earlier results at 5080 A and 6020 A, it is clear that line strength 
is the dominating factor affecting the shifts. In order to make a quantitative 
examination of this effect over all wave-length regions, line strengths are 


expressed in equivalent “ ergs’"’, a unit proportional to the absolute energy 
absorbed by the line if A (rulliangstroms) is the solar red shift at wave 
length A it is shown that there is a linear reiation between A/A* and the 
logarithm of the equivalent “ ergs’’, a relationship which is applicable at 
all wave-lengths. 158 solar red shifts measured in the Allegheny Observatory 
Bureau of Standards work are analysed in exactly the same way as the 20 
Oxford results. ‘These confirm the empirical relationship over a wave-length 
range extended by tooo A, and a greatly increased range of line strengths 


Introduction, —The red shift observed for the solar lines is attributed to a 
combination of the predicted relativity shift and possible low-pressure collisional 
shifts due to neutral hydrogen. ‘The relativity shift is precisely determined as a 
wave-length change of 2-12 parts in 10°, but the collisional shifts, which are of the 
same general order of magnitude, are so uncertain in exact magnitude and sign as 


to prevent any real confirmation of this theory of the solar red shifts. ‘The least 


uncertain collisional shift is that for a resonance line. For such a line we should 
expect a collisional shift to the red (1g), and therefore a solar red shift in this case 
in excess of the relativity shift. 

The work now described extends the earlier measurements (1) to the red in 
a region which includes a resonance line. Section I gives an account of the 
determination of absolute wave-lengths for solar and vacuum-arc lines by the 
method of circular channels, and the results obtained. In Section II the red 
shifts, now determined over a considerable range in the visible spectrum, are 
discussed. It is found that the outstanding observational characteristic of the 
solar red shift is its dependence upon line strength, rather than on any other 
factor. 
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SECTION I 
Wave-length determinations 


1.1. Region.—-\n selecting a region the first requirement is that it should 
contain a resonance line measurable both in the solar spectrum and in the 
laboratory source. A study of Miss Moore’s Multiplet ‘Table (15) shows that 
while there are 60-70 resonance lines of neutral elements in the solar spectrum 
between 4000 and 7000 A, the number of lines possible for the present problem 
is nearerto 10. ‘The resonance line finally selected was the calcium inter-system 
combination line 4'S-4*P° at 657283 A. The line is unfortunately rather 
faint in the Sun (Rowland number 1), but such an inter-system line should be 
intrinsically narrow, and we need not fear broadening due to hyperfine structure, 
since “Ca, for which the nuclear magnetic moment is zero, is 97 per cent of 
naturally occurring calcium. ‘The line is fairly easily obtained in the vacuum arc, 
the only disadvantage being again that it is somewhat faint. ‘The region 6580 A 
has the very great advantage that it is close to the primary standard at 6438 A, 
so that uncertainties in the measured wave-lengths on account of chromatic 
variation of phase change, camera focal length, and shift of pattern centre are all 
negligible. 

In the 6580 A region the lines given in the lower part of ‘Table I have been 
measured in the solar spectrum. ‘The faint water-vapour line and the faint 
iron lines have been included to give an idea of the accuracy obtainable, for lines 
of intensity comparable with that of the calcium line. ‘The seven lines in the 
upper part of ‘Table | have been measured to extend the solar observations in the 
red, and to give some measures of the red shift of calcium lines generally, in case 
of any peculiarities in the behaviour of the element. ‘The two calcium lines in 
this region are the measurable lines closest to the resonance line. Columns 1, 
2 and 3 in ‘Table I are from the Rowland Revision (23). ‘The next four columns 
are from the 1928 1.A.U. table of standard solar wave-lengths (24); they show 
the recommended wave-length, the measures of St John, of Babcock and of Burns, 
Meggers and Kiess in the Allegheny Observatory Bureau of Standards work. In 
column 7 the A.O.-B.S. measures have been augmented by the values in brackets 
from the A.O.-B.S. papers (7, 9, 10) where the 1.A.U. table has no entries. 

| 


Measured lines in the solar spectrum 


Rev. Rowland ‘Table LA.U. St John Babcock A.O.-B.S. 

El. Int. Wave-length 1928 
nA 

6490°798 pan eos (-790) 
Ca 6 6493°796 “788 “788 
Ke 6495 °001 4°994 994 993 ‘995 
Fe 2 64960480 (469) 
Ba! 4 6496°916 7 (-g11) 
re 6498-954 945 ‘045 ‘945 
Ca 4 6499°063 054 655 654 
Fe 5 6569°232 ‘224 ‘224 "223 
Ca i 6572 804 see (°736) 
Fe 6574°243 (234) 
Fe 2 6575°045 (022) 
Fe 6 6592°934 926 926 926 
re 4 6593 
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1.2. Observational procedure.—The method of circular channels has been 
used for the determination of absolute wave-lengths in the spectrum of the Sun 
and for the same lines, as far as possible, in the vacuum-are spectrum. The 
experimental procedure follows very closely that already described in Paper I (1). 

The dispersion of the Oxford prism spectroscope has dropped to o-5 mm A ! 
at 6550A, and moreover we should expect the slit pattern to be broadened 
relative to the dispersion compared with the instrumental performance at shorter 
wave-lengths (6). Both factors increase the range of wave-lengths presented to 
the etalon in the slit image. It is not practicable to work with a spectroscope slit 
width much less than 0-1 mm on account of exposure times and the narrowing of 
the heterochromatic channels in which the absorption lines must be measured. 
This means that a 2:5 mm spacer is the longest it is permissible to use in the inter- 
ferometer without serious scatter from order to order (20); even shorter spacers 
are desirable from this point of view. 

Work was also attempted with a 1-5 mm spacer, and a full discussion of these 
results, which proved to be unreliable, is given in Paper III] (to follow). For the 
reasons there explained, errors in the absolute wave-lengths measured with a 
2°5 mm spacer and longer are likely to be small, and errors in the red shifts only 
accidental. ‘The measurements at 6580 A, the region of the resonance line, 
were repeated after the work described in Paper IIl was done. For the lines 
656g°2, 6592-9 and 6593:-9A which are measurable with fair certainty, the 
differences in the red shifts obtained by the original and by the improved procedure 
were less than o-oo: A. ‘The slight changes in the experimental arrangements 
from those described in Paper I (1) are discussed in Paper III. 

The results now given were all obtained with a 2:5 mm spacer, with inter- 
ferometer plates silvered to give just over go per cent reflection in the red. The 
plate silverings for the 6580 A region were slightly heavier than for the 6490 A 
region. 

Solar plates.—Vhe solar exposures were all obtained at the centre of the 
180 mm image in focus at the spectroscope slit. With slit widths from o-o8 mm 
to the exposure times at 6580 A were 8-12 minutes on Ilford Rapid 
Process Panchromatic plates. The sensitivity of these plates is falling off quite 
rapidly at this wave-length, but there is very little difference in sensitivity between 
the I.R.P.P. plates and Astra 7 or Astra 8, which were also tried. In the 6490 A 
region with a slit width of o°og mm the exposure times were 14-2 minutes on the 
I.R.P.P, plates. Kodak Press Contrast developer was used throughout. ‘The 
plates for the 6490A region were obtained in 1952 August, and those used for 
the 6580 A region in 1954 April and August. 

Vacuum-arc plates..-Whenever possible the lines given in Table | were meas- 
ured also inthe vacuum arc. ‘The vacuum-arc wave-lengths measuredat Allegheny 
Observatory for the lines of Table I are given in Table Il, together with King’s 
intensity estimates. ‘The iron lines have been measured by Burns and Walters 
(11), the calcium measures are those of Wagman (21). 

With our experimental arrangements an iron line of arc intensity 2 is on the 
limit of measurability. Some vacuum-arc measures have been made for the 
line 6575°045A in Table | (this is a line of King’s estimated intensity 2) but 
only a computed value is given by Burns and Walters. With the arc conditions 
we aim at, 1c. a pole separation of not less than 10mm, currents not exceeding 
8 amp. and pressures of a few cm of mercury, exposures of 4-5 hours are barely 
sufficient for such lines. For the faint arc lines in the 6580 A region the exposure 
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times were reduced for some plates by using the spectroscope undiaphragmed 
(17). In the 6490A region only the strong iron line at 6494°981 A has been 
measured, the exposure time required being only 5 min; the fainter line at 
6496464 A was not visible on these plates. 

The definitive calcium arc plates were obtained using graphite or carbon 
pole-pieces, the positive pole cored and filled with calcium chloride. The arc 
so produced would not always draw out to the full 1omm, and for the calcium 
plates the pole separation was 6-10mm. For the strong calcium lines at 6490 A 
the exposure times were less than 2 min, but at 6580 A, for the resonance line, 
with slightly heavier silvering on the interferometer plates, about one hour was 
required with the diaphragmed spectroscope and approximately 20 min with 
the diaphragm removed. Some of the earlier work was done with a rod of 
calcium metal in a copper holder as the positive pole-piece and a copper negative 
pole-piece. ‘This was a much more convenient source to use than the calcium 
salt, but calcium metal in a form suitable for turning is not easy to procure, and 
at 6490 A there is a possibly overlapping copper line. 


Vacuum-arc lines 


"TABLE 


(Allegheny Observatory measured wave-lengths) 


Iron lines Calcium lines 
Wave-length Int. Wave-length Int. 
in A mA 
6494'951 25 6493°750 Bo 
6496°4604 5 6499°649 


65609215 
65§92°91 10 
6593 °575 


A rotary vacuum oil pump was used to maintain the pressure in the are at 5 cm 
of mercury or less, as read by a closed-arm U-type mercury gauge mounted on 
the arc. ‘The readings of such a gauge are subject to uncertainty owing to the 
possible presence of condensable vapours in the closed limb. ‘The gauge was 
therefore tested by means of an absolute aneroid type gauge; this showed that 
the pressures actually used had been about 8 cm of mercury when the arc was cold, 
rising to 14cm when it became hot. 

Ilford HP3 plates and Kodak Press Contrast developer were used for all the 
are spectra. ‘The arc plates at 6490 A were obtained in 1952 August and those 
at 6580A in 1954 July and August. 

1.3. Plate measurement and reduction..The plates were measured on the 
Cambridge machine fitted with a rotary stage as already described (1). A modi- 
heation of the rotating table has been made which considerably eases the subsequent 
reduction of the measures. It enables the cadmium pattern centre O (see Fig. 1) 
to be brought within a few thousandths of a mm over O’, the centre of rotation 
of the table. ‘The plate is clamped to an upper stage which can be turned by 
means of a tangent screw over a lower stage to which it is pinned. ‘The lower 
stage is itself pinned to the rotating table and can be turned over it by a tangent 
screw with a motion perpendicular to that of the upper stage. ‘These motions 


' 
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are used to secure that the four arcs of any one cadmium ring have very closely 
the same setting when they are turned into the line of travel of the microscope. 
In this way the distance OO'(p) can easily be reduced to under oo1mm. This 
means that the difference between the micrometer settings of the two positions 
in which any one arc cuts the line of travel of the microscope is, very nearly, the 
true diameter of the interference ring, assuming no shift of patterncentre. A small 
correction depending on p and the distance of the spectrum line from O may be 
necessary, but the solution of a triangle for each measured arc is no longer required, 


and angular readings for the arcs are not necessary, except for identification 
purposes. 


Fic. 1.--Plate reduction 


The reduction now proceeds as follows. Suppose we have a system of 
rectangular coordinates with origin at the cadmium pattern centre O, and x-axis 
along the arms of the cadmium cross which lie parallel to the spectroscopic 
dispersion. Any spectrum line then lies approximately parallel to the y-axis, 
and we suppose that two arcs of the line interference pattern lie at P and P’, 
(x, y) and (x, —y). From the measures in the cadmium pattern we can determine 
©’, the centre from which the measures have been made. [Let O' have coordinates 
(x’, y’) and let the measured distances O'P, O'P’ be s, and s,. Assuming that O 
is also the centre of the line pattern, as we may for this region so close to the 
cadmium red line, then the diameter D of the ring of which P and P’ are parts 
is given by 
OP’? = 25, 255 + 4(25, — 254)" — 4p? + 2x. 4x". 
Because 0’ is very close to O the second and third terms are negligible, and we get 
= 28, 289+ 2x. 4x" 

Our measures give 2s, and 2s, directly. In the very small correction term, 
4x’ comes immediately from the measures in the cadmium pattern arms lying along 
the x-axis (4% is 2Ar in Paper I, Section 2.1 (1)), and the x-coordinate of the 
spectrum line need only be known very approximately. We note further that the 
correction term is constant for all the rings of any one line, for a given O’. 
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The rest of the reduction follows exactly as described in Paper | with the 
considerable simplification that we have no corrections for the chromatic variation 
of phase change or camera focal length. 

1.4. Observed wave-lengths.—The absolute wave-lengths obtained for the 
absorption lines are given in Table III. The solar lines have been corrected for 
diurnal rotation and the Earth’s orbital motion. ‘The wave-lengths are mean 
values for the 18 mm spectroscope slit height and refer to a line of points extending 
over ot of the solar diameter at the centre of the disk. ‘The measures have all 
been made with a 2°5mm spacer (approximate order of interference 8100). 
At 6490 A the lines fall close together and in the ‘normal’ position they were all 
measured on the violet side of the pattern centre (the same side as the centre of 
the heterochromatic fringes). Plates A412, A413 were obtained with the lines 
shifted to the red side of the pattern centre. At 6580 A the lines in the ‘ normal’ 
position fall on both sides of the pattern centre and the measured plates were all 
obtained with only slight variations about this position. 

Wave-lengths are given to o-oo1 A in ‘lable II1, though all the reductions have 
been carried through to o-ooo1 A; the means in column 3 and the standard 
deviations of the mean in column 4 are based on the four-figure values for the 
individual plates. ‘The faint lines at 6580A were very difficult to measure ; 
we note that the plates obtained in 1954 April give wave-lengths in the mean 
0002 A longer than those from the 1954 August plates (or o-oo1 A longer for the 
better-determined lines), On the basis of the investigations described in Paper 
111, the April values are somewhat less reliable than the August values. ‘The mean 
wave-length in column 3 is the straight mean for the nine plates. 

‘The wave-lengths for the are lines are given in Table [V. ‘These again were 
all obtained with the 2°5 mm spacer. In the 6490 A region the first two columns 
of the individual plate values were measured in the ‘normal’ position, and the 
last three in the ‘shifted’ position. ‘The lines at 6580A were all measured in 
the normal position, ‘The plate values in the first two columns at 6580 A were 
obtained with the diaphragmed spectroscope and those in the last four columns 
with the full aperture. It will be seen that there is no change in the measured 
wave-lengths on this account. 

Comparing the mean values in ‘Tables II] and IV with those in ‘Tables I and 
Il, we see that the solar wave-lengths now obtained are on the average o-oo1 A 
shorter than Babcock’s values, and than the A.O.-B.S. values. The seven arc 
wave-lengths average 0-002 A shorter than the Allegheny Observatory values, 
though the differences are somewhat erratic. We should have expected the rather 
higher pressures now used in the vacuum arc to have resulted in longer wave- 
lengths rather than in shorter ones, 


Section 

A general discussion of the solar red shift in different spectral regions 

2.1. An analysis of the results so far obtained.-We have now measured absolute 
wave-lengths for solar lines and vacuum-are lines over the wave-length range 
5080 A-6580 A, and in this section we will consider the red shifts shown by the 
solar lines. ‘The results from Paper | and those nowgiven in Section! are collected 
together in Table V. Columns 1, 2 and 3 give the element, Rowland number and 
identifying wave-length for the line. ‘The next two columns are from Miss 
Moore's 1945 Multiplet Table (16) and give the multiplet number of the line 
and the excitation potential in volts of the lower state involved in the transition. 
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Veasured wave-lengths in the 


Plate No [Fe A422a Ag2z2zb Agz4a 
(Ca 395a Azosb Agz2 A4gz3za 
Mean obs Standard 
kl wave-length deviation | 
Int. | 
mA of mean 
Ca 6 6493°777 0°0004 | °777 776 777 7799777 
6494°9%0 0°0003 gko go go 
Ca 4 6499°045 | 645 644 645 645 645 = 
| Aa Ay 
Plate No. 4 Fe | A750 A753 A760 A767 A772 A775 
(Ca | A752 4754 A765 A776 


vacuum arc (individual plate values) 


Mean obs 


Standard 
how 
wave-length deviation 
Int 
in A of mean 
ke 5 6569°216 + 00004 218 ‘216 216 "215 216 
Ca 6572°778 00002 778 778 778 ‘779 ‘777 778 
Fe 2 6575°O17 0°0008 O15 O17 ‘018 
Fe 6 65§92°913 00004 ‘O12 git O14 ‘O13 
Fe 4 6593869 0°0003 $70 ‘$68 ‘868 870 


Column 6 gives the fractional angstroms of the measured arc wave-length. 
As we saw in Section 1.2, the pressure estimated for the vacuum arc is of the 
order } of an atmosphere, and this is applicable to the earlier as well as to the 
present measures. ‘The line shift at such a pressure does not often exceed 
ooo! A, but it is nevertheless desirable to correct the results for the small pressure 
shift as far as is possible. In order to estimate the correction we have used 
Babcock’s systematic study of the shifts in iron lines between the vacuum-arc 
and the arc with pressures up to one atmosphere (§). Babcock found that the 
shifts were always to the red as the pressures increased, and he derived an empirical 
relation between the depression d(cm~') per atmosphere of a spectral term and 
the excitation potential V (volts) of the term. For triplet terms he deduced the 
relation d= 1-93V, and for quintet and septet terms d= 1-15 V*—1-27V. 
Ihe shift for any line is given by the difference between the depressions for the 
two terms involved in the transition. Such an observational relation is really 
necessary for each clement individually, but for the small pressure correction 
now required we have used Babcock’s results for the iron lines for the calcrum 
and manganese lines as well. Some pressure shifts have been given by Monk 
(14) for a number of manganese lines: they are of the same order of magnitude as 
those now used. For the calcium lines no such measures appear to be available. 
Corrections ranging from zero to 00018 A have been applied to the measured 
arc wave-lengths to give the corrected vacuum-arc wave-lengths in column 7. 

Column 8 of Table V gives the measured solar wave-length, and the difference 
between column 8 and column 7 gives the observed solar red shift in column g. 
Chis quantity, expressed in milliangstroms, we will denote by A in the discussion 
following. ‘The line equivalent width in column 1o is that given by Allen (2); it 


is his value in mA derived from r, the central intensity of the line corrected for 
blending. 
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The first effect we look for in the results is a change in the red shift with 
wave-length such as is predicted on the relativity theory. An inspection of the 
table shows, however, that there is no correlation of the mean shift in each region 
with the wave-length of the region. ‘The mean shifts for the wave-length groups 
are 10°1 mA at 5080 A, 7°7 mA at 6020 A, 11-6 mA at 6490 A and 9-4 mA at 6580 A. 
Nor is there any obvious dependence of the red shift on excitation potential. 
We note that the calcxum resonance line 4'S—4°P° at 6572°8A, which from 
Spitzer's prediction should have the largest red shift as a result of the relativity 
shift and a red collisional shift, has in fact the smallest value of A measured. ‘The 
dominant feature of the results is the change in red shift with line strength, and 
this shows to some extent within each wave-length group. ‘The number of lines 
in any one group is too few, and the observational uncertainty of the red shift of 
a single line too great, to derive a relation between line strength and red shift 
at any one wave-length; we must therefore discuss the results obtained in the 
different spectral regions simultaneously. 

The equivalent width of a line in milliangstroms (or in cm ') is a variable 
unit in which to express line strengths, because it refers to the continuous back- 
ground intensity which changes with the wave-length region. If we express the 
strength of a line in absolute units of energy absorbed from the continuous 
spectrum, then we have a unit which has at least the same physical significance 
at all wave-lengths. We have therefore expressed line strengths in a unit which 
we will call equivalent “ ergs’’, by allowing for the variation in the continuous 
background intensity. Column 11 in ‘Table V gives the ratio of the continuous 
spectrum intensity /,(0) at the line wave-length A to that at 5000 A, using 
Minnaert’s table (13) of the observed normal monochromatic intensity /,(0) per 
dA=i1em. ‘The equivalent “ergs” in column 12 of Table V are obtained by 
multiplying the equivalent width in column to by the factor in column 11. ‘The 
absolute magnitude of the line strength unit is for our purpose quite arbitrary ; 
we have aimed at using a unit of the same order of magnitude as the familiar 
equivalent width in milliangstroms, and the two units are in fact equal at 5000 A. 

big. 2 shows a plot of A, the red shift in mA, against the logarithm of the 
equivalent “ergs’’. We see that there is a reasonably well-defined increase of \ 
with line strength, though the scatter of the points is quite considerable. Assum- 
ing a linear relation between the two quantities the computed correlation 


coefhcient is +0509. ‘The corresponding significance level is very nearly 
0°02 (12), that is, there is a probability of m that a coefficient of this size would 
result purely from chance. ‘The two regression lines for the distribution are 


shown in the figure. 

A wave-length dependence of \ would show in Fig. 2 by the systematic 
departures of the As in the different spectral regions from the mean curve, and 
there is some indication of such an effect. If A is proportional to wave-length A 
(as on the relativity theory or for a Doppler shift) then we should expect to 
reduce the scatter of the points by plotting A/A against line strength instead of A. 
Fig. 3 shows A A (A also in milliangstroms) plotted against log equivalent “ ergs’’. 
There is obviously some reduction in the scatter of the observed points. ‘The 
computed regression lines shown in Fig. 3 are much closer together and the 
correlation coefficient is now + 0°683, corresponding to a significance level better 
than o-oo1. In Fig. 3 the relativity shift corresponds to \ A= 21-2 x 10 
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We conclude then that A is dependent upon wave-length as well as upon line 
strength, though in an empirical approach we can only guess at the form of the 
two dependences. If we next assume that A is proportional to A® and plot A/A* 
we are expressing the shift in units proportional to frequency; Fig. 4 shows 
4/A*(in cm *) plotted against log equivalent “ ergs”, with the computed regression 
lines. ‘This again reduces the scatter in the observed points from that in 
Fig. 2, the correlation coefficient is +0°737, higher than for the A A plot. 


Fic. 2.--Red shift 4 as a function of log equivalent “ ergs"’. Open circles : lines of 
length less than 5100 A. Barred circles lines of wave-length between 5100 
Filled circles : lines of wave-length greater than 6100 A 


wate 


1 and 6100 A 


2.2. Comparison with the A.O.-B.S. results.—The existence of a relationship 
between red shift and Rowland number in individual spectral regions has been 
several times pointed out in the A.O.-B.S. papers (7, 9, 10). ‘The use of Rowland 
numbers to express line strengths obscures the fact that there is a relationship 
common to all regions, since the strength of lines of the same Rowland number is 
now known to vary considerably in different regions of the spectrum. Burns 
regarded the intensity effect in the red shift as an insuperable obstacle to accepting 
the predicted relativity shift as established for the solar lines (8). ‘The intensity 
effect is also to be found in the work of St John (18) who interpreted it as 
originating from motions in the different photospheric levels. What is new 
in the present discussion is that the relationship we have found is common 
to all wave-length regions, and within the observational scatter is a linear one 
between A/A* and the logarithm of the absorbed energy. It is obviously 
important to see if this relationship is also to be found in the earlier results 
The A.O.-B.5S. papers and the paper by St John both contain large numbers 
of measured red shifts with which the present results may be compared. St John's 
material is, however, more heterogeneous than that of the A.O.-B.S. workers: 
his solar wave-lengths are means of grating measures by himself and interfero- 
metric measurements by Babcock, and the vacuum-arc wave-lengths are obtained, 
some using a vacuum arc and some from the air arc with corrections to a vacuum 


al 
| 
A 
10 
% ° 
oO 
- 
4 16 16 19 20 24 22 
\ 


416 M.G. Adam 


source. Burns, Meggers and Kiess give red shifts from wave-lengths inter- 
ferometrically measured in the solar spectrum and in the vacuum-arc spectrum. 
The same experimental arrangement and procedure are used for all regions of the 
solar spectrum and the same technique for the vacuum-are spectra. Their 
material should undoubtedly be homogeneous and for the purpose of this 
comparison we have accordingly used the A.O.-B.5. results. ‘The solar spectrum 
in the A.O.-B.S, work is that of the integrated disk; there are possible small 
differences between integrated-disk wave-lengths and centre-of-disk values (25), 
so that we should not necessarily expect quantitative agreement with the Oxford 
results. 


y 


s a function of log equivalent *' ergs 


Fic 4/* as a function of log equivalent “' ergs"’. Symbols as in Fig. 2 


From the tables on pp.112, 130, 142 in Publications of the Allegheny 
Observatory V1, VII, VIII, all the unblended lines of neutral elements were 
listed for which a red shift is given in these tables and an equivalent width in 
Allen’s tables. ‘The A.O.-B.S. lines of weight D were omitted since the 
estimated uncertainty in the solar wave-lengths is more than 1 part in 3 x 10°. 
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The equivalent widths were taken from Allen's tables as follows: 3978-4300A 
(4), 4277-6600 A (2), and 6600-6768 A (3). This procedure gives 158 lines 
with measured red shifts and equivalent widths, namely 131 Fe, 12 Ca, 6 Ti, 5 Ni 
and4Mn. ‘The distribution in wave-lengths and in equivalent widths is shown 
in Table VI. ‘These results may be analysed in exactly the same way as the 
20 entries in Table V. The equivalent widths given in milliangstroms were 
again expressed in equivalent “ ergs”’. 


Taare VI 

Lines providing measured red shifts from the A.O.-B_S. tables 
Wave-length No. of Range of 
range in A lines equivalent widths 
3977-4500 39 82-749 mA 
4500-5000 21 65-166 
5000-5500 19 76-132 
5500-6000 14 26-137 
6ooo-6500 59 21-245 
6500-60768 6 74-130 


Our first comparison is for identical lines in the A.O.-B.S. and the Oxtord 
results. Red shifts are not given in the A.O.-B.5S. tables above for the lines 
5078-97, 5079°2, 6572°8 and 6575-0 A; our comparisons are therefore with the 
remaining sixteen lines of Table V which are common to both series of measures. 
For these sixteen lines correlation coefficients are shown in ‘Table VII together 
with the corresponding Oxford results. ‘The abscissae in each of the plots to 
which the coefficients refer are log equivalent “ ergs’’, and the respective ordinates 
are given at the top of the columns. Fig. 5 shows the A/A* scatter diagram for 
the 16 A.O.—B.5S. observations, together with the regression lines 


The signifi 
cance level, for the correlation coefficient +0°643 for the 16 lines, is between 
oor and o-oor, 


Vil 


TABLI 


Correlation coefficients from the A.O.-B.S. and the Oxford observations 


A A/A 


+ 0°50 0°643 
16 lines 


(xtord 


20 lines 3 


The A.O.-B.5. results thus confirm the observational relation deduced from 
the Oxford measures, with a rather larger scatter of the observed points. 

It may be objected that the results we have obtained are due to the individual! 
peculiarities of the lines we have selected, and the confirmation from 16 A.O.-B.S 
results does not meet the point. In order to test the generality of the relation 
between A/A* and line strength we have examined the whole of the A.O.-B.S 
results, which multiplies the number of lines by nearly 10, and extends the 
wave-length range by 1000A. Fig. 6 shows for the 158 points the relation 
between A/A* and log equivalent “‘ergs’’. We see that go per cent of the lines 
fall in the range of line strengths we have already covered. In order that th: 
few stronger lines, which greatly increase the range in log equivalent “ ergs” 
should not influence the result we have calculated the correlation coefficient 
and regression lines for log equivalent “ ergs” only up to 2:2, since the strongest 


line in ‘Table VII falls in the interval 2:1-2:2. ‘lhe calculated correlation 
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coefficient is +0°504, which for the 141 points used indicates for our empirical 
relation a significance level better than o-oo1. ‘The regression lines are shown 
in Fig. 6 (by a broken line beyond the range for which they were computed). 
‘The y on x regression line, which is based on the assumption of errors only in 
\/A*, is probably the significant one in this case ; it will be seen that this line (the 
less steep one in Fig. 6), which has been computed for the weaker lines, also fits 
the strong lines very well. The smaller correlation coefficient for the 141 points 
indicates a larger scatter for them than for the 16 A.O.-B.S. points in Table VII, 
possibly a consequence of unsuspected blends in the solar lines, 


4]: 


as a function of log equivalent © ergs"’. A.O.— B.S. measurements for identical 
lines. Symbols as in Fig. 2. 


The strongest lines in Fig. 6 nearly double the range over which we may 
examine the relation between A and line strength; they should moreover be less 
affected by blends than the weaker lines. In order to make use of them to the 
greatest advantage we have used the g strongest lines from Fig. 6 together with the 
16 identical lines from ‘Table VII which we believe to be reasonably free from 
blends. For these 25 lines (using all A.O.-B.S. results) the plot of A/A against 
log equivalent “ ergs’ gives a correlation coefficient of +0°785, and the plot of 
A A* against log equivalent “ergs” a correlation coefficient of +0882. Fig. 7 
shows the latter plot with the computed regression lines; the dotted line is the 
regression line (y on x) for the A/A* plot of the 20 Oxford lines. 

2.3. Conclusion.—-Yor the change of red shift with line intensity, we have 
found a well-defined functional form applicable at all wave-lengths. ‘The 
differential shift between weak and strong lines is larger than the absolute shift 
we should expect either from the relativity theory or from collisional shifts. 
Observationally the intensity effect is the most important result and the one we 
should try to interpret first. The quantity \/A* appears to be the most 
significant one in the increase of red shift with line strength, and this suggests 
that the solar red shift arises from the displacements of atomic energy levels which 
vary as the ratio of the line absorption coefficients. In so far as the abundances 
and excitation potentials of the different elements involved are comparable, the 
displacements vary as the ratio of the atomic absorption coefficients. 
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In view of the observational results we must conclude that our approach to 
the problem of the solar red shift should still be largely empirical. It is essential 
to determine as precisely as possible the intensity and wave-length variations in A. 
In particular we require more information about the red shifts of the stronger 
lines. The behaviour of the few stronger lines which have been included in our 
discussion is in good agreement with that of the weaker lines, but these stronger 
lines are all confined to the violet region of the spectrum and the strongest of them 
is only Rowland number 10. Observations of strong lines in other spectral 


on 


60 


22 24 28 


Fic, 7.-A.O.~B.S. measurements, 4/* as a function of log equivalent “ ergs’ for identical lines 
and stronger lines, Symbols as in Fig. 2 

regions are needed, and moreover the red shifts measured by Miss Warga (22) 
for the strong magnesium lines are anomalous. For the wave-length variation 
we have not, in the visible region, a very large range from which to detect the 
rather small wave-length change in A. We can only improve the situation by 
reducing as far as possible the accidental errors in the measured red shifts. Some 
of the observational scatter in the results we have discussed arises, no doubt, 
from unsuspected blends in the solar spectrum, possibly from the peculiar 
behaviour of individual elements and factors other than wave-length and line 
strength which affect red shifts. The smaller scatter of the Oxford results 
(Table VII) suggests, however, that accidental errors of measurement are 
responsible for at least some of the large scatter of the A.O.-B.S. results in Fig. 6. 
The method of circular channels apparently gives increased accuracy of measure- 
ment, and the paper to follow suggests improvements in technique which should 
give wave-lengths with an accuracy as great as the nature of the lines allows. 
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INTERFEROMETRIC MEASUREMENTS OF WAVE-LENGTHS 


Ill. A NOTE ON OBSERVATIONAL TECHNIQUE IN THE 
METHOD OF CIRCULAR CHANNELS 


M. Adam 


(Communicated by the Director of the University Observatory, Oxford) 


(Received 1955 August 29°) 


Summary 
While the method of circular channels has unique advantages for wave- 
length measurements, it is more sensitive to errors due to irregularities in the 
interferometer plates, than is the usual method of parabolic channels. ‘The 
experimental arrangement is examined from this point of view, and it is shown 
that by slight changes in the original set up errors in the absolute wave- 
lengths due to plate irregularities, can be reduced to less than o’oo1 A. 


‘To obtain interference fringes of absorption lines in a continuous spectrum 
relatively free from inter-order scatter, ‘Treanor (§) has shown that the slit pat- 
tern of the auxiliary spectroscope should be considerably narrower than the range 
of the interferometer. ‘The high spectroscopic resolution thus required is neces- 
sarily associated with high dispersion, and this in turn makes it impossible, in 
any conventional arrangement, to photograph simultaneously the fringes of the 
cadmium primary standard line and the fringes of absorption lines in a different 
spectral region, except by the use of an auxiliary camera. ‘Treanor’s interfero- 
metric method of circular channels may, however, be adapted (1) to overcome 
this serious limitation in measuring wave-lengths, and all the wave-lengths de- 
termined at Oxtord since 1949 have been measured directly against the primary 
standard, In this sense the solar wave-lengths fulfil for the first time the most 
essential requirement of secondary standards. 

Simultaneity of exposure of the standard and the unknown lines requires, 
however, photographic equality of the illumination from the two sources. In 
the experimental arrangement so far used, the equality has only been achieved at 
the cost of using a larger area of the interferometer plates for the cadmium light 
than for the solar (or arc) light. Fig. 1 (a) shows the optical trains for the cad- 
mium and solar light as used since 1952; the arrangement differs very little 
from that already described in Paper I (1). A diaphragm D, in front of the cad- 
mium lamp is imaged by the lens L, the field lens F, and the collimating lens C 


as a circle 13 mm diameter between the interferometer plates at P. Similarly 


the field lens F, and C image the spectroscope collimator (2) at I as a rectangle 
3°38 mm wide by 0°77 mm high in a plane 5-4cm behind the interferometer 
plates. Each point on any spectrum line uses a patch of this size on the inter- 
ferometer plates, but the patch is shifted laterally and vertically depending on the 
point on the field lens from which the light comes. 


The shaded rectangles in 
* Received in original form 1955 May 5s. 
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Fig. 1 (b) show the total plate areas used by lines in positions 1 and 3 on the field 
lens, while the circle corresponds to the diaphragm D, (14 mm in diameter) in 
front of the interferometer and represents, closely enough, the plate area used by 


the cadmium light. We ignore the spread by multiple reflections at this stage. 


It is thus clear that not only is the plate area used by the cadmium standard 
larger than that used by the solar lines, but also different lines do not use identical 
areas on the plates. Any departures from plane parallelism in the plates will 
introduce errors in wave-lengths. 


)CADMIUM LAMP ‘S 
= D> 


Ah 


2 
Fic. 1.——-Interferometer illumination for solar and cadmium light. 


Errors introduced by plate irregularities vary inversely as the length of 
spacer used. ‘The general agreement between the wave-lengths measured with 
a 2°5 mm spacer and a 5 mm spacer in Paper | indicated that any such errors 
were within the uncertainty of the observations. When, however, a 1°5 mm 
spacer was brought into use at 6500 A with the experimental arrangement above, 
it was found that the wave-lengths measured with it were systematically shorter 
by approximately 0-004 A than those measured with the 2-5 mm spacer, ‘The 
differences were shown equally by the are and solar lines, they were not, appar- 
ently, dependent upon line strength and, for a region so close to the primary 
standard, cannot be ascribed to uncertainty in the chromatic variation of phase 
change, or camera focal length, or displacement of the pattern centre. We con- 
cluded that they arose from non plane-parallelism of the interferometer plates. 

If our explanation is correct it indicates that, as the plates are normally set 
up, the air gap ft between them is greater in the middle than at the edges. In 
determining wave-lengths we use 2t/(p + 


«), p and « being the integral and frac- 
tional orders of interference. 


2t is derived from the cadmium pattern and re- 
fers to the full circle in Fig. 1(b), but « refers only to that area of the plates 
used by the solar line. In order to get a quantitative measure of the plate ir- 
regularities it was decided to photograph the localized fringes formed between 
the plates when they are set at a small angle. ‘The experimental arrangements 
necessary to secure sharp fringes which accurately contour the surfaces are 
pointed out by Rasmussen (3) and ‘Tolansky (4), and Rasmussen has used this 
method to test interferometer plates. ‘The plates were set up in the usual mount 


with three pieces of mica 0-012 mm-o-o1s5 mm thick placed in the positions 
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usually occupied by the spacer studs. ‘The wedge angle between the plates was 
adjusted by means of the interferometer springs to give approximately one fringe 
per mm in cadmium green light. By micrometer measures across the fringe 
photographs, values of 2¢ (with an arbitrary zero) were deduced at 642 points 
over the plate area and these have been used to draw the “ contour map”’ of the 
air gap which is shown in Fig. 2. The figure shows the irregularities of the 
plates as a pair in one particular orientation when they are set out as parallel as 
possible, and 2¢ is arbitrarily set equal to zero on the contour line indicated. 
‘The interval between the contours is 1/100 of a wave-length of cadmium green 
light ; between the stronger lines it is 4/20 corresponding to A/ 40 on each plate. 


FG, 2.~-Variations from planeness of the pair of interferometer plates. 


The plate curvature is of the form we have predicted and the variations in 2¢ 
are comparable with those to be expected from the tolerances given by the makers. 
The curvature is, however, rather more than sufficient to account for the wave- 
length anomalies we have found. It is possible that the curvature is accentu- 
ated in setting up the interferometer with the very thin mica separators and by 
the unequal spring pressure which is used to set the plates at a small angle (3). 
The stud positions at which the pressure is applied are indicated by 8,, S,, S, in 
Fig. 2 and in Fig. 1 (6). If the curvature is in part introduced by spring pressure 
it may well vary with spacer length. 

The elimination of errors due to plate irregularities.-\n order to test the plates 
under normal conditions an Osira cadmium lamp was used as the light source 
for the spectroscope. ‘The wave-length of its red line was measured exactly as 
tor a solar or an arc line in terms of the cadmium red line produced in the side 
arm of our experimental arrangement. By changing the position of minimum 
deviation the spectroscope cadmium line can be placed in any position across the 
field lens F,. We specify the position by the distance in mm on the photo- 
graphic plate from the pattern centre, positive on the violet side and negative 
on the red side. 
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The results of three such sets of measures using the 1-5 mm spacer are shown 
in Table I below, the line positions correspond to the limiting ones normally 
used in measuring wave-lengths. A is the “error” in the wave-length of the 
spectroscope cadmium red line. For the first two sets of measures in Table I the 
experimental arrangements were those of Fig. 1 (a). 


If these differences AA arise, 


Tasie 


| Line position 


AA 


Differences in Cd red measured wave-lengths 


\A 
mm (1) (1) Coine : patches 
| 
16 0°0049 0°0040 


6 0°0044 0°0039 0°'0004 


0'0030 0°0034 0'0004 


0°0046 0°'0007 


as we suppose, from the different areas on the plates used by the two cadmium 
sources, then we should expect to remove them by making the areas coincident. 
By using a 33 cm lens for the field lens F,, the spectroscope collimator is imaged 
between the interferometer plates. This focused patch (3:8 mm x 0-77 mm) 
is shown in Fig. 1 (6) by the dotted rectangle. ‘The aperture in the diaphragm D, 
before the cadmium lamp, was also made rectangular 3:8 mm = 0°77 mm, and 
this is imaged as before between the plates (approx. 1: 1) to coincide with the 
spectroscope collimator image. The last column of ‘Table I shows the values of AA 
obtained with the 1-5 mm spacer with these coincident patches. As we expect, 
the large fairly systematic differences in the wave-lengths have disappeared and the 
remaining small differences are probably within the uncertainty of their deter- 
mination. 


Using the same interferometer plates a similar test was also made of the errors 
to be expected with a 2-5 mm spacer and the experimental arrangement of Fig. 1. 
We should expect smaller differences than were shown by the 1-5 mm spacer, 
but the AAs obtained were in fact all under o-oo1 A and probably not greater than 
their estimated uncertainty. ‘These results certainly suggest that plate curva- 
ture is to some extent introduced in setting up the interferometer. ‘The precise 
amount no doubt varies on different occasions and presumably depends also on 
the relative orientation of the plates. 

With interferometer plates of the quality of the pair examined above (‘Table 1) 
coincidence of the spectroscope and standard light patches is sufficient to secure 
absolute wave-lengths correct, as far as plate errors are concerned, to 0-oo1 A, 
But even with coincident patches, in the method of circular channels multiple 
reflections carry the spectroscopic and standard light into different regions on the 
interferometer plates, and with a less good pair of plates, which were also used, 
wave-length errors of 0-001-0°002 A could still be caused on outlying lines (2-5 mm 
spacer). We have taken care of this spreading-of-beams effect by modifying the 
diaphragm on the field lens F, ; side arms have been added to the simple cross 
which give on the photographic plate two cadmium red “ lines” more remote in 
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the field than any solar or arc line being measured. ‘The diaphragm on F, is cor- 
respondingly modified so that there is no overlap of the solar and cadmium light 
on the plate. From the “ errors” 


in the cadmium lines the small corrections to 
the solar or arc wave-lengths can be deduced. With this procedure the absolute 
wave-lengths should be free (to less than o-oo1 A) from errors introduced by plate 
irregularities. ‘lhe reduction in the plate area used by the cadmium light gives 
a corresponding reduction in light flux; this may be compensated by increasing 
the red light reflectivity of the beam-splitting cube and, if necessary, the intensity 
of the solar light must be reduced by filters. 

We must then point out that the absolute wave-lengths given in Papers I and 
I] (1) may be slightly in error due to the effect of plate irregularities. On the 
basis of subsequent tests on the two pairs of plates used, we estimate the maxi- 
mum value of such errors at 0-002-0°003 A. Other lines of evidence suggest, 
however, that it is unlikely that such a large effect is in fact present in the values 
given. In that the solar and vacuum arc wave-lengths have been determined 
with identical experimental procedures, their values will be affected by plate 
irregularities in exactly the same way. Consequently the red shifts, used for the 
discussion in Paper Il, have no systematic errors due to plate irregularities, just 
as they have none due to the uncertainty in the chromatic variation of phase 
change, or camera focal length, or displacement of the pattern centre. 


Umversity Observatory, 
Oxford 
1955 August 27 
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ON THE FLEXURE OF MERIDIAN CIRCLES 
R. a’ E. Atkinson 


(Communicated by the Astronomer Roya! 
(Received 1955 March 25) 


Summary 

The methods which have been used or suggested for determining meridian- 
circle flexure are without exception unsatisfactory, Very few of them can 
detect a second harmonic, and many of the results actually obtained are 
unacceptable astronomically. The paper analyses in detail the effects of 
uncertainties in the flexure function, and also in the refraction, on the 
solar corrections, giving a more rigorous treatment than has so far been used ; 
this rigour, and also an improvement in the methods of determining flexure 
observationally, both seem desirable if any further advance is to be mack 
in the precision of the results. A new method of directly determining all 
harmonics of the flexure function is proposed, and the necessary equation of 
condition is set up for a particular case which seems the most surtable one 


Introduction.—\t is universally recognized that the flexure of a meridian circk 
is one of the most difficult of instrumental errors to determine. As was first 
pointed out by Bessel (1), a value applicable to the horizontal position can be 
obtained immediately, by observing the two collimators after they have been set 
on one another; but the natural assumption that the flexure at other points on 
the circle will be this “ horizontal flexure’’ multiplied by the sine of the zenith 
distance frequently proves irreconcilable with the astronomical results, Various 
methods of determining the true form of the flexure, as a function of zenith 


distance, have been suggested, and some have been tried very extensively, 


but none seems free from fairly serious objection; since any uncorrected error 
will enter directly into the declinations of every star catalogue based on the 
instrument in question, and will also affect the corrections to the Sun's « lements, 
the position is definitely unsatisfactory. ‘The present paper reviews the different 
attacks which have been made on this problem, discusses in detail the effects 
on the selar elements which different correction functions will produce, and 
proposes a new method of measurement which appears to be free from the 
objections that have been ¢« xplicitly noticed, 

In order to avoid misunderstanding, it is first desirable to clarify a pont of 
terminology. Flexure is sometimes spoken of as positive if it “ increases the 
zenith distances’’, whether these are north or south, and the increase is ther 
meant in the numerical sense and not algebraically; it is, however, simpler in 
some ways to regard the flexure as an ordinary transcendental function of z, 
with z positive north of the zenith (say), and to agree that in all cases the observed 
Z.D. is to be increased algebraically by the amount of the flexure in order to get 
the true Z.D. Such flexure functions as f=/,sin2 must then be called antu- 
symmetrical, not symmetrical, about the zenith, although they correspond to an 


428 R. Atkinson Vol. 115 


arithmetically (or even ‘‘ geometrically”) symmetrical effect; a flexure varying 
as cosz will similarly be spoken of as symmetrical about the zenith, although it 
deflects the zenith point, and although arithmetically it increases the Z.D.’s on 
one side of the zenith and decreases them on the other. In the present case, 
the algebraical approach seems definitely more convenient than the geometrical, 
and we shall follow this convention throughout. 

The nature of“ flexure” .—-\t is well known that if the eye-end and objective-end 


of the telescope suffered identical flexures, their combined astronomical effect 
would be zero. 


The two tubes can be made nearly identical, especially if a 
draw-tube for focusing is dispensed with (2); if their lengths are exactly equal, 
and if the “‘length’’ (flange to second nédal point) of the lens cell is equal to the 
“length”’ (flange to focal plane) of the micrometer box, and if the instrument is 
also well balanced (as it normally is, with some delicacy) both with and without 
the eye-end and cell, interchanging the latter will eliminate the effect of any 
differential flexure of the bare tubes which may still exist, at least if this occurs 
only as a first harmonic (or as a sum of odd harmonics) in the circle readings ; 
but the objective and micrometer themselves are of course quite dissimilar, and 
any internal movement in either will have the same astronomical effect as a flexure, 


with the added complication that it may possibly be non-elastic, and may even 
show discontinuities 


For this reason, and also because it will not in any case 
eliminate any even harmonics, there seems little point in actually interchanging 
the two ends, and the practice has indeed been fairly generally abandoned, even 
in instruments which readily permit it. 


Personality errors should be eliminated 
by the systematic use of a reversing prism, but if this precaution is omitted they 
also will enter into the “ flexure”’ 


indeed, in at least one programme (3), a 
systematic difference in the zenith distances of near-zenith stars was found, 
depending on whether the observer lay head-north or head-south, even though a 


reversing prism was used, and it thus seems that this device cannot be relied on 
implicitly, 


Certainly there is no reason why such personality effects as are not 
eliminated should vary as sin z, or should not have a discontinuity at the zenith. 

The problem is (as it seems) further complicated by the possibility of flexure 
in the circle; the circle readings have to be combined with the Z.D. micrometer 
readings in most (but not all) flexure tests, and of course in all astronomical 
observations, and even if one assumes that circle flexure will vary as sin(z—y) 
there may be no evident reason why y should be small. 


If it is not, the combined 
flexure of telescope and circle will have terms in both sinz and cosz; in the 
Washington g-inch T.C., the “circle flexure’ was found to be considerably 
larger than the other (3), and it seems that the resultant might have any phase. 

It is however a little surprising that observable circle flexure should exist at 
all. ‘The true flexure of an elastic but symmetrical circle, firmly attached to the 
main axis all round but free to droop under its own weight, would produce only 
the same effects as slight constant displacements of the various microscopes, and 
would thus be completely irrelevant even if only one microscope were used; the 
microscope displacements would not even be noticed, since the microscopes are 
initially adjusted on the basis of the circle divisions as actually read. If the 
circle is not equally stiff all round, or is not attached to the axis in a uniform and 
symmetrical way, detectable circle flexure may of course exist; but if local 
weaknesses are postulated, they seem in general more likely to affect the differences 
of the microscope-readings than the sum of all four (or six), and on the whole 
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one would then be inclined to expect observable circle-flexure to occur as a fairly 
high harmonic of z, not merely, as the first, which is what has been found in 
practice. A single local excess of weight would be effective, but it seems somew hat 


improbable that any manufacturing flaw of this sort could be large enough 


. There 
is, however, 


one case where an apparent first harmonic would be found alone - 
a uniform circle which was effectively attached to the axis 


at only two points, 
with one of those attachments stiffer th 


an the other, would suffer rotations about 
some virtual pivot (the resultant of these two partial constraints), and if the 
constraints were elastic these rotations would be proportional to the sine of that 


pivot’s Z.D., and would be equivalent to translations plus similar rotations 


about 
the axis 


the translations would disappear with two or more microscopes, and the 
rotations would be interpreted as first-order flexure. 
in the fact that the normal method of determining 
** movable’ 


The importance of this lies 
circle flexure is to have the circle 
' (i.e. with respect to the telescope), to rotate it every other year (say) 
to a new position, and to analyse the measurements of horizontal ( 
also vertical) flexure in terms of these circle positions ; it may, 
harder to guarantee a fully adequate and symmetric 


and perhaps 
however, be much 
al connection between circle 
and axis when this provision for rotation is made than when the circle is deli 


verately 
constructed as a “‘fixed”’ one. 


It would thus not be altogether surprising if 
first-order circle flexure were an error introduced mainly by the provision made 
for detecting and measuring it, and a well-designed “ fixed” circle might in that 
case be preferable to any attainable excellence in a movable one 

Methods for determining flexure. (a) Bessel’s method. 
flexure of circle and telescope is what is required 
be to determine this directly 


Since the combined 
astronomically, the ideal would 
at a number of different zenith distances. In 
principle, as again Bessel first pointed out, it would be possibl. 
collimators in some other position than the horizontal. or in 
positions, to determine the flexure for 
as a function of 


to set up a pair of 
a series of such 
each position, and to express the results 
This method could not be employed, as it stands. in an 
ordinary pavilion, and does not appear to have been explored at all; 


but it would 
seem reasonably practicable if the lower of the 


two collimators were replaced by a 
plane mirror and the upper one were fitted with a Bohnenberger eyepiece. 
However, even if suitable arrangements could be made, the collimator method 
does not really give the true astronomical fle xure in any 


case, but only the 
difference of two values 180 apart ; 


thus however many collimator positions one 
uses it cannot detect any even harmonics in the flexure func tion 
it cannot detect a second harmonic. 
pointed out previously.) 


, and in particular 
(This weakness seems not to have been 
If it can be assumed that there are no harmonics above 
the first, there is no real need for a large 


number of collimator positions after all; 
a single pair of vertical collimators, in 


addition to the horizontal pair, is sufficient 
to determine the complete function, and the lower one of the vertic 
then be replaced by the ordinary mercury-pool (4). 

‘The vertical flexure” 


al pair can 


, 1c, the flexure (including that of the « ircle) evaluated 
by the use of vertical collimators or their equivalent, is not really 


material in a 
modern instrument, since it is eliminated by reversal 


Indeed, if the flexure 


function is expressed as a Fourier series in z, this is true of all the cosine terms in 


the complete series. If there are no sine harmonics above the first, it must thus 


be sufficient to determine the horizontal flexure /, by means of the 
reverse the instrument regularly, and assume that f=/f 


collimators, 
rotating the circle 
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in addition should be unnecessary, so far as this question goes, whether there is 
circle flexure or not. If the circle is never moved, many of the corrections 
discussed below can also be obtained with greater weight; and since the only 
information which is obtained by moving it is information that is not wanted, 
it would seem better to keep it in one position even if it is designed to be movable. 
(b) Loewy's method. Methods of test of this kind, i.c. employing only 
collimators, mirrors, etc., but not stars, may be classed as “ purely instrumental ”’ ; 
there may also be “ purely astronomical’’ methods (or, at least, checks); and there 
are methods which we may call mixed “ instrumental-astronomical”’ ones. 
‘lwo other purely instrumental methods seem worthy of mention here. The 
first 1s due to Loewy (5); a meniscus lens, with its concave radius of curvature 
equal to half the length of the telescope and its focal length equal to half that 
radius, was very firmly mounted in the central cube, with the concave side towards 
the eyepiece. Using the concave side as a mirror, a real image of the Z.D. micro- 
meter wire could be seen in the same focal plane, and the setting for coincidence 
noted; the variation of this setting, as the telescope was pointed to various 
zenith distances, gave the flexure of the eyepiece half of the telescope. Using the 
lens property, a mark on the telescope objective could also be imaged in the focal 
plane; and the variation of Z.D. micrometer settings on this, as the telescope 
was moved, gave the sum of the flexures of the two halves. ‘The difference of 
the two then followed. The main difficulty, that of determining the flexure 
of the auxiliary lens and mounting itself, was solved (it was claimed) by adding 
weights to this mounting, repeating the observations, and “eliminating” the 
flexure due to the original weight. If the difficulty certainly can be got over 
in this way, the method appears sound; but in the case of the reflection observa- 
tions it is tilts of the mounting, even more than bodily displacements, which 
are to be feared, and it might be difficult to be sure that the loads were always 
quite suitably applied, from this point of view. Even a small miscalculation 
here can have very serious effects indeed. ‘The results actually obtained were 
claimed to be very consistent, and it is interesting to note that they did indicate 
the presence of an appreciable second harmonic; the flexure was found to be 


+1"04sin +017. C08 + 0"-308in 22 — COs 25. (1) 


‘The method evidently yields only the telescope flexure, without the circle 
flexure. However, it was argued that the latter would be expected to be zero; 
ana the horizontal collimators, which, of course, measure the combined (horizontal) 
flexure, gave +1°:15, in fairly good agreement with the + 1°04, above, for the 
telescope only, ‘This is perhaps the only case in the literature where an important 
second harmonic was directly determined; but the possibility that some flexure 
of the mounting was not fully accounted for by the process of varying the loading 
cannot be entirely ruled out. The method is inconvenient, since the meniscus 
prevents all astronomical observations while it is in place, and it does not seem 
to have been adopted elsewhere. 

(¢) Bonsdorff's method. ‘The other purely instrumental method is due to 
Bonsdorff (6); a plane mirror was supported in front of the objective, by a 
mounting designed to guarantee that flexure of the mounting would move the 
mirror strictly parallel to itself, while flexure of the telescope would not affect it 
at all, Auto-collimation measurements were then made with this mirror, in a 
variety of telescope positions, so as to derive the complete flexure function of the 


pa 
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telescope. It was stated that adding extra weights to the mirror verified that it 
was functioning correctly, the reading at any one Z.D. being unaffected, as it 
should. It was also stated that the results indicated a variation purely with 
sin z, but this appears to be incorrect; the treatment given to the individual 
readings would necessarily have removed any term in cos =z, and if these readings 
themselves are examined they indicate not only that there was such a term, but 
that it was actually about three times as large as that in sins. It is unfortunate 
that no comment was made on this; it seems rather improbable that so large a 
cos z term could be real, for in this case also the circle flexure evidently does not 
enter. ‘The method is believed to have been tried at several other observatories, 
but to have been found difficult in most cases. It is clear that perfection of 
design, and high precision in manufacture, are essential if the mirror’s move- 
ments are to be sufficiently parallel; a slight difference in stiffness as between two 
members that were nominally identical, or a slight difference in friction as between 
two joints that were both nominally free, would be enough to impair the action, 
and it would perhaps not be very surprising if two superficially indistinguishable 
mountings gave somewhat different results. 

With both these methods, if the results did not agree with those obtained from 
the collimators, it would be difficult to know what to do with them. ‘The 
difference might indeed be ascribed to circle flexure ; but since the circle readings 
must be used in the astronomical work itself this would mean that the collimator 
results had to be relied on after all, and the others rejected (so far, at least, as the 
first harmonic is concerned), unless the form of the circle flexure could be 
discovered independently. 

(d) The “R. and D. method”. By “mixed” instrumental-astronomical 
methods we understand those in which stars are observed, but nothing is assumed 
about their positions except the fact that for any one star the (mean) declination 
is constant; its absolute value does not enter. In addition, an instrumental 
arrangement is made which permits the same star to be observed in two or more 
different positions of the telescope, related to one another by some condition 
which is known with precision ; and a comparison of the circle readings actually 
obtained, with those to be expected on the basis of the known condition, gives a 
condition which the flexure must satisfy. ‘The best-known of these methods is 
that of observing the star both directly and by reflection in a mercury pool. ‘The 
precisely-known condition in this case is that the sum of the two observed circle 
readings, together with the appropriate flexures, must be exactly 180° ; moreover, 
half the difference of the two readings should give the star's altitude, without 
relying on the ordinary nadir pool at all. Since the nadir observation is of a 
different kind from all observations of stars (and is also made in a different position 
of the observer, and of the micrometer), this possibility of eliminating any 
systematic error that may affect it appears at first sight valuable 

The method was tried extensively at Greenwich (7), but the results were never 
very satisfactory. It seems probable that one reason for this lay in the unsuitable 
nature of the pavilion in this particular case (8). (‘The aperture of the Airy T.C. 
pavilion is only 3 ft. wide, and is 16 ft. above the axis of the instrument; and the 
instrument itself is also partly sunk in a pit, so that the mercury pool is in many 
cases below ground level, which may perhaps add to the refraction dangers inherent 
in the pavilion.) ‘The method appears to have been found more successful at 
Washington (9), though a mercury pool can of course be used with a wide 
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pavilion-aperture on very still nights only. It has not been employed recently at 
all; and actually, even in ideal circumstances, it is defective in principle: it can- 
not detect or eliminate any component of flexure which is antisymmetrical about the 
horizon. This weakness does not appear to have been pointed out in the literature, 
at least in this precise form; it seems to deserve explicit consideration. 

Let c be the circle reading, reckoned from the zenith, positive for Z.D. north, 
when the telescope is set on a star at zenith distance z (as refracted). (For 
convenience, we speak of the combination of Z.D. micrometer reading and circle 
reading as the “ circle reading c”’.) Let f(c) be the flexure at the circle reading c, 
positive if its effect is such that the true Z.D. is greater (algebraically) than c by the 
amount of the flexure ; so that the zenith distance of a star (including refraction) is 

z=c+f(c). (2) 
Let cy be the reading for a given star if it is observed directly, and c, the reading if it 
is observed by reflection in the mercury pool; cy is supposed corrected for the 
necessary small change of latitude of the pool, though this precaution has perhaps 
not invariably been observed in practice. 

Then we have 


s=¢n+f(cp) (3) 
and z=1—Cy—f (cp). (4) 
Let f now be divided into a portion, /,, symmetrical about the horizon, and a 
portion, f,, antisymmetrical about it; then (ep) and (cy) (ep). 
Inserting these in (3) and (4), and adding and subtracting the results, 

23 = ly (Cp) (5) 
and t+ Cy 7+ 2fg(Cp). (6) 


Equation (5) shows that any systematic error affecting the zenith (or the nadir) is 
indeed eliminated when the true Z.D. is inferred in this way (since it will enter 


equally into cy and cy); and equation (6) is an equation of condition for determining 
fa; we may assume for example 


f[a(c) = X [a,, cos + b,,,,, sin(2n + t)c] (7) 


using 48 many terms as seem reasonable, and on taking a sufficient number of 
observations we shall be able to solve by least squares for the a’s and 4's and so 
determine fg completely. But there is no similar equation of condition for /,, 
and f/, is present in (5) and also, of course, in (3) and (4): thus one cannot say that 
‘the combination of a direct observation with a reflected one eliminates the effect 
of flexure ’’, nor can one infer a total flexure for use with single direct observations 
(which must be the majority), unless /, is known on other grounds. 

The terms which are antisymmetrical about the horizon are, of course, the odd 
cosines and the even sines of ¢ (or, as we may equally well say, of z); neither odd 
nor even sines are cancelled out by reversing the instrument. ‘Thus even if the 
instrument is reversible the mercury-pool method will leave undetected and 
operative any term in sin 2s (or in sin 42, etc.) which may exist; and in addition it 
will fail to detect or eliminate any term in cos z (or in cos 32, etc.) if the instrument is 
not reversible, as the Airy T.C. itself was not. Higher harmonics than the second 
are perhaps unimportant, at least unless they occur in the circle rather than in the 
telescope, or are used to describe discontinuities, but it seems a little dangerous 
to assume that even the second harmonics are absent. It is an unfortunate 
coincidence that neither the collimator method nor the mercury-pool method can 


= 
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say anything about a term in sin 22, nor can such a term be eliminated by inter- 
changing the lens and eye-end, or by reversing the instrument. 

(e) The “astronomical method” (flexure and refraction). In practice, no 
matter what methods are used for determining flexure to begin with, they are 
fairly often discarded again at a later stage of the same programme, the final choice 
of a flexure function being then based exclusively on the astronomical criterion that 
the distance from the pole to the equator is go". This is, of course, sufficient to 
determine one parameter only, and the flexure is therefore assumed to be of the 
form f=/,sinz, with f, so adjusted that the astronomical criterion is satisfied. 
All preliminary determinations of flexure, e.g., with the collimators, are thus 
reduced to the status of mere devices for keeping the arithmetic compact; an 
accurate determination is pointless, except perhaps as a check that no unexplained 
change has occurred, and any arbitrarily assumed value of f, near the final one 
would be sufficient. Occasionally, when the original determinations of horizontal 
flexure and of refraction led to positions of the pole and equator which, together, 
did not satisfy the astronomical criterion, the flexure determinations have been 
kept and the refraction ones amended instead ; but in either case the situation so 
created is rather unsatisfactory. It is always undesirable to throw away apparently 
good observations without discovering why they are wrong, and particularly so 
if it is only done in order to preserve others which may themselves be erroneous. 
The Washington 6-inch results for 1925-41 (10) and 1941-9 (11), which can 
certainly be classed as among the best available, show this difficulty very plainly ; 
at the latitude of Washington, flexure and refraction are definitely hard to separate, 
on the basis of above-pole and below-pole declination comparisons, and when the 
distance from pole to equator proved to need correction it was necessary to 
sacrifice either the preliminary determinations of flexure or those of refraction with 
little to go on. In the former series, the refraction coefficient was changed by 
about 0”: 16 although the probable errors of the rejected values (slightly different at 
different dates) had been + 0"-016; in the latter series, the flexure coefficient was 
changed by 0”-04 while the probable errors of the rejected determinations (12) were 

+0"-03. ‘The conflict in the second case is superficially less, but only because the 
equator-point error which had to be corrected was much less; if it had been as 
large as the other, the flexure coefficient would have had to be changed by 0:24 

In both cases, of course, everything that was known about possible systematic 
uncertainties, as opposed to the visible probable errors, was considered before 
making the decision; but it does now appear at least formally possible that the 
provisional determination of refraction and the provisional determination of 
first-order flexure may both have been substantially correct, and that the error of 
the equator-point may have been due to a second harmonic in the flexure, which 
could not have been detected by the methods used. And at least it is clear that a 
considerable advance in the reliability of the final results might be hoped for, if the 
different possibilities could be more surely separated. It is the case, unfortunately, 
that even though two alternative correction functions may not be clearly distin 
guishable over a moderate range of near-polar stars, they may have appreciably 
different runs over the 47 degrees of declination, centred on the equator, which 
are involved in solar observations; planetary observations necessarily involve a 
somewhat larger range still, and the range of declinations observed in preparing a 
star catalogue may be larger again. It seems worth while now to investigate 
just how serious these effects may be ; we shall see that they can in fact appreciably 
exceed the probable errors claimed for some of the solar corrections. 
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Analysis of the solar corrections.—I\t is well known that if the residuals 
(Obs.-Tab.) of the Sun’s declination are analysed, they yield corrections to the 
Sun’s mean longitude L, to the obliquity «, and to the two quantities A and k from 
which the eccentricity and the longitude of perigee can be inferred; in addition, 
there usually remains a constant term, — Ad, (13), which is commonly regarded as 
‘the error of the equator-point”’. Ad, would be zero if the instrumental errors, 
and the refraction, were correctly known and allowed for ; if it is not zero, and is to 
be removed, this must be done by applying some specifically assumed correction- 
function, such as 

=p +qsins (8) 
(if the first harmonic of the flexure is held to be responsible) or 
\.6=p+qtans (g) 
(if the refraction is); as we have now seen, it may be desirable to consider also 
4.6 =p +qsin 2z. (10) 
In each case, if the latitude has already been so adjusted that declinations above 
and below pole are consistent, the constants p and q are to be so chosen that the 
function vanishes at the pole, i.c. we must have 
\,6 = K(cos¢ — sin z) (11) 
or = K(cot¢d— tan z) (12) 
or \.6 = K(sin 24 — sin 2s), (13) 
where ¢ is the latitude ; and in current practice A is then determined so as to make 
\,.0f (14) 
if \.A(F) is the “ equatorial”’ value of 4,4, i.e. its value when 6=0. Whichever 
of the functions (11), (12), or (13) is selected, it will act as a systematic correction 
to the DS, which characterize a star catalogue obtained with this instrument 
(as compared with others), and this has indeed been the principal purpose of 
deriving it; in addition, as was remarked by Hough (14), it will in general affect 
the solar corrections too, 


The equation of condition in its usual form is 
Ad, + sine cosa AL + sinade + 2 cos sin dAh — 2 sine cos* x cosdAk = Ad (15) 


where Ad are the declination residuals; and since the coefficient of Ak never 
changes sign it is clear that — Ad, 1s not strictly the mean value of these residuals ; 
it is also not the value of the left-hand side of (15) when 6=2=0; and neither 
of these values is strictly the one to which (11), (12), or (13) ought to be adjusted. 
Although, as we shall see, the practical error is not as yet really serious, it might 
easily become so if experimental accuracy improved a little further; we will 
therefore now establish what the rigorous condition for A,6 is. 


We may eliminate 4 from (15) by means of the relation 
sina; (16) 
if we also insert the numerical value of «, we obtain 
Ad, + 0° 398 cosa + sin ade + (0°413 2x — sin 4%) AA 
(0° 389 + 0° 398 cos 2% + 0-009 cos 4a) = AO (17) 


(neglecting higher harmonics); and this form shows explicitly that the coefficients 


are for the most part very conveniently orthogonal if we can assume equal weights 
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at equally-spaced values of « (1§). We shall make this assumption for the 
present; the question of weighting will be further considered at a later stage. 
An additional simplification, so far as A,5 is concerned, follows from the fact that 
any function depending only on z is the same in the autumn as in the spring, at 
the same values of 6; since (for the Sun) « and 4 are connected by (16), we may 
express \.6 by a Fourier series in «, and it must therefore be by one that is 
symmetrical about «=7/2 and —7 2. It can thus only involve the even cosines 
and odd sines of «, and these are all orthogonal to the coefficients of both AL 
and Ah, so that it is only the quantities \d,, Ae, and Ak which can be affected by 
the kind of correction we have to consider. 

If we difference two least-squares solutions, both having identical coefficients 
on the left, and employing “ observational’ quantities on the right which differ 
by the correction function 4,4, the resulting equations are exactly what we 
would get if we solved the equation 

+ sin — (0° 389 + 0° 398 cos 2% + 0-009 cos 44) Ak = AS (18) 
as if it were an equation of condition, to obtain the three additional corrections 
and Ak by the regular least-squares process. In the “corrected” 
solution, these will be added to the original Ad,, Ae, and Ak, which have already 
been obtained by solving the “ uncorrected”’ equation (17); and the rigorous 
criterion for the choice of the magnitude (but not, of course, the form) of the 
correction function 4,4 is that there shall then be no absolute term left, i.e. 


+ Adg = 0. (19) 


Writing \.6 as a Fourier series in « (with the restriction already noted), 
putting 


+ C08 24 + ¢,C08 40+... +d, +d, sin 4a 4 (20) 

in (18), we may readily derive the normal equations 
A. + 0° 3894.2 fo (21) 
bd, (22) 


+ (0°389* + 0°0792)A.k = — 0°389¢, — 0° 199¢,— 0045¢4; (23) 
and separating (21) and (23) 


AR 2°5 (24) 

+ 0°97 + (25) 
Equation (19) thus leads to 

— = Ady (26) 


as the rigorous condition for adjusting the “amplitude” of the selected A,5. 
The “equatorial” value of is given by 


= + + (27) 
and it is, of course, —4,4(E£) which is strictly “the error of the equator-point”’. 
The traditional criterion for adjusting 4,5, namely (14), is clearly not identical 
with (26). 

We may replace (26) by an approximation which climinates the necd for 
explicitly deriving the Fourier series (20). ‘The summer-solstice value of AA, 
A ACS) say, is 


NAS) = Cg ~ + d,~ dy (28) 


30 
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while the winter-solstice value is 
AAW) = + dy + dy; 


(29) 


we thus have 

+ = ~ + (30) 
and the right-hand side of this is very similar to that of (26); indeed, if c, and 
¢, happen to be in the same ratio as the coefficients of cos 2% and cos 4a in (17), 
i.e. if + 0°023¢y, (30) becomes 

b[4,5(S) + = 0°977¢2 
while (26) becomes 

= ~ 0°97 76» 

so that the correct adjustment of 4,6 is obtained when 


+ 4,8(W)] = Ad,. 


‘This equation may be regarded as a “rule of thumb’ 


(31) 
condition for adjusting 
Ad, and it will always be a fairly good one, since c, and 2-3 per cent of c, will 


both be small even if they are not equal; thus, so long as we rely on declination 
measurements only, we may say that the correct adjustment of A.6 is obtained 
by making the mean of its two “tropical” values, and not its equatorial value, 
equal to Ad,. ‘To take a particular case, if we wish to remove a Ad, of +0"-40 
(say) by means of a 4.4 of the form (12), in latitude 40°, use of the equatorial 
value will lead to 


= 0"235 — tan z (32) 
while the tropical mean leads to 
=0":204 — 0°°171 tan z; (33) 


the difference amounts to 0°-03 at the zenith, to o”-o5 at the equator, and to 
o”"-10 até 30 , and the two resulting star catalogues would differ systematically 
by a Dd, of this magnitude. Actually, as we shall see below, these differences 
should normally be reduced by anything up to 50 per cent; but the traditional 
use of (14) still appears a little undesirable, more especially since the amount of 
labour which it saves is trivial. 

It may be useful to tabulate the results of applying the harmonic analysis (20) 
to the three functions (11), (12), and (13), for three representative ¢-values. 
The five coefficients of (20), all normalized to c,—0°977¢,—0°022¢,= +1, Le. 
adjusted to cancel a Ad, = + 1”-00, and the corresponding values of K to insert 
in (11), (12), and (13), are shown in ‘Table I, and values for other latitudes may 
be interpolated graphically. 

Taste | 


Correction (11) Correction (12) Correction (13) 
(sin 2) (tan z) (sin 22) 

d= 4o 40 50 30 40 50 30” 40" 50°” 
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We may use these figures to study the effect of removing any particular Aé, 
by any one of the three assumed correction functions. 
tions in latitude 50° have yielded a Ad, of +0"-30 (after correcting the pole), and 
it is decided that an erroneous first harmonic in the flexure is responsible, then since 
0”-30K = 0"-223 in this case, the correction function is 


If, for example, observa~ 


A.6= +0°:223(cos — sin z) (34) 
and we see from (22) and (24) that the solar corrections previously obtained will 
be modified by the subsidiary corrections 


Aye =0"30d, = —0"-06 (35) 
and Ak = — 0°: 30(2°5 1c, + 0°06c,) 0” (36) 
If we decided to correct the same Ad, by means of a second harmonic in the 
flexure, we should have similarly 


AS = +0°:181(sin — sin 22) (37) 
= +0"'02 (35) 
AR 0”-07, (39) 


while the decision to correct it 


by modifying the refraction instead (16) would 
result in 


+0"-110(cot — tan (49) 
Ne (41) 
A +0"-10. (42) 


‘The differences between these solutions appear appreciable, in view of the fact 
that the internal probable errors claimed for Ae and Ak may be of the order 
+ 0-03 or less, and it is somewhat embarrassing if one has to choose between them 
arbitrarily. 

Table I also allows us to put the effects of an inability to decide between 
different correction functions in another way. For any given latitude, we may 
difference the figures for two different correction formulae (or may make any 
linear combination of all three such that the sum of the three amplitudes is zero), 
and the result wil) be an expression which will leave the separation between pole 
and equator unchanged. For lat. 40°, for example, subt: 


racting the “sins” 
figures from the “‘tan z’’ ones and multiplying by Aé, 1”-000, we get 


Aj = +0°-428(cot ¢ — tan 2) — 0°-738(cos ¢ — sin 2) 

-0"-055 —0"-428 tan z + 0°°738 sin z, (43) 
and this expression can of course be multiplied by any other amplitude-constant 
instead without altering its property of leaving the pole-equator separation 
unchanged. (It vanishes at the pole, of course, but it does not and should not 
vanish at the equator; for this reason it might be a little difficult to construct it 
without using the above harmonic analysis.) By direct trial, one can see how 
large an amplitude it may be given without seriously affecting the below-pole 
minus above-pole declination differences, 4’—5, on which the adjustment of the 
pole itself depended ; if the original choice, as between a flexure adjustment and 
a refraction adjustment, was demonstrably wrong, the application of (43), with 
a suitable (perhaps negative) amplitude, will result in an improved run of these 
8’ —3, but even when this does not happen, the run will not be made worse over 


an appreciable range of amplitudes. ‘This is, of course, merely a quantitative 
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version of the qualitative statement that flexure and refraction are not clearly 
separable. (In general, applying (43) will also slightly change the pole itself, 
owing to the fact that a finite range of circumpolar declinations is involved , 
this effect can be dealt with by a further approximation if necessary.) 

The harmonic analysis of (43), taken directly from the appropriate coefhicient- 
differences in Table I, is 


Ad = — 0°:086 — 0":086 cos 2% + 0" 002 cos 44 — 0"'121 +0"-O18 sin 3%, (44) 


and the corresponding solar corrections, from (22), (24), and (25), are (for this 
amplitude, i.e. 1°00) 
Ave Ak@= +022, and (necessarily) A.4,=0; 
even though the largest amplitude likely to be permissible (so far as the circum- 
polars go) is probably not over 0-3 or o”-4, the uncertainty thus introduced is 
somewhat unwelcome. 
Actually, the above analysis is an over-simplification. In general, corrections 
to the solar clements are derived from residuals in « as well as from those in 4, 
and it is usual to form combined normal equations and treat both sets of residuals 
together. ‘The equation of condition for the right ascensions, if (as is usually 
done) we give equal weights to all Aa rather than to all Az cos, 1s 
E + cose — cos« tan dAc + 2 sin «sec SAA 
2cos« sec 6 coseAk = Az, 
and if we eliminate 6 as before this becomes 
+ (1004 — 0°086 cos 2a) AL — 0-217 sin 
+ (2°1368in « — 0°044 sin 3a +...)AA 
(1°877 cosa + 0°042 cos 3a +...)AR= Ax. (46) 
‘This equation shows much the same kind of orthogonality as (17); and when 
they are combined, with equal weights for equal intervals of « as before, almost 
all the non-diagonal terms in the normal equations still vanish. If we add (17) 
and (18) and solve them together with (46) in this way, remembering that (17) 
and (46) themselves have already been solved together, and giving equal weights 
to an observation in R.A. and to one in Dec. (as is usual), we find that equations 
(21), (22) and (23) reappear with their right-hand sides and the terms in A,d, 
unchanged ; the coefficient of A.« in (22) is increased from 0-500 to 0-524, and the 
coefficient of Ak in (23) is increased from 0-231 to 1-994. Equations (25), (22) 
and (24) are thus replaced by 
Abo = Cg 00420, O'0010, = + Ady (47) 
Ae = 09544, (48) 
A Ak = — 0° 108¢, — 0°002¢,. (49) 
‘Table | can readily be re-normalized by dividing each column by the value 
of ¢9~0°042¢,~—0°001¢, appropriate to that column. ‘The “¢,”’ which result 
are now all very close to 1-0, while the new “¢,”’ (and “‘c,’’) are not substantially 
altered; thus the rigorous criterion (47), for adjusting 4,5, which replaces (26), 
is now only slightly different from 
Co = Abo, (50) 
and falls almost half-way between (14) and (31). A,& is in all cases now very 
small, while A.« is changed less than 15 per cent from its previous value; these 
last two statements reflect the fact, well known qualitatively, that in a “ combined” 
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solution Ak is determined mainly by the «-residuals and Ae by those ind. Except 
where the original c, (Table 1) differ markedly from unity, such difference- 
equations as (43) and (44) reappear without any major modification, and the 
star corrections are also little affected. 

The above conclusions have been obtained on the basis of equal weights for 
normal points, both in « and 4, which are evaluated for equally-spaced «-values. 
Weights are sometimes assigned in proportion to the numbers of observations 
made in these ranges, and this will in general involve reduced weights for the 
winter months, and will somewhat upset the orthogonality and separability which 
have much facilitated the present analysis. 
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For any moderately uniform distri- 
bution of observations, the results will not be very greatly altered; but there are 
additional weighting considerations, some of them difficult, and the total amount 
of uncertainty which these import may be appreciable. 

In the first place, there is some evidence that residuals in « (seconds of arc) 
should be replaced by residuals in «cos; this would reduce the weights near 
both solstices systematically, but not by very large amounts. Secondly, the usual 
assumption that observations in « and 6 have equal weights is probably based 
on little more than the difficulty of coming to any decision at all on really objective 
grounds. ‘The systematic errors of the two methods are certainly quite different ; 
but they are very hard to assess. ‘The question how far (22), (24), and (25) should 
or should not be replaced by (47), (48), and (49) depends entirely on this relative 
weighting, and it is difficult to feel much confidence about any particular choice, 
The resulting uncertainty is not great except for Ak, but is clearly not reflected 
in the probable errors. In the third place, it would seem at first sight logical 
to underweight all low-altitude declination measurements, and this includes all 
those made on the Sun in winter, and not only those made on sub-pole stars ; 
at Greenwich, the winter Sun is only a few degrees higher than the limit at which 
sub-pole declinations are rejected altogether. ‘The temptation to give full 
weight to every measurement of the Sun’s declination which can be obtained in 
the winter months, even though the altitude is low, is a strong one, since for 
several months on end no observations except low-altitude ones can be obtained 
at all; but if the full weight is not really there it would be better to accept the fact. 
The weights at low altitudes can be adequately assessed for night-time observa- 
tions, both in the north and in the south, and it is improbable that the weight at 
noon could be better. Against this must be set, however, the possibility that 
noon observations in summer may be so much worse than those on high-altitude 
stars as to be no better than the winter ones; it is the scatter in the solar residuals 
themselves which should really decide their formal weights. ‘The same arguments 
apply in R.A. too, in principle, though the actual variation of weight may be less 
here. Whenever it seems even reasonably correct, taking everything into 
account, to assign equal weights to normal points that are equally spaced in «, 
the resulting simplification of the analysis is so great that it should probably be 
done, but in some cases it might clearly be improper. 

However this may be, there can be little question that so long as the real 
form of the flexure function is unknown it will be impossible to eliminate most of 
the uncertainties discussed above, either in the solar corrections Ak and Ae, or 
in the systematic declination corrections of the star places which may be obtained 
with the instrument in question; a detailed and reliable determination of the 


complete flexure function would clearly still be worth < considerable effort, if 
it could be achieved. 


2 
§ 
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The “ pentag’’ method of determining flexure.—I\t seems possible that this 
might be done, more completely than it has in the past, by making use of an 
optical square, or “ pentag’’. ‘This consists in principle of two plane mirrors 
very rigidly mounted together so as to enclose an angle of 45°; it has the property 
that a ray of light incident on one mirror, in the plane perpendicular to the line 
joining them, and then reflected in the other mirror also, is deviated through 
go” whatever its angle of incidence. If the angle between the mirrors is 45° + 0, 
the ray is deviated by go°~2v. For our present purpose, the exact value of v is 
quite unimportant (so long as it does not exceed 1‘ or so), but its constancy must 
be guaranteed within a small fraction of a second. Small “ squares’’, as used in 

ange-finders etc., are commonly made as solid prisms, and the most economical 

form is then pentagonal. When larger ones are required, fabrication from plates 
1s more suitable; there are obviously several ways of constructing a “hollow 
pentag’’ out of fused quartz only, and it should be quite practicable to achieve 
the necessary permanence even with a fairly large effective aperture. 

If such a pentag is mounted in front of the objective of a transit circle, by a 
counterpoised system supported only on the east-west axis, or on the ball-races 
through which the main counterpoise system acts, the circle reading for any star 
will always be go” — 2v greater, or less, than that for a direct observation of the 
same star, except for the difference of the two flexures. If c, is the reading with 
the pentag, and ¢ the reading without, and if the upper sign refers throughout to 
the arrangement in which c,~c +7 2, we have 


z=c+f(c) (51) 

Fm/2+20+f(e,) (52) 
and so 

(53) 


Writing f(c) = X(a, cos nc +b, sinnc), we obtain the equation of condition as 
f(e)—f(e+ + 2X[sin sin n(c + 4) 
b,, cos n(c + 7/4)}] 20—c,—cF 7/2; (54) 


or, putting in successive values of n and transforming so as to remove half the 
double signs 


\V 24, C08 (¢ 7/4) + 2a, C08 + 2a, 7 4)+... 


+ 2b, sin (¢ 7/4) + 2b, sin 2¢ + 2b, 8in (3c +7 /4)+...F2 


~ 


=¢,—¢cF n/2. (55 

The coefficients vanish identically for harmonics of order n = 4m, so that terms 
of these orders in the fiexuire function cannot be determined by this method. 
Although one may not feel entirely certain that a fourth harmonic will be negligible 
merely because the third and fifth ones may be found to be so, it appears somewhat 
improbable that any harmonic (of c) as high as the fourth can really be appreciable, 
unless it is used to describe a discontinuity. ‘The only likely discontinuity 
would seem to be in the zenith, and a discontinuity there is readily discoverable 
and can with advantage be explicitly removed before proceeding to any Fourier 
analysis; the indeterminacy of the fourth harmonic should then be unimportant. 
We must, however, consider the case n =o. 

A zero harmonic in the “ flexure" is the same thing as a constant error in all 
the circle readings used in flexure work or on stars; it may be put down as a 
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44! 
systematic error in taking the nadir, or as a “‘ wrong”’ division error of the nadir 
division. 


As is well known, this is eliminated by reversal, if it really is a wrong 
division error; if it is a systematic error in taking the nadir which is not reversed 
by reversing the instrument, it is in effect a w rong co-latitude. In either case, it 
seems substantially irrelevant, so far as the astronomical results are concerned. 
As there is no absolute term in the trigonometrical series in (55), the quantity ¢, 
the intrinsic error of manufacture in the pentag, can be obtained as one of the 
unknowns. Formally, it might always be convenient to do this; but it seems 
improbable that its value could be obtained in this way with as much accuracy 


as should be attainable in a pure laboratory determination. If not, and if the 


least-squares value always agreed with the laboratory one within their probable 
errors, one would of course gain by substituting the laboratory value, and 
transferring the term 20 in (55) to the right-hand side. 

The unknowns in (55) would evidently all separate completely if observations 
of stars could be made all round the circle. The total range of ¢,-values is much 
greater than the range of z-values of observable stars, since the pentag allows 
stars to be observed when the telescope is pointing very considerably below the 
horizon; but the total range of c-values is of course only the same as the range 


of observable z, say from — 80° to +80", and the range over which three different 


telescope positions are possible for the same star is appreciably more restricted 
again, since when the telescope is too far below the horizontal the incoming beam 
will foul a collimator. ‘There is, however, no star which cannot be observed in 
at least two positions, and the separation of the unknowns should be quite 
adequate, for a moderate number of terms. 

There seems no need to construct a pentag large enough to cover the entire 
aperture of the instrument. A four-inch beam will allow observations on a large 
enough number of stars for a flexure programme, and one would then presumably 
stop the beam down to this figure whether or not the pentag was in. ‘The various 
coefficients in (55) can be evaluated for each star in the programme once for all. 

If the indeterminacy of the fourth harmonic is considered undesirable, a 
“pentag”’ giving some other deflection than go could be constructed. ‘The 
modifications in the analysis are straightforward, but the laboratory determination 
of the exact angle might be difficult and the astronomical one might have to be 
accepted ; there are also some other disadvantages 

Although errors in positioning the pentag which affect the angle of incidence 
are completely immaterial, errors which tilt the plane of incidence are not. A small 
rotation from the correct position, either about the direction of the incoming beam 
or about the axis of the telescope, will mean that the star is observed somewhat 
before or after transit, and therefore not quite at culmination. ‘The error should 
however be casy to allow for, since the error in transit time (which will be known) 
is of the first order while that in zenith distance is of the second 

It thus appears that use of a pentag should permit a better determination of 
the flexure than has as yet been possible by any method. It i 
mercury pool in permitting a second harmonic to be detected; in permitting 
observations over a much greater range of z-values, including all values near the 
zenith itself; and also in being far less liable to be interfered with by a breeze. 
(Since the mercury pool demands extreme steadiness in the air surrounding it, 
it introduces a strong temptation to reduce the pavilion aperture dangerously 
at all times.) The pentag is superior to all “ pure” instrumental devices (even 


superior to the 
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if these were, in their own field, free from objection) in using stellar observations 
exclusively, so that no systematic error due to the rather different nature of a 
collimation, or auto-collimation, observation can enter; in addition, of course, it 
is superior in that it includes the effects of any circle flexure (which no pure 
instrumental method except the collimator one has yet succeeded in doing), 
and especially in being almost independent of the precision with which the 
auxiliary apparatus itself is mounted. 


Royal Observatory, 
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Summary 


A statistical model of the variation of latitude is formulated and a 
detailed analysis of its properties is made. Difficulties 
necessity to account for observational errors are discussed 

\ computational procedure is obtained and applied to a 35-year series 
of observations in which the ¢-term has not been taken into account. It is 
found that the damping factor is 


arising from the 


kK = 0°322 + 0'103 year™', 


giving a relaxation time of about 3:1 years; and the free period is 


T y~ 1°267 + 0°039 years 


7 
462°8+ 14°2 mean solar days 
The value of « is several times that of the only previous determination, and 
the free period is longer than has hitherto been determined 


These results are provisional pending the completion of a full numerical 
analysis 


Part |. Introduction 
1. If 4, m are direction cosines of the instantaneous axis of rotation of 
the Earth with respect to right-handed 


axes fixed in the Earth, the dynamical 
equations of motion may be written 


/ +ym= — A, (1) 
m . yl = A, 


the solution of which, if /,, m, are the 


initial values of / and m, is 
l= 1, cos yt — sin yt +1, (2) 
m = ly sinyt + my cos yt + my. 
Here, 27 y is the so-called Chandler period of the (circular) free component of 
the motion, w, is the mean angular velocity of the Earth, (A, A, C) are the moments 
of inertia and F and G are the products of inertia of the Earth. Assuming that 
the products of inertia are periodic with periodicities of a year and sub-multiples 
of a year, the forced motion, /,, m,, will have elliptical components corresponding 
to each period in F and G. 
This dynamical theory suggests that harmonic analysis of values of / and 
m, which are observed as the variation of latitude, should yield the value of the 
Chandler period and the various periods of the forced motion. However, very 
many investigations over different intervals have been made by a number of workers 
and have all yielded different values of the Chandler period. As early as 1917, 
Kimura gave an empirical formula for 27 y containing components with periods of 
22, 26°67, 50 and 8o years from which Melchior (1) has calculated a table of values 


of the Chandler period during the interval 1g00-1960. The maximum value 
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is 1°2066 years in 1914 and the minimum is 1°1123 years in 1929. Dhifferent 
authors have found, by harmonic analysis of different intervals, Chandler periods 
within this range of values. Melchior has developed a theory relating variations 
of the Chandler period to changes in the axes of the elliptic forced motion. 

In 1940, Harold Jeffreys (2) rejected the validity of harmonic analysis as 
a means of determining the Chandler period on the ground that undamped free 
vibrations in the sense of elementary dynamics probably have no physical 
existence. Irregular disturbances of the motion may well repeatedly change the 
phase of the motion; despite this there exists a natural free period. 

The appearance of the series of observations suggests that statistical fluctua- 
tions are not negligible, and harmonic analysis cannot take such fluctuations 
adequately into account. Jeffreys’s paper, in fact, contains the only example so 
far published of an analysis which ts based on a proper statistical model. _Jeffreys’s 
treatment is re-examined in this paper in the light of recent developments in the 
theory of time series analysis. ‘lhe problems of estimation and statistical inference 
associated with his model are discussed in detail, account being taken of the fact 
that the effect of observational errors may not be negligible. ‘The theory of the 
statistical analysis is given below in Part II]. ‘The principal steps of the com- 
putational procedure based on the statistical analysis are outlined at the beginning 
of Part [1], which contains the results of carrying out this analysis on a certain 
scries of observations. 


An analysis of the residual motion is given in Part LV. 


Part Il. The Statistical Model 
2. Specification of the statistical model.—The equations of motion (1), 
including a damping term, may be written 


ym \ (3) 


where Fy = a,k A, G,=aegG A and « is the coefficient of damping. 
Now let 
Py(t)= Flt) (4) 
G(t) = Gt) + 
where /’y, Gg are the systematic parts of F', and G, respectively and /,(t), g,(¢) are 
the fluctuation terms, | /,(7), g,(t)} being in fact a stochastic process in two variables 
with zero means. We make the usual assumption that F's and Gg have a yearly 


period, and can be represented to a sufficiently good approximation by the first 
few terms of their Fourier series, e.g. 


X cos vrt + by, vrt), 


where p is a small integer and 27/v=1 year. 
The simplest statistical assumption as regards the fluctuation terms is that 
i/(t)} and {g,(t)} are independent stationary impulse processes. ‘This means 


that only the integrals /(t)= | f/,(u)du and J(t)=| g,(u)du strictly exist, {J(t)} 


and |J(t)} being statistically independent homogeneous additive processes, 
characterized by the property that the change of J (or /) in any time interval 
depends only on the length of the interval, and changes in non-overlapping 


by 
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tume intervals are statistically independent of one another (see e.g. Bartlett (3), 
Chapter III, Section 6, and Moyal (4), page 165). If the impulses are dis- 
crete, the changes in / or are discontinuous 

rom (3), #yand Gy give rise to forced components /,(7) and m,(t) of the same 
torm (the damped components may be ignored since the motion may be assumed 
to have started a very long time ago); and for the fluctuation terms we have 


x(t) + ry(t) “(dd (u) — dI(u)), (5) 
where 
x(t) = l(t) — L(t), 

) ) (6) 
y(t) = m(t)— m,(t). 


Thus x, y are the contributions to the motion of the statistical fluctuations in the 
products of inertia of the Earth. 


If we now consider observations at times separated by a constant interval A, 
and write 


Utyp+nh)=1 +x 


m, = m(t, +nh)=m,, +4 
we tind that 


where 
a=e€ cos yh, 
B=e “sin yh, 
“h 
e™\cos yudI(n+ th—u)+sinyu dd (n+ th—u)}, (8) 
and 
h 
Bess | e™ yu dJ(n+ th—u)—sinyu dI(n + th—u)}. 
6 


‘The residuals {e,, 7,} form an independent stationary process, i.e. its properties 
are invariant under translation of the time axis, and {e,,,,,,,}, (70, 1, 2...) 
are all statistically independent. If we make the further natural assumption 


that the rates of increase of variance per unit time for / and J are equal, we have 
also 


vare,, = vary, 
and (9) 
cov (« ar In) =O. 
Thus we have arrived at a model identical with that postulated by Jeffreys 
(2, page 143) apart from his additional assumption that the residuals follow a 
normal distribution. 


Our assumptions for the fluctuation terms can of course give only a very 


approximate representation of the actual mechanism, but an analysis based on 
them may nevertheless still be sufficiently accurate to be useful. ‘That Jeffreys’s 


odel follows from our statistical treatment may be regarded as an argument 
for using this as a working hypothesis. In (2) the model 1s introduced by analogy 
with Yule’s treatment of the problem of a pendulum bombarded by randomly 
arriving particles, but this has since been shown to be incorrect (see e.g. Bartlett 


(5), page 35). 
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The assumption of normality of the residuals is more difficult to justify. 
‘This is approximately true only provided that the mean number of impulses 
during the interval of observation, /, is large, a condition which is perhaps not 
very realistic for values of A which are usually considered in practice (of the order 
of a month). Jeffreys even claims (2, page 148) to have identified a period of 
more than 12 years during which the effects of fluctuations were negligible. 
Hence it is advisable to dispense with the normality assumption if at all possible. 
It is not difficult to show that we can do so provided that the effect of observational 
error can be ignored, so that equations (7) can be used as they stand. However, 
as Jeffreys points out, observational errors are not necessarily small compared 
with x, and y 


,»» 80 that they should be allowed for in the model. Unfortunately 
when this is done an important part of the analysis becomes very complicated 
unless the residuals are normally distributed, and indeed is tractable only when 
both the residuals and observational errors are normally distributed. 

We therefore re-define (x,, y,,) to include the effect of errors of observation and 
assure that 


(10) 


where both {x,’, y,’}, which satisfies (7), and the errors of observation, {u 
are normal. It seems reasonable to assume that 


nie 


varu, =Varv, say, 

(11) 
cov (u,, U,)=0, ) 

and (0, 1, 2,...), are all statistically independent and also that 

|v,» U,} 18 statistically independent of the series {x’,, y’,}. 

These properties, together with 


(a, cos vrnh + b,, sin vrnh), 


(a,,,cosvrnh sin vrnh), 


au, + BU, 44, SAY, 


(13) 


= t — BU, —2U, nay Bay, J 


completely specify the model. Equation (12) specifies the true dynamical 


contributions to the motion as opposed to the statistical fluctuations ; equation (13) 
results from inserting (10) in (7). 

3. Correlational and spectral properties of the model, —-\t is not difficult to show 
from (13) that the autocovariance matrix for lag s, defined by 


are 
V(s)= , 


ate cossyh sinsyh 
V(s)=: mal (s>o), | 


sinsyh cossyh 


(14) 


is given by 


+ 
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‘The autocorrelation matrix for lag s, (s 0), obtained by replacing x,, y, in V(s) 
by x,,/o(x), y,/o(y), where o*%(x) = var x, o*( y) = var y, is therefore 


cossyh sin 
R(s)=Ae 

sinsyh cos syh (15) 
where 


A=a? {a*+a'%1—e 


(16) 
Equation (15) summarizes the main correlational properties of the model and was 
in effect given by Jeffreys. ‘The theoretical correlograms for the autocorrelations 
of the |x} and {y} series are identical, consisting of an exponentially damped 
simple harmonic oscillation, and a phase shift of go” yields the theoretical 
correlogram for the cross-correlations. 

This result expresses the physical fact that, apart from the damping, the 
expected motion is given by a simple rotation of the pole about its mean position 
through the angle syh in time sh. ‘The effect of the observational errors is seen 
in the reduction of all the lag correlations by the factor A. . 

We also note that the spectral density matrix F(w)=(1 27) L Vise”, 


i.e. the Fourier transform of the autocorrelation matrix, is given by 
F(w) =a" +a? + 24 cosw)* — 48" sin? w} ! 


+a? + B*— 29 cosw 218 sinw 


(17) 
216 sin w + — 24 cos uw, 


4. Analysis of the systematic components of the variation.-We now have to 
consider methods of estimating the unknown parameters in the model, and of 
calculating the precision of our estimates, given an actual series of observations 
\L., m,}(r=1,2,...N). There is also the associated problem of testing whether 
the model is consistent with the observations. ‘lhe parameters fall naturally 
into two groups, the first consisting of those which occur in the equations (12), 
Dips Ome, Aefining the systematic components of the variation, 
and the second of those associated with the equations (8) and (13), ie. a, B 
(or « and y), o* and o*, which define the behaviour of the fluctuations. ‘The two 
groups can be dealt with to a large extent independently of each other; we must 
however begin by considering estimates of the parameters in the first group, 
since these are required before we can calculate the estimates for the second group. 
We shall assume from now on that N is sufficiently large for an asymptotic theory 
to be used; this may not necessarily be true in practice (see Section 6 below), 
but no exact theory has yet been developed. 

It can be shown by using an approach due to Whittle (6, 7) that the maximum 
likelihood estimates are asymptotically the same as those obtained by standard 
harmonic analysis i.e. 


2 24 3 2 
c*, /,cos vrth, b°,, sin vrih, 
N inl \ i—1 


They are asymptotically unbiased (and are in fact exactly unbiased if the 
observations cover a number of complete years). ‘Their asymptotic variances 
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and covariances are given by the inverse of the information matrix as in classical 
maximum likelihood theory, but are easily obtained directly. We find that the 
estimates fall into asymptotically uncorrelated groups (c*,, c*,,), (4%), ny, 
b* na)» (2% may 49), etc., and that if we write =(c*, c*,,), = 
(4* 1, b*,’ =(b%,, b*,,,), etc., these being row vectors and c*, a*,, b*,, etc. 
the corresponding column vectors, 

cov(c*, 

cov (a*,, a*,’)~cov (b*,, (19) 

cov (a*,, b*,’)~ — 24 
where F(w) is the spectral density matrix given by (17), and the symbols # and 4 
denote real and imaginary parts. Thus c*,, c*,, are asymptotically uncorrelated 


and have variance 
2 
ao 
a*-+ 
1+a*+ =t/ 


and the asymptotic covariance matrix for the group (a*,,, a*,,,, 6*,,, 6*,,,) is 


mr? ir mr 
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a* ny b ir b mr 


Y, 


b° Vv, P 


r 
where a*(1 + a* + B® — 20 cos vrh) 

N | (1 +a" + — 2a cos vrh)* sin® vrh | 


and O 4 o*8 sin vrh f 
|(1 + — 2a cos vrh)* — 467 sin* 


(21) 


af 


All these estimates are also distributed asymptotically in the multivariate 
normal form, so that we can not only calculate their estimated standard errors, 
once estimates of «, 8, o*, a’*, with standard errors which tend to zero as N >, 
have been obtained (estimates which have standard errors proportional to 1/\//N 
are given in Section 5), but can also derive approximate tests of significance and 
confidence intervals for the true values of the parameters. For example, a*, 
would be regarded as differing from zero at a level of probability approximately 
equal to o-gs if it was greater in magnitude than 1-964/ P*,, where P*, denotes the 
estimated value of P,. Again the significance of the whole set of estimates 
6*,,,) Of the (r —1)th harmonic components could be judged by 
referring the expression 

{P* (a*,,? + 20° ng ‘ (P¢,3 = O*,*) 
to a y® distribution with four degrees of freedom; this follows from the fact that 
the inverse of the matrix (20) is 


I 0, o 
0, 
0, 


‘The latter test is useful in the problem of deciding on the appropriate value of 
p, i.c. on the order of the highest harmonics which should be included in the 
systematic variation. ‘his cannot be specified exactly theoretically although 
it may be assumed that the amplitudes of the harmonic terms will, on the whole, 


d 
| 
|, 


No. 4, 1955 The variation of latitude 449 


decrease as their order increases. ‘Thus one might proceed by first allowing for 
fundamental terms only, i.e. taking p= 1, and testing whether their amplitudes 
differ significantly from zero; then if they do so (which should certainly be the 
case here) taking p= 2 and testing whether the amplitudes of the first harmonic 
terms differ significantly from zero, and so on, the value of p finally adopted 
corresponding to the order of the last harmonic terms whose amplitudes were 
significant. ‘This method is somewhat laborious, since at each stage the estimates 
a*, B*, 0°", o** occurring in P*, and O*, have to be recalculated, and will of course 
be misleading if some high-order harmonic has a fairly large amplitude; but it 
becomes more satisfactory if we increase the value of p by more than unity at 
each stage, and also continue as a check for at least one stage beyond that which 
first yields no significant harmonic terms. However, if the interval between 
consecutive observations is a rational fraction of the period (i.e. a year), there 
will be a maximum number of distinguishable harmonic terms irrespective of 
the number of observations, and if this number is not too large, it will be best to 
include them all in the specification (12) and carry out a complete harmonic 
analysis, even though some terms are found subsequently not to be significantly 
different from zero (see Section 9g). 

It may be noted that the above analysis can be shown to be asymptotically 
valid even when the residual fluctuations {e,, ,} and the observational errors 
[Un U,} are not normally distributed. It is sufficient that the statistical 
independence properties stated in Section 2 should hold. ‘The estimates (18) 
are of course then no longer maximum likelihood estimates. 

There remains the question as to whether equations of the type (12) give 
a satisfactory representation of the systematic variation. Here, however, there 
is the difficulty that to test whether the observations are consistent with the model 
we must in general postulate some alternative hypothesis for the form of the 
systematic variation: for example, that we require additional terms dm and dim 
in /,,, and m,,, respectively to take account of a linear trend (corresponding perhaps 
to variation with a very long period). ‘The only practicable way of avoiding this 
difficulty seems to be to divide the series into several consecutive sub-series, 
from each of which estimates of the form (18) are obtained. ‘Then if the sub-series 
are sufficiently long for the asymptotic formulae (19) to apply, and for the correla 
tions between the estimates for the different sub-series to be neglected, it is fairly 
easy to test whether corresponding estimates differ significantly. It is of cours 
useful to examine the graphs of the estimated / and m residuals 


{a*,, cos vrnh + b*,, sin vrnh}, | 
rel | (22) 

m,* =m, —c*,,— {a*,,,cosvrnh + sin vrnh}, 

rel 
plotted against n, in case these reveal any unusual features (but see Section 12) 
5. Analysis of the fluctuation components._-Estimates of the second group of 
parameters a, £, o* and o” are obtained from the auto-covariances and cross 
covariances of the series {/,*, m,,*}, which are defined by 


n 


1 
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The effect of the differences between /,* and x,, and between m,* and y, will 

become negligible for sufficiently large \, so that we can regard c,,(r) as equivalent 

to I 


and so on. 


If we do not take account of observational error, i.c. assume that o'* =o, 
the procedure is straightforward. Maximum likelihood estimates of « and £ 


are asymptotically equivalent to those obtained by the method of least squares, 
i.e. by minimizing 


—al,,* + Bm, * +(m,..,* — Bl,* —am,*))}, 


which gives 


f vis (24) 
- Also o** is obtained as the mean square residual, equivalent to 
= + — — B**}. (25) 


It is easily shown that «* and 8* are asymptotically unbiased and that 
var a*~var — a? — (26) 
cov(«*, B*)~o,whence with estimates of « and y obtained from 


these being also asymptotically unbiased, we find that 


var (x*h)~var (y*h)~{e™ — 1} 2N. (28) 


although the estimates are not then equivalent to those obtained by the method of 
maximum likelihood. 


‘These formulae remain valid when the «, and 7, are not normally distributed, 


In the general case, however, the maximum likelihood approach becomes very 
complicated. It is possible to proceed by the method of Whittle (7) but this 
involves very lengthy algebraic calculations, and yields a set of equations which 
are difficult to solve. A more practical method consists of using covariances for 
lags h and rh, as suggested by Jeffreys. ‘The above estimate of y, given by 


cot y*h = + 1)} = say, (29) 


where 


P(r) = + Coal”), = Cyl) (30) 
is easily seen (remembering equation (15)) to be consistent, i.e. to differ from y 
by not more than any pre-assigned quantity «, however small, with probability 
arbitrarily near unity, if N is taken to be sufficiently large. Similarly 
| p(r)}® + {q(r)}* is found to be a consistent estimate of and so 
a consistent estimate of « is given by 


MD p(1)}* + + (gr) (31) 
‘The best value of 7 to take is not obvious, Clearly for large r the precision of «* 
will be low since the expected value of the denominator in (31) tends to zero 
exponentially with r, but there is also the effect of the final division by r—1. 
Jeffreys states that r should be chosen so that e¢ mh is about 0°45, but does not 
explain how he arrives at this result, and it is incorrect at least in the case of small 
xh, for which it can be shown that when A~1, i.e. the observational error variance 


: \ 
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In view of this, and the 
fact that the precision of «* can be calculated most accurately for small r, it seems 
best in general to use (31) with r= 2 


is small, r should be as small as possible (see Section 6). 


It is easy to see that y* and «* are asymptotically unbiased, and that their 
asymptotic standard errors are proportional to 1 /N. 


Formulae for these may be 
obtained by expressing y* — y,«* 


y « a8 asymptotic linear functions of the quantities 
Op(t), Og(1), Op(r), dg(r) (where 5p(1) = p(1)—p(1) and p(1) is the mathematical 
expectation of p(1), etc.), and then using the expressions for the asymptotic 
variances and covariances of the p’s and q’s. ‘The latter follow from formulae 
of the type 

I 


N +5) + + 27 +5)}, (32) 
where e.g. 
cov {X,(t), Xg(t + = (33) 


(In this paper the sign ~ means “ is asymptotically equal to ”’.) 
For instance we find that 


(N —r) var p(r)~2{o" + o7/(1 ™)}* 
x | cos aryh{ar + 1 
+ cos |. 


After some rather heavy algebra, we arrive at the results 


oa I A 
I 4 lmh 2+ (2r 
var {«*h(r —1)}~ oN | 
(35) 
In particular for r= 2, (35) gives 
var kth = — 1 4+ 2(1 (36) 
‘To estimate A, we may use the equation 
Aste p(1)}* + + (37) 


From (34) and (35) we obtain estimated standard errors of «* and y* on sub- 
stituting «* and A* for« and A, and hence approximate confidence limits for « and y, 
since «* and y* can be shown to be asymptotically normally distributed. 
asymptotically normally distributed, with mean A and variance 


var — 1)* + 2A(1 — 1) + (1 — + + 2))/2N. (38) 


A* is also 


We can therefore assess approximately the significance of the difference of A* from 
unity. If A* was considered not to be significantly lower than unity at an appro- 
priate probability level, we should regard the observations as consistent with the 


4 
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hypothesis that observational errors can be neglected. Finally a consistent 
estimate of o*, which is required in the estimation of the variances and covariances 
(19) of the harmonic components, may be obtained from the equation 


+ = + 099/(1 = — (39) 


It should be noted that the normality assumptions of Section 2 are really 
used only in the derivation of equations (34), (35) and (38). Without these 
assumptions it is in general very troublesome to derive expressions for the 
asymptotic variances of y*, «* and A*, and the estimation of the variances is no 
longer straightforward. 

The validity of the correlational structure given by the model (13) may be 
examined by comparing the sample covariances ¢,(r) with the theoretical 
covariances V,,(r) defined by (14), or the sample correlations (obtained by dividing 
the corresponding covariances by +/ {¢,0)c;(0)}) with the theoretical correlations 
defined by (15). We can, for example, use the goodness-of-fit test associated 
with the ‘‘ G-operators” introduced by Bartlett and Rajalakshman (8); this can 
be shown to be applicable here although the model is more general than those 
considered in their paper. In addition we can test whether c,,(0)/c..(0) and 
¢,,(0), which are asymptotically uncorrelated with the ‘ G-forms” and with each 
other, differ significantly from 1 and o respectively, using the expressions for their 
asymptotic standard errors obtained by means of the general formula (32). 
It seems, however, difficult to construct a suitable test of the assumption that the 
residuals are normally distributed. Some indication of the appropriateness 
of the assumption may be obtained by constructing frequency distributions of 
Liar” + Bem,* and m,,,*—B*l,* —a*m,,*, but it is not known how the 
significance of deviations of these from the normal form should be assessed ; 
the standard y? distribution certainly cannot be used. 

6. The effect of a small damping factor.—Difficulties may occur in the above 
analysis if xh, the product of the damping factor and the interval between 
consecutive observations, is small. We then have from (35), provided that rxh 
is also small, 


var {x*h(r — 1)}= [A*%(r — + (1 + (1 —A)eh{2A(r — 1) + (1 + 1)}], 


(40) 
so that if V(«*) denotes the coefficient of variation of x*, 
(1 ~A)? + (1 —A)wh{2A(r 1) +(1 + 1)} 


If rxh is not necessarily small, the first term on the right-hand side of (41) 
is replaced by (c?”—1—2y)/2N«hy*, where y=(r—1)ch. Since this increases 
with y, V(«*) will be approximately minimized by taking r to be as small as possible, 
i.e. r=2, when 1—A is sufficiently small for the remaining terms to be 
neglected. If the latter are taken into account, the choice of y to minimize V(«*) 
depends on « and A in a rather complicated way; however, as was indicated in 
Section 5, the accuracy of the approximation (35) tends to decrease as r increases, 
and the value r = 2 seems suitable in general. 

Since, to a good approximation, V(«*) >(N«h)-** it is desirable that Nxh 
should be appreciably greater than unity. Indeed, when this is not so, the 
approximate lower confidence limit for x, given e.g. by «*{1 — 1-g6V'*(«*)} for a 
confidence probability of o-95, where V’*(«*) is the estimated coefficient of 
variation, is likely to be negative, no matter how small the observational errors 


the 

Nive 
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may be. Thus, since negative values of « are inadmissible, we can only consider 
the upper confidence limit; the precision of «*, or the estimate 1/«* of the 
“ relaxation time”’, which has to a good approximation the same coefficient of 
variation V(«*), cannot then be regarded as satisfactory. 

There is, however, a more serious difficulty when the condition Nxh> 1 is not 
satisfied. ‘This is that the asymptotic formulae for var («*h) and var (y*h), the 
latter being approximately (1/2NA*)[2«h +(1 —A)(1 + 3A)] for small «hk, may no 
longer be good approximations. For the representation of y*—y and x®—« as 
linear functions of the 45p(r) and g(r), used in their derivation, is valid only on 
the assumption that the coefficients of variation of p(r) and g(r) are small, and it can 
be shown that (provided A is not small) these coefficients of variation are both 
approximately equal to(N«h)**. (This assumes that the asymptotic formulae for 
var p(r) and var g(r) are still applicable ; it is difficult to estimate the relative error 
of these, but a rough argument suggests that they are of the same order of magni- 

N 


tude as that caused by replacing Le by Le ic. as which is 
v=0 

e.g. less than o-o5 when Nexh>1°5.) ‘The same difficulty arises in connection 

with the formula (38) for varA*. A further point is that the distributions of 

«*h, y*h and A* may then no longer be approximately normal. 

These difficulties are really due to the fact that when «A is sufficiently small 
the elements of the autocorrelation matrix decrease very slowly as the lag increases. 
In particular, when the effect of observational errors is small, so that we have A~1, 
this can be shown to imply that the observations are equivalent to a very small 
number of independent observations, and hence we should not expect the 
asymptotic theory to apply. 

A more refined method of statistical analysis which avoids these difficulties 
does not seem to be available at present; in any case it is conjectured that when 
Nx*h is of order unity, i.e. the time Nh covered by the observations is comparable 
with the estimated “ relaxation time’, it is not possible to obtain much information 
about the true values of the parameters in the model. 

Jeffreys (2) estimated « to be about 0-06, the unit of time being a year. Part 
of the evidence for this estimate is provided by a statistical analysis of a series of 
154 observations at intervals of 0-3 year. However, with N=154, h=0-3, 
x= 0°06, Nich is only 2°8, so that the above difficulties certainly arise. In fact if 
we substitute these values of N, x, A together with the value A=0-933 given by 
Jeffreys (2, p. 145) in (41) we obtain V(«*)=0-72, which is sufficiently large to 
indicate that the analysis is not satisfactory _If Jeffreys’s estimate were accurate 
we should have the same situation arising in the analysis of the series of observations 
used in Part II] below, for which N= 420, h= 4. However, our estimate of « 
turns out to be much larger, i.e. «* = 0-322, giving N«*h = 11-3, and an estimated 
coefficient of variation V*(«*)=0°31. This value of N«*h is not particularly large, 
and some doubt may therefore be felt as to the accuracy of the asymptotic formulae 
for standard errors and confidence limits. Despite this, however, the full 
analysis, the results of which are summarized, was thought to be justified. 


ve 0 


Part I11. Computational Procedure and Preliminary Results 

7. Computational procedure.—Estimates c*,, €*,,,, 4% etc., defined by (18), 
are first obtained, it being supposed that the question of the order of the highest 
harmonic to be included has been settled. 
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These estimates are used to calculate the estimated residual series {/,*} and 
{m,*}| from (22). From these series, auto-covariances and cross-covariances for 
lags s = 0, 1 and 2 are calculated and substituted in formulae (29), (31), (37) and (39) 
of Section § to give estimates of x, y, 0? and o*. Further auto-covariances and 
cross-variances for higher lags are computed at this stage if the goodness-of-fit 
test mentioned at the end of Section 5 is to be applied to the results. 

The estimated standard errors of «*, y* and A* are then calculated from (34), 
(35) and (38), and finally the estimated standard errors of c*,, c*,,, etc. from (19), 
which also gives the covariances required for carrying out the approximate test 
of the significance of a set of harmonic coefficients of given order described in 
Section 4. 

8. Data used in the analysis.—As yet we have only applied the results of this 
analysis to a series of 420 ‘“‘monthly’’ values of / and m covering the 35 
years 1900-1935 taken from the various publications of the International Variation 
of Latitude Service. In some cases data have been readjusted after original 
publication ; the latest available values have always been taken in such cases. 

For the reason discussed in Section 6 the longest possible series of results 
are needed. (Unfortunately the data available have not been obtained in a uniform 
manner since 1900. Firstly, the number of stations at which observations have 
been made declined throughout the period. From the inception of the 
programme, at the beginning of 1go0, to 1915 January, six stations (Carloforte, 
T'schardjui or Kitab, Mizusawa, Ukiah, Cincinnati and Gaithersburg) were 
in operation. In 1915 January, Gaithersburg ceased to operate, as did Cincinnati 
in 1916 January and Kitab in 1919 May. From 191g May onwards only three 
stations were used in the determinations, this being the minimum number 
necessary for determining /, m and the so-called z-term. Secondly, the method 
of reducing the observations was changed at 1922-7. 

Observations are made on twelve star fields, consecutive pairs of which 
overlap by about a month, and the values of / and m pertain to the mean date of 
each group-combination. For example, in the early years, star Group | was used 
for observation from September 23 to December 6, Group II from November 2 
to January 4 and the published values of / and m for the group combination 
I-I1 refer to the mean date of the overlapping period November 2 to December 6, 
This gives 12 irregular intervals per year. After 1922-7 the group combinations 
were altered to make the intervals equal. A good account of the actual procedures 
and programmes is given by Melchior (10, Chapter 1). ‘The irregularities in 
the intervals may be regarded as leading to slight, irregular phase changes in the 
motion and their effect will be to increase the uncertainty of the results. We feel 
that these ser «' results are better than the 4; yearly series which are usually 
quoted and which j -ffreys used in his analysis. ‘The reason for this preference is 
that the most essential requirement of the series examined is that successive 
terms should be in no way correlated except by the physical processes underlying 
the motion. ‘This is true of the 4 yearly series but the 4 yearly series are 
derived from the former by graphical interpolation and in the process of drawing 
smooth curves through the rather irregular looking curves of the +] yearly 
series, correlations between successive runs of 3 or even 4 values must be 
introduced. Comparison of the graphs of the 4 and } yearly series do in 
fact indicate a high degree of smoothing. ‘here must therefore be bias introduced 
into the auto-covariances for small lags. In the fuller analysis which we hope 
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to make later this will be examined more closely, particularly as it has a bearing 
also on a point raised by Professor Jeffreys (private communication) regarding 
the effect of errors introduced by changes in the star fields occasioned by their 
proper motions, and the effects of precession. 

We have used the results from which the so-called z-term has not been removed, 
as the results so obtained appear to be possibly more homogeneous than those 
obtained when the z-term is taken into account. They are certainly unsmoothed. 
This point will require examination in our fuller analysis, as will also the effect of 
the important change in the method of reduction made in 1922-7. 

g. Preliminary results.—VYor ease of computing, the 35-year means for each 
month have been calculated and analysed by the usual simple method of numerical 
Fourier analysis. ‘This yields all terms up to the cosine term for the fifth harmonic. 


The results are given in Table I. As expected, the fundamental is the dominant 
term. 


Tasie I 


Results of harmonic analysis 
(Unit 0”-o1) 


c* a*, b*, a*s be, 
11062 §5°105 70'286 2°138 | 6-914 
m 10°057 74°790 45°276 5'545 3°105 


Estimated 
standard 


error 


2°345 4°O14 1°757 2°086 o'771 

m 2°235 + 3°943 1°024 1 ‘000 1°257 1°§62 
Estimated 

standard t+ 1°70 + 1°48 1°39 + 0°97 


error 


The series {/,*} and {m,,*} were easily obtained by deducting from each value 
of / and m the 35-year mean for the appropriate month. ‘This of course implies 
that all the components given in ‘Table | have been removed whether or not they 
later turn out to be significant. 

The auto-covariances and cross-covariances, ¢,(r), 7 = 0, 1, 2, calculated from 
the series {/,*} and {m,*} are given in ‘Table Il (taking o”-o1 as the unit of 
measurement). 

Taste Il 


Covartances of and 


Auto-covariances Cross-covariances 
(0) It 289°114 3 4 §33°685 
9 972'436 € 192°782 3 
€,;(2) 7 8 


12 377°847 6 4 480°670 6 
Caaf) 10 716°248 2 8 071636 4 
7 718°382 8 
From these were obtained the estimates and standard errors given in Table III. 
Finally we calculated the estimated errors of the harmonic coefficients. 
These are given in the bottom line of Table 1. ‘The first and second harmonic 
terms, though relatively small, clearly differ significantly from zero, but it appears 
that the higher-order terms could be omitted. ‘This is confirmed by applying 


+ 3°05 t 10°77 + 2°40 
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the x’ test described in Section 4, the results of which are given in Table IV. 
Ideally, the numerical analysis should now be repeated omitting the insignificant 
harmonics, but as their contributions to {/*} and {m*} are trivial this has not 
been considered necessary. 
Taste 
Constants for free motion 
Estimated 
standard error 


Tasie IV 
Significance of harmonic constants: y* test 
Order of harmonic I 2 3 4 


Value of y* 10°95 12°30 3°54 3°94 
Degrees of freedom 4 4 4 4 


PO?) <0'05 <0’02 >0'30 >0'20 
10. Free period and damping coefficient.—From Table III it is seen that 
(i) the estimated damping coefficient is 


K* = 0°322 + 0°103 year 
giving a relaxation time of 3:1 years ; 
(ii) the estimated free period is 


2n/y* = + year 
462°8 + 14°2 mean solar days. 


The period is better described as a 15-monthly rather than 14-monthly period. 

11. Discussion of the results. —We have already remarked in Section 6 that the 
estimate of the damping coefficient « is much larger than that given by Jeffreys. 

‘The estimated free period is greater than those obtained by previous analyses, 
although it is in fact still compatible with the traditional value 27/y= 1-2 year. 
From the analysis of data for the years 1892-1933 Jeffreys obtained 
2m y = 1°223 + 001g years, which is in quite good agreement with the present 
result. He preferred a result, however, based on the analysis of data for the 
much shorter period 1908-3-1921°5 from which he obtained 27/y = 1-202 + 
year. 

The need to take account of observational error in the model is clearly shown 
by the value of 1 —A*, 00197, which is more than 5-6 times its estimated standard 
error and certainly differs significantly from zero. When the analysis of Section 5 
ignoring observational error is carried out, we obtain the value «* = 0-560 for the 
damping. 

There appears to be only one other occasion on which so long a period has 
been suggested. Jeffreys (9, pp. 213-4) has made a solution, using Takeuchi’s 
theory of the shell and an approximate treatment of the core of the Earth. It 
predicts a period of 460 days for the free motion—in close agreement with our 
present determination. We understand, however, that this theory gives an 
unsatisfactory value for the nutation and that further work is still in progress. 

The standard deviation of the observational errors seems to be quite large, 
our estimate, 15-3 units, being more than four times that suggested by Jeffreys. 


y*h=23° 40 44 
A® =0°9803 + 0°003447 
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This may be due in part to the uncertainties being increased by the slight irregu- 
larities of the intervals as suggested in Section 8. The increased uncertainty 
may be due however to the non-removal of the z-term. 

In Sections 4 and 5 we have suggested other points which bear examination. 
We hope to report on these later. ‘They deal mainly with the adequacy of the 
model but we may mention that the application of the Bartlett-Rajalakshman 
goodness-of-fit test requires a very much larger amount of numerical computation 
than has even been involved in the analysis so far completed ! 


Part IV. The Residual Motion 


12. The distribution of the residuals,—Reverting to equation (13), it is seen 
that when «* and £* have been evaluated, estimates «'*, 1'* of «' and 7’ can be 
obtained by replacing x and y by their estimates /* and m*. We do not see how 
to separate the fluctuation components (¢, 7) from the observational error (u, v) 
and so cannot test them separately for normality and randomness. However, 
we have computed the series {e'*} and {7'*} and examined them. ‘They are both 
symmetrically distributed with zero means and have standard deviations which 
do not differ appreciably from those calculated from the values of the estimated 
parameters. Both distributions are flattened compared with the normal 
distribution, but not much so, ‘The residuals are displayed in Fig. 1. 

We conclude therefore that the evidence, as far as it goes, gives no reason for 
disbelieving the hypotheses of Section 2, namely, that both {«, 7} and {u, v} are 
normally distributed and form independent stationary processes. 

Jeffreys (2, p. 147) inspected the series corresponding to our {/*} and {m*} 
and concluded that there were abrupt changes in amplitude and phase at certain 
times. He chose the period 1908-3~1921°5 for further analysis since he found 
it to be one in which the motion was almost undisturbed. Inspection of Fig. 1 
shows similar phenomena, the period up to about 1915 being apparently 
disturbance free. ‘Thereafter the {/* } and {m*} series seem increasingly disturbed. 
To what extent this is due to physical causes, or to the smaller number of stations 
employed, or to the change in the method of reduction made in 1922 is not clear. 
The point is, however, not that the {/*} and {m*} series should follow a certain 
pattern but that the mean number of “impulses” during the intervals 4 should 
be large (Section 2) and this cannot be decided by a rather subjective examination 
of the actual series. ‘The real test of departure from the assumptions made in 
the model rests in the analysis of the residuals {«, 7}. Although these cannot be 
isolated, examination of the series {«'*} and {»'*}, as stated above, does not reveal 
any non-random effects. ‘There are nowhere any systematic runs of values longer 
than one would expect in random series. For this reason we do not regard as 
satisfactory any procedure for isolating periods for preferential analysis solely 
on the basis of the appearance of the series {/*} and {m*}. 
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ERRATUM 
M.N., 115, No. 1, 1955: 


D. G. Ewart, On the relation between the stream and the ellipsoid constants. 
P. 50, line 3 from bottom, for y=7,t read y=ar,t. 
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